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Abstract 


The analysis has shown that if the stratification is unstable for ascending motion over an 
infinite area and stable for descending motion, and small random perturbations of various 
scales are introduced, the final disturbance evolved will have the dimension of a cumulus 
cloud. The stable stratification in the descending region has the effect of increasing the critical 
lapse-rate and narrowing the ascending region and making the descending motion widespread, 
so that the centers of the ascending regions are further apart. The distance estimated from the 
theoretical analysis agrees favorably with those obtained from satellite observations. 

The distributions of the various perturbation quantities pertaining to a circular symmetrical 
disturbance have been obtained. They show many resemblances to the corresponding distri- 
butions in tropical storms. Since the main source of energy of the mature tropical storms is 
the latent heat of condensation, it is really tempting to attribute the formation of these large- 
scale disturbances directly to convection in conditionally unstable air. However, the analysis 
shows that some organizing mechanisms are needed in order to channel the release of the 
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energy of unstable stratification into a large scale circulation. 


I. Introduction 


Free convections occur very often in nature, 
but are usually governed by different physical 
processes on different occasions. In earth’s free 
atmosphere such free convections are maintain- 
ed not by the mean temperature stratification, 
. which is usually stable, but by the water vapour 
distribution. If the temperature and the water 
vapour distribution in the atmosphere is such 
that when an element is lifted, enough latent 
heat may be liberated by condensation so as to 
raise the temperature of the element above that 
of the surroundings, then the stratification is 
unstable for ascent. This combined influence of 
temperature and water vapour stratification is 
most conveniently represented by the vertical 
distribution of the mean equivalent potential 
temperature Og. Instability requires saturation 
and an upward decrease of Og. The correspond- 
ing temperature lapse-rate is the moist adia- 
atic. 

Since most of the condensed liquid water or 
snow fall out as precipitation, a change of 
state in reverse direction occurs only in a much 
smaller extent in the descending branch of the 


1 The research was sponsored by the U.S. Weather 
Bureau under contract No. CWb-9939. 
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circulation. Therefore the mean stratification 
is usually stable for descending motion. Such 
a system which is unstable for saturated 
ascension and stable when unsaturated and for 
descension is said to be conditionally unstable. 
In such a system, we must treat the ascending 
motion and the descending motion separately. 

The problem of convection in a conditionally 
unstable atmosphere has been analyzed by 
BJERKNES (1938) and PETTERSSEN (1939) by the 
slice method, which takes into consideration 
the influence of the stable stratification in the 
descending region but neglects the pressure 
perturbation and friction. The results are that 
for a saturated atmosphere, whenever the 
lapse rate y is greater than the moist adiabatic 
lapse rate y,,, the stratification is unstable, but 
only those perturbations with a very small 
ascending region can grow. 

More recently, HAQUE (1952), Lizzy (1960) 
and Kuo (1960) applied the perturbation 
technique to this problem, and also take into 
consideration the Coriolis force because of 
their interest in the large scale motions. The 
first two authors treated the motion as inviscid 
and adiabatic, therefore were also unable to 
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determine a true critical value for 005/02 
or a preferred scale of the perturbation. In 
Haque’s discussion, only the largest pertur- 
bation with zero growth rate is considered, and 
he applied his result directly to the initiation 
of tropical storms. A similar interpretation has 
also been given by Lilly. However, such inter- 
pretations are highly questionable because it 
disregards all the other growing disturbances. 

The aim of this paper is to determine the 
critical 90£/9z for the initiation and the 
maintenance of convective motion in a 
conditionally unstable atmosphere, and to 
obtain the preferred scale and the nature of 
these convective motion. For these problems 
it is necessary to take into consideration the 
viscous effect and that of heat conduction, or 
some equivalent effects. It will be shown that 
in an atmosphere which is at rest and is either 
absolutely or conditionally unstable, only 
cumulus cloud type convection can be expected 
to evolve if random infinitesimal perturbations 
are introduced. In case of stable stratification 
in the descending region, the distance between 
the ascending centers will be relatively large, 
and may be as far as 20—40 kilometers apart. 
Such convective systems also tend to arrange 
into long bands when a shearing mean wind 
is present, as have been discussed by KUETTNER 
(1959), SoBHAG (1931), FRITZ and WEXLER 
(1960) and by WINSTON (1960). 

The second aim of this paper is to explore 
the possible mechanisms which will produce 
an organizing effect on the convective motion 
for the initiation and maintenance of a large- 
scale circulation, and help to control the small 
scale convection. A number of such factors, 
which seem to be important for the initiation 
and development of the large scale tropical 
cyclones, will be discussed in the last section. 

In addition, section two of this paper will 
be devoted to the derivation of a set of dynamic 
and thermodynamic equations which, even 
though simplified already, still govern almost 
all of the important atmospheric motions. 
Further simplifications can be obtained by 
introducing additional restrictions regarding 
the nature of motion. 


2. General equations 


In this section the equations appropriate for 
the class of motions in which the distributions 
of absolute pressure p and of density @ are 
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approximately as that of a normal atmosphere 
at rest shall be derived. 


We define the normal atmosphere as an 
atmosphere which has the same overall 
stratification, the same total mass and the same 
total internal energy as the actual atmosphere, 
and which is at hydrostatic equilibrium. 
Denoting by po, Qo, 19 the pressure, density, 
and absolute temperature of this normal 
atmosphere, then these quantities are functions 
of the vertical coordinate z alone, and p, and 
0, are related by the hydrostatic equation 


(2.1) 


where g is the gravitational acceleration. For 
simplicity, we assume that the atmosphere 
satisfies the equation of state of a perfect gas. 


Since Ty is the equilibrium temperature 
when the fluid is at rest, its vertical distribution 
is determined by nonadiabatic processes such 
as heat conduction, radiation and addition of 
latent heat of water vapor by condensation and 
evaporation. Heat conduction alone will give 
a linear dependence of Ty on z, while the 
other nonadiabatic processes may result in 
more complicated vertical variations. Here we 
merely consider Ty as a known function of z. 


In dealing with a gas medium it is more 
convenient to use the quantity s= log © as 
a dependent variable. Here © = T (P/p)!~v/ 
is the potential temperature, therefore s is 
proportional to the entropy. 

Denoting the departures of the various 
quantities from their normal values by small 
letters with a prime, and assuming that 


(2.2) 


we then obtain the following relations from 
the perfect gas equation 


P'IPo < 1, 0/00<I 


0"/00 = P'lPo — T/T, 
piles! ake 


ae (2.3) 
yP 
where y = c/c,, c, and c, being the specific 
heat at constant pressure and at constant volume 
respectively. Thus p’/po, @’/09, T’/T, and 
s’(= ©'/9,) are of the same order of magnitude. 
We mention that these primed quantities 
may have non-zero horizontal and vertical 
averages. 
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Hereafter the following symbols shall be 
used: 


— 

V = velocity vector with components 
u,v,w. 

>> 

i, j, k = unit vectors in the directions x, y 
2, with z increasing upward. 

a Se wee 4 
= = i —+ j — denotes the horizontal 

Lay 
delta operator. 
> 9 | 
v = %, k= denotes the three-dimen- 
dz 

sional delta operator. 

ie is the Coriolis parameter 

= = 

© = v x Visthe relative vorticity vector 
with components &, 7, ¢, given by 
£ Ow dv ou dw 

dy az" az ax’ 
- ov du 
= dx dy 

T = P'Î00 

v2 = 92]ox?+92/9y? and v?= v2 + 92/922 
denote the two-dimensional and the 
three-dimensional Laplacian. oper- 
ators. 

v = kinematic viscosity coefficient or its 
equivalence. 

k = thermometric counductivity or its 
equivalence. We shall approximate 
part of the radiation process by an 
equivalent conductive process. 

1 90 OT, 4 

ig) ee ÎT, is the sta- 
Oy dz Ainge 
bility factor of the normal dry at- 
mosphere. 

6 ” I 200 _ 4% „To ime = & 

Oy dz IN RIZ EHRT, 
+s. is the density stratification 
factor. 


Neglecting higher powers of p’ and 0’, the 
equations of motion can be written in the 
following vectorial form 


> > 
+k (s.u+gs')+v v2V 
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For carth’s atmosphere, the value of s, is 
of the order of 10-5 m=! whereas that of o and 
g/yRT, are about ten times larger. Since for 
most problems the scale height of the motion is 
usually smaller than H = 1/0 (= 10 km), we 
may neglect s, against either o, g/yRT, or 0/dz. 
Using this approximation, eq. (2.4) then 
reduces to the following 


=S 
EB 
T+ (fle+ 0) x V= -v(n+— Vi) + 


> = 
Dos ck a vere (2.4 a) 

The assumption that o’/o, is always small 
implies that 0/0, is approximately constant for 
a moving particle of fluid, i.e., 


I doe 5 20 ACC TOUL 
On dh NO, N dz 


Consequently the approximate form of the 
continuity equation is given by 


2 
veV=ow (2.5) 

This approximation implies a limitation on 
the time scale of the variations. 

Through the use of this approximate con- 
tinuity equation, sound waves and high fre- 
quency gravitational waves whose period is 
less than 8 minutes are excluded from the 
present consideration. We mention that when 
the vertical scale of the motion is much smaller 
than the scale height H( = 1/0), the equation of 
continuity reduces to that of incompressible 
flow. 

The first law of thermodynamics can be 
expressed more simply in terms of s: 

Zu 
2. +V-vs+s,w-k pages OK (2.6) 
ot CLS 
where Q is the accession of heat to unit mass 
per unit time by non-adiabatic processes other 
than by conduction. 

It may be mentioned that the equations (2.4), 
(2.5) and (2.6) are very similar to the equations 
in the (x, y, p, f) coordinate system, except 
that now the vertical equation of motion need 
not be hydrostatic. 

The variable x can be eliminated from eq. 
(2.4) by an application of the curl operator v x, 
which yields the vorticity equation. Combining 
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it with the continuity equation (2.5) it then 
takes the following form 


AU Lure 
(= )-Z vty 28225 
dt 00 00 Qo 
Pas at 
a el (2.7) 
0 


Equation (2.6) together with the three 
equations (2.7) and the appropriate boundary 
conditions suffice to determine the dependent 
variables u, v, w and s’. 

To determine x, it is convenient to obtain 
another equation by applying the delta 


operator v. to eq. (2.4 a) and making use 


of eq. (2.5), which result in the following 
equation 
nn (So) + 
dz dz 
2 3 ou; (2 8) 


It may be mentioned that if the horizontal 
scale of the motion is much larger than the 
vertical scale, the hydrostatic relation will be 
valid. Equation (2.8) then splits into two 
equations, one is the hydrostatic relation, the 
other is equivalent to the so called balance 
equation, and relates v?r to the wind field. 

For the determinations of the critical lapse 
rate and the preferred scales of motion for the 
relatively complicated cases of convection 
when the influence of the stable descending 
surroundings is taken into consideration, we 
shall find the energy integrals are extremely 
useful. Therefore we shall derive the energy 
equations at this point. ae 

Multiplying eq. (2.4 a) by @9V scalarly! and 
integrating over the entire volume, making use 
of eq. (2.5) and assuming that the normal 
velocity vanishes on the boundary we obtain 


à 
ama Jus dn =D (2.9) 


where K = — j V2 dm is the total kinetic 


1 If it is desired to use eq. (2.4), we may multiply this 
equation by 0,V@/@; integration over the total mass 
then yields a modified kinetic energy integral in terms of 

1 (VO 


1 Saar ar dm, where © is the mean value of ®.. 
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energy, D = - ‘i oe dm is the total rate of 
k 


dissipation, 7; being the viscous stress, and is 
v; 
given by tx = Ze where 1%) is the viscos- 


ity coefficient in the direction of x,, including 
that due to smaller eddy motions. The sum- 
mation convention over repeated indices i and 
k has been used in 7; and in D. 

Similarly, multiplication of eq. (2.6) by 
£00s /s- and integration lead to the following 
thermal energy integral 


JE Be s'Q F 
ref im = © [RS dm D, 


= £ ike dm and D,= ef ses dm. 


Sz 


(2.10) 


These two equations show that for inviscid 
and adiabatic motions, the total energy (K + E) 
is conserved. 

For the special case of adiabatic motion in 
an isentropic atmosphere, eq. (2.9) alone suf- 
fices to express the energy conservation rela- 


: a 
tion because we then have [ ws dm = = [ zsdm, 


which indicates again that the increase of ki- 
netic energy must be accompanied by an up- 
ward transport of s. A necessary but obvious 
modification of eq. (2.10) for this case is to 
remove the factor 1/s, from this equation. 


3. Differential equations for small perturbations 


Since the main concern of this paper is to 
determine the critical lapse rate and to investi- 
gate the nature of the convective motion near 
the critical point under the influences of earth’s 
rotation and the release of latent heat of con- 
densation, the governing equations (2.4 a) and 
(2.6) can be linearized. It is then possible, and 
also advantageous, to obtain a partial differen- 
tial equation for a single dependent variable. 

Since some kind of current is usually present 
in the atmosphere, we shall assume the 
existence of a constant mean velocity u,. The 
influences of a variable uy and of horizontal 
baroclinicity will be investigated later. 

Elimination of s’ between the third equation 
of (2.4 a) and eq. (2.6) yields 

(gs2 + LL’) w= a pasa £Q (3.1) 
924) body 
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d 0 à d 

= —— 2 4 = — — 

where L rs PY and L > t+ Mo 5 
-kv?. 


Note that U, occurs only in the operators 
s Pr: ) 

L and L’ in the combination U, = It can be 
x 


removed by taking a coordinate system which 
moves with the constant speed U,. Therefore 
for the disturbances ich are moving with 
Up, it has no effect.1 This term disappears also 
when the disturbances are independent of x. 

Applying the horizontal delta operator v,. 
to eq. (2.4 a) and making use of eq. (2.5) we 
obtain 


D) 
L(+ eu = —fE+vin (3.2) 


Lo 


Elimination of £ between (3.2) and the third 
of (2.7) yields 


I Joow 
(f? +18) Salve (3.3) 
Finally, eliminating x between (3.1) and 
(3.3) and introducing the new dependent 
variable y = 0 )we?*/? we obtain the following 
equation in y: 


: 9? 2 ; 
L'(f? +L?) (3-2 ) +L (gs. + LL)v?y = 
_ 89001? i 
QT, LviQ (3-4) 


For the earth’s atmosphere the value of 
02/4 is about 2.5 x 10-9 m? which is much 
smaller than 92/92?, therefore the term with 
o? in (3.4) may be neglected. 

In this study we shall be concerned only 
with the addition of latent heat due to conden- 
sation, which is considered to be determined by 
vertical motion alone. Thus we shall assume 

dm, om; 


Sel de Br: 


motion, 


for saturated upward 


(3.5) 


Q=o for unsaturated or descending motion. 


Here m, is the saturation mixing ratio of the air 
and L, is the latent heat of condensation. With 
this Q we may combine the term on the right 


1 The equation for nonviscous flow in a shearing mean 
current will be given in section 7. 
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side of eq. (3.4) with the last term on the left by 

introducing the stability factor S, defined by 
Meme teat AL al 5 Ar 7 

S=-;, Dear à (y — Ym)= -Se for 


(3.6) 


ifs 
S= -s,= 9 (y — ya) for unsaturated or de- 
0 


saturated ascending motion, 


scending motion, 


where Sz= 0,190;/9z, Or is the equivalent 
potential temperature and y,, and yg are the 
moist and dry adiabatic lapse rate. It may be 
mentioned that according to this assumption, 
S is always discontinuous across the plane 
separating the ascending and the descending 
motions, if the descending air is unsaturated. 

With this simplification, equation (3.4) 
reduces to the following 

2 2 

L(f? +12) (3 3 Ty) NS 


(3.7) 


The vertical distribution of s, and Sg in the 
West Indies region is given in figure 1, which 
are deduced from the mean soundings given 
by JorDAN (1958, 1959). 

We mention that if the hydrostatic approxi- 
mation has been used, the second L?L’ term will 
disappear from this equation whereas if a 
geostrophic balance is assumed for the x-equa- 
tion of motion, then the first L?L’ term will 
disappear from this equation. It has been 
shown in another paper (Kuo, 1960) that in 
discussing the stability problem for motions 
produced by unstable stratification we must not 
use these two approximations simultaneously; 
otherwise the possibility of instability will be 
completely lost. 

In order to avoid too much mathematical 
complexity, we shall only consider two- 
dimensional systems, either in the form of long 
rolls or with circular symmetry. In either of 
these two cases the variables depend only on 
one horizontal coordinate. The long rolls 
may be used to represent the long cloud 
strips, while the circular system may be used 
to represent an isolated system such as a single 
cloud, a tropical storm or a tornado. It may 
also be considered as an approximation to the 
hexagonal or even square patterns. This latter 
interpretation is particularly desirable because 
of the extreme complexity of the three- 
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Fig. 1. Vertical distributions of the stability factors S 

and the relative humidity h in West Indies during the 

hurricane season (July— October). The numbers in the 

parentheses along the h-curve indicate the lift (in mb) 
needed to produce saturation. 


dimensional solutions in case when the de- 
scending region is stable. 

Equation (3.7) must be solved with respect 
to certain specific boundary conditions. The 
external boundaries of the system under con- 
sideration are horizontal planes at z = o and 
Z= co. Sincep = 0 we 2 and 05, = 000€ 
where 099 is the mean sea level density, we 
shall require the vanishing of y both at the 
lower and at the upper boundaries. The latter 
is assumed in order to have the kinetic energy 
per unit volume approach zero! as z approaches 
infinity. Therefore two of the boundary 
conditions that must be satisfied by y are 


w=0 at z=0 and as z +00 


(3.8) 


1 For problems concerning motions in the upper atmos- 
phere, it will be more appropriate to require y to remain 
finite as z + ©. 
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For simplicity, weshall use y = 0 at z = o and 
z = h instead. 

When the viscous effect is taken into 
consideration in full, additional conditions 
must be imposed because the differential 
equation is of higher order. For simplicity, we 
shall consider the lower boundary as very 
smooth and therefore require the vanishing of 
tangential stresses and of w and s (the so called 
free boundary conditions). These conditions 
can be satisfied by harmonic solutions. 


In the horizontal directions the fluid is 
assumed to extend to infinity. The conditions 
we impose in the horizontal direction is 
either finiteness or periodicity of the solution. 
Specifically, for the cases of circular symmetry 
or long rolls we shall require the vanishing of 
the radial velocity u at the axis of symmetry. 
We shall also require u to vanish at r= fo, 
which may be finite or infinite. In addition, 
dw/dr is also assumed to vanish at these two 
radii. Thus the external boundary conditions 
in the r-direction are: 


y=oand vy=oatr=oandr=r, 


(3-9) 


For inviscid flow these two conditions are 
identical and therefore there is only one con- 
dition for each given value of r. 


Since in the problems treated below the 
stability factor varies discontinuously (either) 
in the horizontal (or in the vertical) directions, 
certain conditions must be imposed on the 
internal boundaries. The necessary require- 
ments are the continuity of the normal 
velocity v, and of the pressure across these 
boundaries. In terms of y, the first requirement 
is equivalent to the continuity of y, while the 
second depends upon the S variation, and may 
be derived from the third equation of (2.4 a) 
and eq. (2.6). 

Thus for inviscid flow the continuity of p 
across a vertical surface requires the continuity 
of (gs - q®)w, therefore the condition is given 


re 


by (from eq. [3.1]). 


where q is the time rate of amplification. There- 
fore w changes its sign abruptly across the 
surface, as has been found by Haque (1952). 
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On the other hand, for steady state viscous 
flow this requirement is given by 


vviw +91 = VV Wa +95; for x=x, (3.11) 


In the stable descending region the motion is 
usually very weak, therefore the hydrostatic 
approximation should hold in this region. The 
first term on the right then drops out and we 


have 

vviw, +85 = 853 (3.278) 
Since the velocity must be continuous in 
viscous flow, we must also require w to vanish 
on this surface. Hence we must require 


Uy =U, w=w=oatx=x, (3.11 b) 
It should be pointed out that the internal 
boundary x = x, which separates the ascending 
from the descending current is a free boundary. 
The position of this boundary is not given, 
but must be determined by the solution. This 
makes the problem somewhat more difficult 
as compared with the other problems which 
involve only fixed internal boundaries. 

When S changes discontinuously in the 
vertical direction, the inviscid equations require 
the continuity of both the horizontal and the 
vertical velocity components, therefore con- 
tinuity of both y and dy/dz. 


4. The growth rate and most preferred scale 
in an unstable layer of constant S. 


The nature of convective motion in a single 
layer of constant S is relatively wellknown from 
the previous studies of HAQUE (1952), CHAND- 
RASEKHAR (1953), LILLY (1960), and Kuo (1960). 
Here we wish to recapitulate some already 
known results to facilitate the discussion of the 
convective motions for a variable S. For this 
case, the solution of eq. (3.7) that satisfies the 
free boundary conditions is of the form 


(4.1) 


where F (x, y) is a harmonic function satis- 
fying the equation 


cake 
y = Ae sin a F (x, y) 


Ta 
DEREN: 


2 (4.1 a) 


viF= 


The length scale being used is a depth h, and 
& is the horizontal wave number and k is the 
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vertical wave number. In dealing with at- 
mospheric motions we should use the eddy 
coefficient of viscosity, which in general has 
different values in different directions. We 
shall denote the eddy viscosity coefficient in 
the horizontal directions by », and that in the 
vertical direction by y,. 

Different solutions of (4.1 a) can be obtained 
to represent different plane forms of the 
convection cells. For example, the solutions 
that represent long cloud rolls and circular 
cells and rectangular cells are given by, 


respectively 
ATX arır 
Ba. Bel) 
Fe = sın ae sin i (4.1 b) 


When the solution (4.1) is substituted in eq. 
(3.7) we find 


20 _ L2fe\ % 2 
len a 2 (x? », +k? v,) 


k? + a2 h? (42) 


Thus q has its maximum for k = 1, corre- 
sponding to the first mode in the vertical. We 
shall use this value of k only. Then q is given 


by 
20 _f2\% 2 
= (© 2 £) U (vz + ayy) 


Ita? 


Eq. (4.2 a) shows that for the smaller distur- 
bances with a horizontal wave number « 


larger than f/V. gS, q is real and therefore the 
disturbance will grow if the first term of eq. 
(4.2 a) is larger than the frictional effect, and 
will be damped out gradually if the second 
term is numerically larger than the first term, 
which is the case for the very small disturb- 
ances. On the other hand, q is complex but has 
a negative real part for the larger disturbances 
whose horizontal wave number «is smaller than 


f/V.gS. Hence these larger disturbances behave 
like the stable gravitational-internal long waves 
in a stable atmosphere, and may propagate 
with a certain speed. To demonstrate this let 
us put 


(4.3) 


BRAZOS 
wy = Ae”tt sin sin > (ax + &y — nt) 


Substituting in (3.7) we find 
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Fig. 2. The growth rate q as a function of the horizontal 

wave number «& for three different values of S. Dashed 

curves given by inviscid approximation, solid curves 
given by viscous equation for v = 10% m? sec7}. 


fr — 02S 778 
n? = SEE, b= m (k®v, + 0®v,) (4.4 a) 
where «? = «+2. Therefore the phase veloc- 


ity is given by 


Cre n _(f?-g025)% 
Tr Roe 


(4-4) 


Since the phase velocity depends on both 
frequency and direction, the propagation is 
both dispersive and anisotropic. 

From equations (4.2 a) and (4.4 b), we see 
that for any given positive S there exists a 
limiting horizontal wave number «, given by 
as = f?/gS which separates the amplifying 
shorter waves from the decaying long waves. 

The physical interpretation of this limiting 
wave number is that for the perturbations with 
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smaller than 9,~the kinetic energy of the 
horizontal motion produced by the Coriolis 
effect through a displacement is greater than 
the potential energy released during this 
displacement, while the reverse is true for 
perturbations with « greater than «4. ; 

Many meteorologists have interpreted this 
limiting wave number as representing the 
preferred scale of motion for the given S, 
and applied it directly to the hurricane for- 
mation (see HAQUE 1952, SYONO 1953, and 
Lıtry 1960). However, since the perturbations 
with «>a» are growing with time while the 
limiting perturbation is not, there seems to be 
no reason to expect this limiting perturbation 
to become the dominant one, at least not 
according to the linearized equations. So far as 
linearized equations are concerned, the most 
preferred scale of motion must be the one 
whose growth rate q is a maximum, or the 
one with a lowest critical stability factor 
S = Sy for marginal stability (q = 0), if we 
consider S has been raised gradually from 
below to above S,,. In this section we shall 
examine the preferred scale from these two 
different interpretations. 

The values of q as given by eq. (4.2 a) for 


f= 5 x 1075 sec”! (corresponding to ¢ = 20 


deg. lat.), h = 10 km and for three different 
values of S which often exist in the tropics 
are plotted in fig. 2. The solid curves correspond 
to v, = » = 10% m? sec-!, while the dashed 
curves correspond to inviscid motions. It is 
seen that when friction is neglected, q has 
its maximum when « = co (L = 0), whereas 
with friction the maximum occurs at about 
& = 1.6, corresponding to a horizontal scale 
r, = 0.76 h for long rolls and r, = 0.48 h for 
a circular cell. This preferred wave number 
shifts very slowly toward a higher value as S 
increases. Thus the higher is S, the smaller is 
the preferred scale of motion. 

Comparing the dashed curves with the 
corresponding full curves, we see that the 
effect of eddy viscosity on q is small for & <1. 
We also mention that except for very small 
values of x, the effect of the Coriolis parameter 
fis also small. 

To obtain the lowest value of S for the 
maintenance of convection we set q =0, 
k = 1,9, = % =» in eq. (4.2). 

Introducing the nondimensional quantities R 


and T, defined by 
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we then find the following relation 
Roe@=(1 +0284 7 (4.6) 


which is identical to the equation obtained by 
CHANDRASEKHAR (1953). From this equation 
we find that for large T' the critical value of R 
and the most preferred « are given by 


"a 
Rn 23 (2) 
2 
(2) 
ce = 
2 


Numerical values of R,, and «, for this case 
have been computed by CHANDRASEKHAR 
(1953). We note that for fixed », the critical 
temperature gradient S,, increases with f" 
while for fixed f, S,, is proportional to »"» 
if we keep the Prandtl number constant. Thus 
in this case, the horizontal scale of the motion 
will shrink indefinitely as the value of T 
increases. 


(4.7) 


5. Convection in conditionally unstable atmos- 


phere 


In order to avoid matkematical complications 
in this section we shall assume that the con- 
vection cells are of the form of long rolls, with 
their axis in the direction of y, so that all the 
variables are independent of y. Circular cells 
shall be treated in the next section. 

The stability factor S is supposed to be 
independent of height while its horizontal 
variation is given by 


Ses, co for O <7 = x, 
StS, = -l28, for x ee; 


where x = x, is the boundary between the 
ascending and the descending branches and 
2x, is the width of the cell. The function y 
: : . TZ 
is assumed to be given by y = ®(x) sin are 
therefore it is only necessary to solve for the 


horizontal function D(x). 


A. Nonviscous solution 


The nonviscous solution for the case with 
a stable descending region and q = o has been 
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discussed by Haque (1952) and by Lırıy 
(1960). The two most important results 
obtained by these authors are (i), the vertical 
velocity has a finite discontinuity across the 
surface which separates the ascending from 
the descending region, and (ii), the ratio of the 
area of the ascending motion to that of the 
descending motion descreases as the ratio 
l?— — §,/S, increases. These conclusions also 
hold approximately when the viscous effect 
is taken into consideration. Because of the 
importance of these conclusions, we shall 
redevelop the nonviscous solution briefly and 
obtain some additional information. 


When » = 0, the solution ®(x) is then given 
by 
D, = À sin «& foro<é<é, 
“>, sing y sinh Pl&s-&) 
sinh BE, — &) 


(5-1) 
for &, <&<&, 


P, 


where & = a and B are defined by 


ra % 
8S1 — @ Prt in 


C= 


Substitution of (5.1) in the boundary condition 
(3.9) yields the compatibility equation 


tanh BAE&=I' tan af, (AE=& — &) (5.2) 


which is similar to the one obtained by Litty 
(1960), except that here I’, « and B all involve q 
and that &, appears explicitly. A similar 
equation holds for circular disturbances. 

Equation (5.2) has an infinite number of 
roots of which only the first is of interest 
to us. This root can easily be computed for 
given values of BAË. The results show that 
as BAË increases from 0.5 to 1.5, the value of 
a&, (which measures the area of the ascending 
motion) is doubled, whereas for BAË > 1.5, 
a&, changes only slightly. 

When BAË is large, i.e., when BAE > 1.5, 
eq. (5.2) reduces to 


x 
= tan? af, 


li 
ji (5.2 a) 


where x = q?/gS,. On the other hand, when 
Baé is large but «é, is small, we obtain 


ver (+1) 


(5.2 b) 
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which shows that the amplification rate 
decreases with increasing / and increasing axé. 
When both BAË and «£, are small, eq. (5.2) 


reduces to 
x= 1 /?/(AË/E) (5.2 c) 


which shows that the area ratio AË/E, must 
be greater than I? in order to have amplification 
“and that the larger is the ratio A&/&, the 
larger is x. 

Since eq. (5.2) contains « and &,, we need 
another relation between q and the other 
parameters in order to determine q and &; as 
functions of «. Such an equation can be obtain- 
ed by substituting the quantities u, v, w and s 
corresponding to the solution (5.1) in the 
kinetic energy integral (2.9), taken over the 
volume of a single cell. For the case of large 
B46 this relation is given by 


sin 2a& 
x (1 +02) - 2 + a2 = — 
20.1 


le) = pale x) (5.3) 
where a = f?/gS,. According to (5.1 a) the 
right side member of this equation is zero. 
Therefore we have the remarkable result that 
q depends neither on / nor on &, and that x 
is a function of the horizontal wave number 
a of the ascending region alone, given by 


a = ae 


et + œ2 (5-3 a) 


This equation is identical with eq. (4.2 a), 
except that the frictional effect has been 
neglected. Even though this equation is 
obtained under the assumption that BAË is 
large such that tanh BAË is nearly equal to 1, 
it is believed that it also holds approximately 
under more general conditions. 
Substituting (5.3 a) in (5.2 a), we obtain 


chef sail! BiB coke a ena 
tan és = | Dana (64 


This equation determines &, uniquely as a 
function of « for given values of I. The values 
of &, as given by this equation for various 
values of « and for / = 1 and! = 2 are given 
in the following table. It is seen that for the 
same &, €, decreases with increasing /, indicating 
that the influence of the descending motion 
in the stable region is to reduce the area of the 
ascending motion. 


EIS ERS.) 


af as 
Table 1. Values of &, | =) as given by (5.4) for 


different values of & and of I. 


RG 0.01| 0.02 0.05| OnLy 0215.05] | ie O| aan 
ı | 78.5] 39-3] 15:7| 7-8] 3-83] 1-37] 0-52] 0.24] 0.16 
2 | 46.4] 23.2] 9.3] 4.6] 2.32| 0.83] 0.32] 0.15] 0.10 


B. Viscous solution and stability criterion 


When the viscous effect is taken into con- 
sideration by the Navier-Stokes equations, it is 
difficult to obtain the growth-rate explicitly. 
On the other hand, the concept of marginal 
stability becomes meaningful, we shall there- 
fore derive the stability criterion according to 
this concept. For this case, the function ®,(£) 
for the ascending region is the solution of the 
following equation 


(D?- 1)? ®,- R D? ,- TG,=0 (5.4) 


where R is the Rayleighnumber for the unstable 
ascending region, and is defined by (4.5), and 
D? = d?/d& for long rolls and .D? = = > 4 
Ss 
for circular cells. 

In the descending region it is reasonable to 
assume the validity of the hydrostatic approxi- 
mation. The equation for.®, is then given by 


(D? — 1)? ®6,+ T@,= R,D? ®, (5.5) 


where R, = gs» h*/kva*. For convenience 
sake, we set Ry, =/?R so that only one 
stability parameter R occurs in the solution. 

Exact solutions of these two equations can 
be obtained, and when the external and internal 
boundary conditions are applied and the 
arbitrary constants have been eliminated, we 
shall find a transcendental equation relating 
R to the scale of the motion. However, such 
an approach is too complicated for the present 
problem and therefore will not be attempted 
here. Instead, we shall only find approximate 
solutions of the equations (5.4) and (5.5) and 
make use of the kinetic energy integral to 
determine the critical value of R and the pre- 
ferred scale of motion. Thus, this method is 
similar to the variation method. 

In this section we deal only with cloud rolls. 
The function ®, is assumed to be given by 


(5.6) 
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If D, is to satisfy eq. (5.5) exactly, then «, must 
be a root of eq. (4.6). However, since only an 
approximate solution is to be used, we do not 
take eq. (4.6) exactly but merely as a guide for 
choosing the three &, in ®,. Therefore we 
write the equation for &, as 

(r+o2)8+T-Ra2=0 (5.4 a) 
Let us arrange the three roots of this equation 
according to their order of magnitude. For 
relatively large R and not too large T we then 
have 


Hence the magnitudes of — «2 and of a? are 
much larger than that of «2. For example, for 
atmospheric conditions we find R is of the 
order of 10? to 10°, while T is from ı to 100. 
Under these conditions the approximations in 
(s.7) are certainly satisfied. We shall assume 
this to be true even for near critical values of R. 

Since &, is real and is much larger than «, 
when R is large, A,/A, must be very small in 
order to have ®, remain positive for o < & < &.. 
For simplicity we take A, = 0. On the other 
hand, x is imaginary and sin «, & increases 
exponentially with & Therefore the solution 
in the ascending region is given by 
ÉCOLE. 


@,=A sin a& — — = sinh us| (5.8) 


u cosh wé 


where & is used to represent «,. This solution 
satisfies ®, (0) = o and w,(&,) = 0. 

It may be mentioned that the hydrostatic 
approximation cannot be used in this region 
when friction is included because it gives a 
rapidly oscillating solution. Similarly, we 
find the solution of eq. (5.5) that satisfies the 
conditions (3.9) at £ = o and & = &, and the 
condition (3.11 b) at & = £, is given by 


_ By cosh B,aë 3 
Ba cosh Bag 


ct Ure 


Dy = BI hf (EE) 


sinh B, (é - a) 


ey Cs ad, By mw tan «&,— x tanh B,AË 
“cosh B,AË u B, tanh B,4é - By tanh B,aé 
(5.9) 
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where A&=£&,-£,. Here B, and B, are the 
roots of the following equation 

(Bis 1) -PRA,+T=o (5.10) 
Since R, is defined by the normal stability 
factor sz, of the dry air, it is usually of the 
order of 105. Therefore the two roots are 
approximately given by 


83 T+ı a 


c= Res pn 


3=-PR+2>1ı (5.102) 


On the other hand, if R, is very small, we 
shall have 

(ae ep. tT a gto) 
This case will be excluded from the present 
consideration. 

Note that when R is large, both «/u and 
B1/B2 are much smaller than 1. Then the second 
terms of the solutions (5.8) and (5.9) contribute 
very little except very near to the surface of 
discontinuity. Since the first terms of these 
solutions are the same as the nonviscous 
solution (5.1) except that in (5.1) « and B 
involve the amplification rate q, we may call 
these first terms the nonviscous approxi- 
mations. 

The variations of w and s at the level 
z = h/, in the two regions are illustrated in 
fig. 3 by the dashed lines for the cases of | = 1 
and { = 2, for a large R and an infinite aé. 
The fully drawn curves represent the non- 
viscous solution. This example represents a 
circular symmetrical disturbance, which we 
shall discuss in section 6. It is seen that for / = 1, 
strong descending motion takes place close to 
the internal boundary and decreases rapidly 
with increasing &, whereas for larger ! (more 
stable surroundings), the descending motion 
is much weaker but more wide spread. It is 
also seen that as / increases, the ascending area 
diminishes. The contribution by the second 
terms of (5.8) and (5.9) are very small except 
very close to the internal boundary & = §). 
However, as will be shown in the next section, 
the viscous terms are important when R is 
relatively small and that they contribute signi- 
ficantly to the dissipation. 

For comparison, the variations of w, s’ and 
pas given by the nonviscous solution for the 
case of unstable stratification everywhere has 
also been represented by the dashed curve w, 


Fig. 3. Variation of w with & for / = 1 and / = 2 in con- 

ditionally unstable atmosphere. Solid curves given by 

inviscid solution. Dashed curves given by viscous 
solution. 


in fig. 3. It is seen that for this case s’is negative 
in the descending region whereas in the con- 
ditionally unstable atmosphere s’ is positive 
everywhere. 

In order to apply the boundary conditions 
(3.11 a) we must obtain s, and s, by solving 
eq. (2.6). When the variation of 0, is disregard- 
ed and ß, and B, differ from 1, these solutions 
are given by 


ANSATZ 
= Rn 7 
cos «Ë cosh u& cos af, 
i Fa: eT cosh ra 


= 


Dhs 0) 2 ee BalË — &,) I 
Koi taal 15 HA © Blox 
cosh B,AË 


cosh sé cosh B, (§ - :)| (5.11) 


Substituting these functions in (3.11 a) and 
neglecting the rather small terms involving 
a/u and B,/B, we obtain 


CLIC 
BA re 

Dass» utan ad a tanh pea 

u By tanh B,AË - B, tanh B,AË 
We are particularly interested in the cases 
when R is relatively large so that the approxi- 


mations (5.7) and (5.10 a) are valid. Eq. (5.12) 
can. then be simplified to 


ee 
Mat Ruh 


(5.12) 


Dep 
———tanh B,aé= 
I + a? Pde 


I (tan &&, - nb ué:) Gina) 


HPEARUO 


This equation differs from eq. (5.2) by the 
factor (1 — B?)/1 + «2. Note that this equation 
involves the two unknowns « and &, besides 
the ratio | which is supposed to be known. 
Therefore another relation is needed for their 
determination. However, since both «? and 6? 
are usually much smaller than 1, a first approxi- 
mation can be obtained by omitting them 
against I, we then have 


(5.12 b) 


The first root of this approximate equation 
are given in table 2a for various values of 
B1AË and I. It is seen that for a given |, a, 
decreases with decreasing B,AË, indicating 
that the placing of a fixed wall at a finite 
distance reduces the area of the ascending 
motion. However, this effect is prominent 
only when B,AË is less than 1.5. For B,AË 
greater than 1.5, «&, is almost equal to its 
optimum value corresponding to B,AË =. 
It is seen that a&, decreases rapidly with 
increasing |. 


tanh B1AË= Î tan ad, 


Table 2a. Values of «£, as given by (5.12 b) for different 
values of B,AË and of I. 


1 
so 


2 3 5 Io 

oo 0.786] 0.464| 0.322] 0.197| 0.100 
2.5 .780 .460 .318 .196 .099 
2.0 .765| .450| .310| .190| .096 
1.5 .735 .425 .295 .I180| .090 
1.0 .650 .360 .250 .150 .076 
0.5 .431 .230| .155| .093| .047 


In order to make use of the kinetic energy 
integral (2.9) we must obtain v, and v, by 
solving the y-component equation of (2.4 a). 
The results are 


_ An T% TZ 


V1 07 
sina& asinh wé cos a&, 
= + es u mall 
( + (Pag) cosas, (5-13) 
BuT2 nz 
Va = the cos 7% 


(ae Ba (E — 2) Fs Bisinh B,(£ - &,) cos ati) 
La B2(B2-1)  coshB,aë 
Substituting the velocity components u, v, 


w and the quantity s given by these solutions 
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in eq. (2.9) and setting 2K/0f to zero we obtain 
a rather complicated equation relating R to 
the various parameters. Only a much simplified 
version of this equation shall be given here, 
which is based on the assumptions that u>«, 
Ba> By, and 6,46 is larger than 1.5, so that 
we may put tanh B,AË = tanh ßza& = 1. Under 
these circumstances the ratio B/A of (5.9) 
reduces to 


Bs ae” I+a* cos aé, 
A sinh B,4E I(t- 2) cosh B,AË 
by (5.12 b) 


and the equation resulting from the kinetic 
energy integral takes the simple form 


sin 2&Ë; 
24Ë: 
where F= (3 + 5a? +1) +(1 + 0?)a® w+ 

+2 Bıße- Bi- Bi 


According to the equations (5.7) and (5.10 a), 
we have 


we R=(1+0)8+T+F 


(5-14) 


T+1)* 
au Mr 1) ; 


R% 


6,8. = (T+ 1) 


Therefore eq. (5.14) is the other equation relat- 
ing R to the quantities « and «&,, besides the 
known quantities T and /. Together with 
eq. (5.12), we may solve for R in terms of «, 
and the preferred horizontal wave number 
%» and the minimum R,, can be obtained for 
every given set of T and /. 

Because of the complexity of the computa- 
tion, only the first approximations of «,, and 
R,, have been obtained. The procedure 
adopted is to compute the first approximation 
of «&, from eq. (5.12 b) for B,AË = 00. These 
are given in table 2 b below, together with 
the corresponding function I= sin 2 «&1/2«&,. 


Table 2b. Values of «&, for B,AË>1 : 5 and various | 
as given by eq. (5.12 b), and the corresponding I. 


I | I 2 | 3 5 Io 
adı 0.785| 0.462] 0.320] 0.198] 0.100 
vf 0.637} 0.863] 0.933} 0.975] 1.000 


The approximate values of «, and R,, have 
been obtained and are given in table 3 for 
various values of | and T. The «,, values change 
Tellus XIII (1961), 4 


only slightly with I, therefore only the mean 
is given in the last line. These results show that 
Ry increases both with T and with I. For 1> 1, 
the R,, values are about twice as that required 
by a layer of fluid of constant R. 


Table 3. Values of Ry as functions of T and 1. 
REEL 


12.5| -23.0| 20,5, 02.0 | O2-21) 225) | asia 
15.21 20.0) 35.517 74.7) 120-0) 270 | 412 
18.6] 35-6] 46.2} 97.0] 139.0} 349] 540 

O07 O:95|| 00) 7730 77.5316. 1.931.219 


50 100 500 | I 000 


Sun. 


Because of the presence of the factor 
(x — 62)/1 + a on the lett hand side of eq. 
(5.12 a), the true values of «&, are somewhat 
smaller than that given by eq. (5.12) in table 
2 a. On the other hand, due to the precence 
of the second term on the right side of eq. 
(5.12 a), the correction is actually small in 
most cases. 

When / is greater than 3, the function I is 
very nearly equal to 1. The R,, and «,, values 
for | = 5 in table 3 are computed from this 
approximation. 

For very large values of T and of I, we have 
the asymptotic formula 


1 
R=4(1ta®)’+ ; (5.15 a) 


From this we find for large T and large |: 
Russen 


I 
A? ea Th 
2 


Thus these asymptotic relations are similar to 
that given by eq. (4.7), but with different 
coefficients, therefore giving a larger R,, and 
a smaller &,. 

It may be mentioned that if the approximate 
relation (5.12 b) is used in the kinetic energy 
integral, it will result in the following equation 
(Kuo, 1960) 


24 2 oe + fa 
Ra 1° = FE (1+@2)°+ r(: ee an) 
+ (Ir a) I {304+ 202 = B4+ let + 28,89} 
(5.17) 


This equation differs from eq. (5.14) in two 
important aspects. First, the horizontal wave 


454 ala Abe 
number « must be greater than zero but less 
than 1, secondly the asymptotic behavior for 
large T and large lis Ry, 5.83 T, %n >0.643. 

We mention that the critical values of R 
given in table 3 are rather small when compared 
with the R values observed in the atmosphere 
when the unstable layer is more than one 
kilometer deep. Under those conditions per- 
turbations of different scales will be able to 
grow simultancously and the motion must be 
turbulent and the linear theory will not be 
quite valid. 


6. Circular disturbances in conditionally un- 
stable atmosphere 


The next mathematically simple plan form 
of convection is that of circular symmetry. 
This form may be used either as representing 
an isolated convection cell or as an approxima- 
tion to the system of hexagonal cells, therefore 
it is of interest to discuss this case in more 
detail, especially in regards to the space 
distribution of the various quantities. For cells 
of this form, the function ® for the two 
regions are given by 


Gee (ees) MORE, LES 
rer LS 
[eT (ime) " (as) 


,= BE |Z, (ie) - Ze (Bit) 


| Bs Zo (iBaés) 


Zy (i)| 
(6.1) 


where 7, and r, are Bessel functions of order 
zero and order one, respectively, while « and 
u and ß, and f, are given by eqs. (5.7) and 
(5.10 a). The functions Z, is given by 


A 1(%2) 


Ar icon 


FH (ee (6. ra) 
where H! and H? are the two Hankel functions 
of order one. Z, is obtained from Z, by 
differentiation. These solutions satisfy the con- 
ditions ®,(0) = ®,(£;) =o and w,(é,) = 
Wa(£:) = 0. The ratio B/A is determined by 
the condition D, = ®, at & = &,. When &, is 
very large, Z, reduces to H! alone and Z, = 
tH. 

Following the same procedure as employed 
in the preceding section, we obtain s and the 
tangential velocity v in the two regions by 
solving eq. (2.6) and the tangential component 
of eq. (2.4 a). The perturbation pressure is 


KUO 


then obtained from the vertical equation of 
motion. The results are 


dSo1 F, ie CR Lite a (x81) 
dz T+a2 u2-ır, (iuË) 


Sy == Aa 


To (us)} 


ge .n Poe nn. I 
m Dr la 


Zo (iB1&1) 2a 
Za (ine) Pe} (629 
_ ap (tle) __« tala) 


2 5] Zu lie), A a 
a F,| b Ba (B3 - 1) Zo (iB 2&1) 


Bees i 
ZUBE) + Gre (a) (6.2b) 
HOUR ge] ey Tr 
F= sn CE OS (6.2 c) 


The as yet unspecified constant C in vg is 
determined by the continuity ofvat& = &.. 
Note that this term does not contribute to 
the dissipation. 

Substituting DED pss and S2 in eq. (ra 
neglecting the small terms involving 81/8: and 
assuming that &, is large, we obtain 

2 ‘BE 
LIL Br Hi (és) 
271 (IRE 
1+ où (HG (1038) 
& (xé:) 7 IT: Ke) 


To (xé:) UT (iu$,) 


tanh B,A&=]1 


(6.3) 


which is the counterpart of eq. (5.12 a). 

The values of the first root of this equation, 
with the factor (1 - ß})/t +a? omitted, are 
given in table 4 for various values of B,AË 
and I. These results again show that «a&, 
decreases both with increasing B,AË and with 
increasing I. Comparing with the values in 
table 2a we notice that here x€, decreases more 
slowly as / increases. However, since for this 
case the area is proportional to &®, we find that 
the ratio of the area of the ascending region 
to that of the descending region is descreasing 
more rapidly with / than for long rolls. For 
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example, for | = 3 and B,AË = 1.5, this ratio 
is 15 for long rolls and is 26 for circular cells. 


Table 4. Values of œ£, as given by (6.2) for different 
values of B,AË and I. 


re 
| 


|“ 


2 3 
oo TOOLS|) 71.28 Te 0.9781..0.83 
2.5 1.0605] 1:27 ET 0.96 | 0.83 
2.0 1.59 1.26 I.IO | 0.95 0.82 
TR 257 | 2.25 | 71.09, 0.94 | 0,81 
1.0 1.45 1.70, 11. 0:907 | 0:84 | ovr 
0.5 7700 0,822 17.0.7251 .0:657110.54 


Another relation between R, x and E can be 
obtained by an application of the kinetic 
energy integral. The form. of this equation is 
similar to that for the long rolls and it yields 
nearly the same numerical results and therefore 
will not be given here. 

As has been mentioned at the end of the 
preceding section, when a deep unstable layer 
exists in the atmosphere, the actual R is usually 
many orders of magnitude larger than the 
eritical value. Under such conditions not only 
the disturbance with wave number «,, but a 
large band of disturbances of widely different 
scales may grow simultaneously; the motion 
is more likely to be turbulent. The upper limit 
of the horizontal scale of this band of amplify- 
ing perturbations is approximately given by 
Gate le Slot To >see Les = 
SotGn 4, „we. find. &45>=:0:007.... The 
radial dimension of the ascending region of 
this limiting disturbance is given in table 5 for 
various values of /. These values compare 
favorably with that of the rain areas of the 
tropical storms, which suggests that convective 
process is relevant to the dynamics of these 
storms. 


Table 5. Radius of ascending region of the limiting 


perturbation. 
I | I | 7 | 3 | 5 | Io 
rime?) | 362 288 | 252 218 | 187 


Now we shall discuss the distribution of the 
various quantities in space as given by the 
solutions. For this purpose we shall choose two 
different examples, one is super-critical (R> 
R,,) and the second is nearly critical. 
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Fig. 4. Radial variations of the mid-level vertical velocity 

(W), the surface tangential wind (V), the surface radial 

wind (-U), the midlevel entropy (S) and surface pressure 

(P) perturbations for a disturbance in a super-critical 
stratification. 


Example tr .R =! 1000, [ =-100, a = 0,31, 


u = 5.62 
Re = Ry, [= I, By = 0.31, 
Boe 306: 


Since R is super-critical, « and u are simply 
given by eq. (5.7) whereas B, and Pf, are 
given by eq. (5.10 a). For this case &/u and 
Pı/ß, are small and the second terms of the 
equations (6.1) and (6.2) are much smaller 
than the corresponding first terms except near 
the surface of discontinuity £ = &,. This is 
also true for all the large scale disturbances, 
particularly for the limiting disturbance « = 
Na. since a" and D ate small. ue 
root of eq. (5.12 a) is nearly the same as that 
of eq. (5.12 b), viz, «6; = 1.614 for ByAE= 
1-5: 

The radial variation of the mid-level (z = 
h/2) vertical velocity for this case is represented 
by the curve W in fig. 4, while the curve 
(—U, V) represents the ground level inward 
radial and tangential velocities. If we take 
W-o as unity, then the unit of the radial 
velocity u is «=! while that of the tangential 
velocity v is T*/a. 

These curves show that for this case w 
changes rapidly across the surface £ = &, 
while u and v are continuous, have their 
maxima occuring here and their radial deriv- 
atives are discontinuous. Since both u and v 


s TZ, 
are proportional to cos 7 the lower part 
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Fig. 5. Radial variations of W, V, -U, P and S for a 
disturbance in nearly critical stratification. 


of the cell has an inward flow and a counter 
clockwise motion while the upper part has 
an outward flow and clockwise tangential 
motion. The vorticity is anticyclonic at all 
levels for large r. 

The radial variation of s at the mid-level and 
that of the perturbation pressure at the top 
are represented by the curve marked with 
(P, S). Because of the smallness of &® and ß?, 
the pressure distribution is nearly hydrostatic 
and the two curves P, S coincide. Since p’ 
TZ 

: 


is proportional to cos = we have low pressure 


in the lower part and higher pressure in the 
upper part. On the other hand, the temperature 

eparture is positive at all levels and all radial 
distances, with its maximum. occurring at the 
mid-level and at the center. 

It may be mentioned that when « and B; are 
small, the radial distributions are nearly the 
same as that given by the nonviscous solutions 
except very near the surface of discontinuity. 


Example 2. R = 16, T=0, & = a, = 0.707; 

LH = 2. 

R, = 64, | = 2, By = 0.3535, 
Da = 8. 


For this case the root of eq. (site b) is «Ë, = 
1.28. On recomputing from eq. (5.12 a) we 
find the corrected value of «£, is 1.24, therefore 
the correction is quite small. 

The radial variations of the various quantities 
are represented in fig. 5. Comparing with the 
fully drawn curve in fig. 3 for 1 = 2, we see 


IeO Mes NIOKO) 


that the present solutions deviate very much 
from the corresponding non-viscous solutions. 
For this case w decreases more gradually 
across the surface £ = &,, and U and V are 
represented by two different curves. The 
maximum radial velocity occurs inside £ = &, 
while that of the tangential velocity occurs 
just on this surface. The pressure deviates 
appreciably from the hydrostatic pressure in 
the ascending region, therefore p and s are 
represented by two different curves. Note that 
s varies only slightly in the ascending region. 

It may be pointed out that even though the 
solutions give different space variations for 
different «, all of them have the same overall 
character, that is the temperature departure is 
positive everywhere and the low level pressure 
is low compared to the undisturbed surround- 
ings. Because these features bear so much 
resemblance to that of tropical storms, one is 
tempted to consider that the origin of the 
tropical storms is simply conditional instability, 
and to identify them with the “limiting 
perturbation” discussed above. However, since 
the most preferred wave number is %, which 
is much larger than «, some other physical 
processes must be operative in order to produce 
a distinct large scale system. In section 8 we 
shall discuss some such factors which favor a 
large scale circulation. 


7. Plane forms and horizontal spacing 


When the stratification is unstable every- 
where, both the critical value R,, and the 
preferred horizontal wave number «,„ are 
given by the vertical solution while the 
horizontal solution determines the plane form 
and the spacing of the cells. For the case of 
smooth bottom and top boundaries and very 
small T, &. = 0.707. 

For long cloud rolls, the width A of a cell, 
which is composed of two Opposite vortex 
rolls, as is drawn in fig. 6 a, is 


A= 2V2h (7.1) 


It can easily be shown that the distance 
between the centers of two adjacent square 
cells is given by the same expression. For this 


on reset xed lee 
case we have F=sin — sin 14 the 
arrangement of the cells are as shown in fig. 


—* and 
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Fig. 6. a) Schematic illustration of streamlines in a cross- 
section of a long roll convection cell. b) Regions of as- 
cending and descending motions in a square cell as given 
by the simple solution F(x, y). c) Ascending and descending 
regions of a physically more probable square cell. d) 
Ascending and descending regions of a hexagonal cell. 


6 b. It is seen that a unit cell ABCD as represent- 
ed by this expression is composed of a square 
central region and four equilateral triangles. 
Sincesl = 2, and AX = 2, h/{l = V2A,.A is.again 
given by (7.1). 

It may be remarked that even though the 
mathematical expressions for the square or 
rectangular cells represented in fig. 6 b are 
simple, this form of the cells seems to be 
dynamically inefficient, judging from the fact 
that the radial motion in the vertical planes 
OP and OQ must be very weak because the 
vertical velocity does not change sign. 


On the other hand, the arrangement in fig. 
6 c seems to represent a more efficient system. 
Here the counter motion takes place in the 
region between the two parallel squares whose 


sides are À/V2 and À. As is indicated in this 
figure, strong radial motion can take place in 
any meridional plane in this system. It is a little 
surprising to find that this cell form is not 
represented by any simple solution of eq. 
(4.1 a). 

The other regular equilateral polygons that 
can fill up the whole space between them are 
triangles and hexagons. The mathematical ex- 
pressions of these polygons have been obtained 
by CHRISTOPHERSON (1940) as solutions of 
eq. (4.1a). From fig. 6d it would seem that 
hexagon is a dynamically efficient configura- 
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tion because the cell wall is nearly at the same 
distance from the boundary of the Central 
Core which separates the ascending and the 
descending regions. It is seen that the boundary 
of this central core is practically circular and 
that the cell wall itself never depart very much 
from. a mean circle. 

The distance between two hexagonal cells is 


(7.2) 


where L is the length of a side of the hexagon. 

As to the question of the relative stability of 
the various plane forms in the presence of a 
mean current with vertical shear, the answer 
may be obtained by an investigation of the 
differential equation for nonviscous flow with- 
out Coriolis effect, which is given by 


(5 Ow 


A=V3 L=3.226 h 


Ow 
a NA =) - (U,+ oU,) L— + 


+ (L2+ g5,)v?w=o 


ox 
(7.3) 


It can be shown that for this case the long rolls 
(2/0x = 0) are the most preferred form of the 
perturbation. Detailed discussion of this prob- 
lem will be given in another paper. 

Now we shall turn to the cases when the 
stratification is unstable for ascending motion 
only. Our problem is to determine the 
dimension of the preferred convection cells 
and of the ascending region from the results 
obtained in the two preceding sections. For 
these convection cells the eftect of earth’s 
rotation is negligible, we shall therefore take 
100004097; 

In determining R,, and %,, it has been as- 
sumed that 6,4 is large enough such that a 
further increase of its value has little effect 
on &,. The results in tables 2a and 4 indicate 
that Baé= 2.0 may be taken as the proper 
value for long roll clouds and BAË = 1.5 
for circular cells. Using these values for Bag 
we can then determine the dimensions of the 
ascending region and that of a unit cell from 
the values of wé,. The results for these two 
cases are given in tables 6 and 7, together with 
the respective area ratios of the descending to 
the ascending regions. 

We note that while the area of the ascending 
region diminishes as / increases, the dimension 
of the cell as a whole becomes larger when the 
descending region is more stable. Thus for 
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1 = 5, the distance between two cloud rolls is 
about 9 h. Take 4 km as the effective depth? 
of the unstable layer, we then obtain 2 x, = 
36 km. On the other hand, if the depth is only 
about 1 km, the horizontal spacing will be 
less than 10 km. These values are comparable 
with the spacing of the long roll clouds 
obtained from satellite observations (see 
KUETTNER, 1958, WINSTON, 1960). However, 
more detailed comparison with observations 
awaits detailed measurements of both h and I. 


Table 6. Values of 2x,/h, A,/A,, and 2 x./h for long 
rolls. œm = 0.707, BAE 22.0. 


I I | 2) | 3 | 5 | Io 

ak, 0.7605} 0.450] 0.310] 0.190} 0.096 
2x,/h CON | Cezar || G28 1.081724 || Os 
2%#o/h 12.498 24.005 RS COIN ONE ST 
AA; 2.6 8.9 19.4 | 52.6 | 208.0 


Table 7. Values of r,/h, ro/h and A,/A,, for circular 
cells. œm = 0.707, BA& =1.5. 


l I | 2 | 3 | 5 | Io 
agı moe 10981004 NO. 
v1lh 0714 Mo SONO 40110 42110820 

2v2/h A| x | CA re 
A,|A, 2.8 |10.6 |25.3 |79.6 | 380.0 | 


8. Factors favorable to development of large 
scale motions 


As has been pointed out in section 6, the 
ascending area of the “limiting disturbance’ is 
of the same order of magnitude as the rain area 
of the tropical storms and that the distributions 
of the various quantities given by our solutions 
bear many resemblances to that observed in 
tropical cyclones. It is therefore very tempting 
to identify these tropical storms with the 
limiting perturbation and to attribute their 
origin to simple convection. However, since 
the most preferred scale of motion in an 
unstable atmosphere is that of cumulus cloud, 
corresponding to « = &, the motion must be 
dominated by cumulus cloud scale convection 
if it evolves from random small-amplitude 
perturbations. Therefore the initiations of the 


1 The effective depth is usually larger than the depth 
of the unstable layer because of the penetration of con- 
vection into the stable layers. 
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much larger and very distinctive disturbances 
such as easterly waves and tropical depressions 
must be attributed to other physical factors 
other than simple convective instability, such 
as barotropic or baroclinic instability of the 
mean wind, or shearing instability. The evi- 
dence that the central parts of many tropical 
disturbances are cold before they transform 
into hurricanes (see RIEHL, 1954) lends support 
to the supposition that these perturbations are 
barotropic in origin, therefore they are totally 
unrelated to the present study. 

However, some of these small amplitude 
traveling tropical disturbances do transform 
into deep tropical cyclones, and the energy 
required for such transformations are definitely 
supplied by the latent heat of condensation. In 
order that there shall be a gain of kinetic 
energy during ascent, it is evident that the 
stratification must be conditionally unstable. 
The question before us is: what are the con- 
ditions that must be satisfied in order these 
large scale convective systems shall develop 
and shall not be destroyed by the much stronger 
smaller scale systems which are most likely to 
develop in the unstable atmosphere? The 
following few points and conditions seem to be 
very significant in this respect: 


(i) A deep layer of moist but unsaturated, 
conditionally unstable atmosphere 
(Ym < y < ya), Which is stable with respect 
to small perturbations, but possesses 
enormous amount of energy that can be 
released by organized, finite amplitude 
perturbations. 


(ii) Rapid removal of moisture by condensa- 
tion and precipitation from the ascending 
air, so that the descending air is relatively 
dry, drier than the normal surface air. 


(iii) A warm ocean surface which supplies 
heat and moisture to the lower atmosphere 
continuously but slowly. The lower at- 
mosphere will remain unsaturated unless 
it has stayed a long time or traveled a long 
distance over the ocean. This condition 
will insure the growth of the large scale 
convective motion and avoid to have the 
cumulus scale convection from developing 
spontaneously. Furthermore, the small 
scale convection systems which do develop 
in the general ascending area of the large 
scale motion will cut off their ow nenergy 
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supply by their descending motion, there- 
fore will be short-lifed. 


(iv) Pre-existence of a large-scale distur- 
bance with weak convergence and diver- 
gence and weak vertical motion, and some 
kind of instability which initiates large 
scale convection. For example, when the 
absolute vorticity of the mean motion is 


ihn o) Se Pri Rei 
© 


definitely initiate large scale perturbations. 


negative, i.e., 


(v) Modification of the stratification in the 
general ascending region through upward 
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transports of heat and water vapor by 
both large and small scale convections 
makes the lapse-rate in this region to 
approach the saturated adiabatic gradually. 
After this has happened, the sporadic 
cumulus scale convection will stop growing 
too rapidly. On the other hand the large 
scale motion will be maintained by the 
horizontal temperature gradient set up by 
the motion. 


Since a quantitative study of the nature of 
the convective motion under these varying 
circumstances requires numerical integration, 
it must be left for further investigation. 
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Abstract 


The design and performance of simple atmospheric models is briefly reviewed as an introduc- 
tion to the question of energetic consistency. The two-layer model of Lorenz appears to be the 
simplest fully-consistent formulation, and the freedom of this model from the constraints upon 
the Coriolis parameter and static stability characteristic of the usual quasi-geostrophic models is felt 
to be particularly important. As a prelude to actual numerical integration, the baroclinic stability 
properties and the energy transformations of this two-layer model are examined. It is found that 
the perturbations are here significantly more unstable than when the static stability is constrained 
in the usual manner. For the finite-amplitude disturbances, it is found that some of their available 
potential energy is transferred to the mean state through the stabilizing action of the systematic 
rising of the warmer air. These results suggest that the variation of the static stability both pro- 
motes the growth of new disturbances and limits the growth of the mature disturbances—eftects 
which would appear to remedy some of the conspicuous performance errors of the conventional 


two-level models. 


I. Introduction and review of baroclinic 


modeling 


During the past decade a number of simple 
(two-level or two-parameter) models of the 
atmosphere were designed for the purposes of 
numerical weather prediction. Asis well known, 
these models are built around the familiar 
quasi-geostrophic approximation, and make 
use of a simplified vorticity equation for fric- 
tionless, hydrostatic and adiabatic flow (Extas- 
SEN, 1952; CHARNEY and PHILLIPS, 1953; 
THOMPSON, 1953; SAWYER and BUSHBY, 1953). 
Several of these models have been tested under 
a variety of conditions, and considerable ex- 
perience has been accumulated with respect to 


1 This research has been sponsored by the National 
Science Foundation under Grant NSF G-14079. 


their behavior with actual atmospheric initial 
data (BusHBy and HINDS, 1955; THOMPSON 
and GATES, 1956; STAFF, JNWP unIT, 1957). 
The principle characteristic errors attributable 
to the physical shortcomings of these models 
appear to be: (1) a tendency to underestimate 
the growth of relatively new large-scale dis- 
turbances, (2) to overestimate the growth of 
the relatively mature disturbances, and thereby 
to produce fictitiously large gradients, and (3) 
to display a spurious growth of anticyclonic 
circulation, especially in the subtropical lati- 
tudes. 

The complete diagnosis of these errors is a 
difficult and continuing task, although certain 
analyses strongly suggest the probable effects 
of certain of the common approximations 
made. The isobaric vorticity equation used in 
these formulations may be written 
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where {= k - V x Vy is the vertical compo- 
nent of the relative vorticity (k being the vertical 
unit vector and V the isobaric gradient opera- 
tor), f the Coriolis parameter, w= dp/dt the 
vertical motion relative to an isobaric surface, 
and Vy the horizontal wind. We note at once 
that the neglected terms representing the ver- 
tical advection of vorticity (wdf/dp), the so- 
called twisting term (k- Vo x 2 Vy /dp), and 
the divergence term ( - £2w/0p) may account 
for the systematic errors noted above. The 
contribution of the twisting term in particular 
has been estimated by REED and SANDERS 
(1953) and by ELIASSEN (1955), who find that 
it evidently plays a significant role in the main- 
tenance of the intense gradients associated with 
a frontal zone. For the larger-scale disturbances 
as viewed with the usual prediction grids, this 
term and the vertical vorticity advection appear 
to be less important than the other terms, 
although their effects often appear on a large 
scale (ARNASON and CARSTENSEN, 1959). The 


replacement of the divergence term — (¢ + f) ° 


dw/0p by — fdw/dp in (1) has also been noted 
as a probable cause of the models’ tendency to 
overdevelop the anticyclonic circulation and 
to distort its characteristic asymmetry relative 
to the cyclonic circulations (Starr, JNWP 
UNIT, 1957). 

In addition to these simplifications of the 
vorticity equation, the horizontal velocity 
field in these models is usually represented by 
the quasi-geostrophic relations 


Vi =f kx a 
Bolles ” 
where &= gz is the geopotential. This approxi- 
mation leads to a systematic overestimate of the 
vorticity in both cyclonic and anticyclonic 
circulations, especially in the lower latitudes; 
it is evidently responsible for much of the 
“spurious anticyclogenesis” observed with the 
geostrophic barotropic model (GATES, 1957; 
SHUMAN, 1957), and contributes to similar 
errors in the simple baroclinic models (CHar- 
NEY and PHILLIPS, 1953; THOMPSON and GATES, 
1956; SHUMAN, 1956; Putts, 1958). The 
approximation (2) and its consequent neglect 
of the vorticity advection by the divergent 
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wind has also been cited as responsible for 
errors of type (2) noted above (THompson, 
1956; PHILLIPS, 1958). 

The adiabatic equation accompanying (1) 
and (2) in the baroclinic models may be written 
in isobaric coordinates as 

oT 
a a a (3) 
where I” = 9T/op-«/c, is a measure of the 
static stability, & is the specific volume, and c, 
is the isobaric specific heat. Here the tempera- 
ture T is in turn related to the geopotential 
through the hydrostatic equation 
Ib 
oe +a&=0. (4) 

In the conventional two-level model the 
stability J” is assigned a constant value in both 
space and time. In view of the generally more 
unstable stratification in the vicinity of devel- 
oping cyclones as compared to anticyclones, 
this approximation may systematically over- 
estimate the effects of vertical motion on the 
temperature distribution, and hence underesti- 
mate the cyclonic development. As in the vor- 
ticity equation, use of the quasi-geostrophic 
approximation in (3) entails the neglect of the 
advection of temperature by the divergent 
portion of the horizontal wind, and this may 
lead to further systematic error. 


2. Consistency and baroclinic modeling 


The efforts to reduce the characteristic errors 
of the conventional two-level models have taken 
a number of directions, all more or less focused 
on the removal of constraints inherent in the 
system (1)—(3). As a first step the geostrophic 
approximation (2) may be replaced by the 
more general balance equation, 


, = th Py yp Py \? 7 
v2d-v-.(fvy) [5 - (>) ]- 
6. (s) 


In the barotropic model this step results in the 
virtual elimination of the anticyclogenesis 
errors of type (3) noted in À ı (SHUMAN, 1957). 
A similar although less spectacular improve- 
ment is also made by use of the balance equa- 
tion (or at least a non-divergent wind) in the 
simple baroclinic models (BUsHBy and HUCKLE, 
1956; PHILLIPS, 1958). 
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Further improvement in the models’ perform- 
ance might also be expected by the inclusion 
of the neglected vorticity and heat transfer 
mechanisms. In particular, the inclusion of the 
divergent-wind vorticity advection and the 
vertical vorticity advection in (1) has been pro- 
posed by THompson (1956), and the inclusion 
of the twisting term has been proposed on a 
number of occasions. Similarly, one might 
consider temperature advection by the diver- 
gent wind in (2), and allow the static stability 
to vary. 

Viewed as independent “corrections” to the 
system (1)—(3), these modifications lead to a 
large number of possible models. Due both to 
the relative difficulty of incorporating these 
effects and to the inconclusive results of the few 
such tests that have been made (CHARNEY and 
PHILLPs, 1953; CHARNEY, GILCHRIST, and 
SHUMAN, 1956), an extensive program of term- 
by-term generalization of the simple system 
(1)—(3) has apparently not been undertaken. 
This is perhaps fortunate in view of the consist- 
ency requirements to be placed upon the mod- 
els, as first discussed by HOLLMANN (1956) and 
in a more comprehensive fashion by LORENZ 
(1957). . | 

These consistency requirements may be most 
easily seen if we consider the complete vortic- 
ity equation for frictionless flow, 


a À 0 
Zr Vu" r(C+ ror Wis + 
+k-Vox a. (6) 


By rearrangement and use of standard iden- 
tities, this equation may be written in the form 


Gt Wet +e (+f) ex) + 


“ern CES TE 


270) 0 
= +A se (k-vox (k x vy)) + 
0 d 

+v.kx (Gkx vw) eZ | <0. (7) 


Here use has been made of the isobaric conti- 
nuity equation 


PAR RE (8) 
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where y is the poténtial function of the diver- 
gent portion of the horizontal wind, and y is 
the stream function of the rotational portion of 
the horizontal wind. The four groups of brack- 
eted terms in (7) represent the corresponding 
four terms in (6), with the last two members of 
the first bracketed group representing the vor- 
ticity advection by the divergent wind. 

If (7) is now integrated over the atmosphere’s 
total mass, assuming w = 0 at both the upper 
and lower limits of pressure, we find 


ILES (9) 


where ds is a fractional element of area on an iso- 
baric surface. Thus the total vorticity of the at- 
mosphere is conserved in the absence of friction 
(and with © = 0 at p = po), and (9) is conse- 
quently an integral requirement to be place 
upon simplified versions of the vorticity equa- 
tion if spurious sources and sinks of vorticity 
are to be avoided in the mean. It is evident from 
(7) that this consistency condition can be met 
in several ways. For example, if either vorticity 
advection by the divergent wind or the diver- 
gence term is neglected, then the other should 
be neglected also, or if only the former is neg- 
lected, the latter should be replaced by — fow/dp 
with f taken as a constant. (Another possible 
selection is considered in $ 3.) Further, the ver- 
tical vorticity advection and the twisting term 
must also either both be included or both neg- 
lected. Viewed in this manner the simplified 
vorticity equation (T) is consistent in the sense 
of (9) only when Vj; is considered non-diver- 
gent and when f = f, in - fdw/dp. Use of the 
geostrophic approximation (2) in this connec- 
tion is therefore a consistent procedure. In some 
numerical integrations f has not been consid- 
ered constant in the above manner (GaTEs, 
Pocınkı, and JENKINS, 1956). This question of 
consistency in simple baroclinic models has 
also recently been considered by’ ARNASON and 
CARSTENSEN (1959) and by Wın-NIELSEN 
(1959). 

We may apply similar considerations to the 
adiabatic equation (3), and for this purpose 
write it in the form 


IT d 
at J. 1)+ 9 (Tex) + (er) - 
oF 14/0 ao | _ 
» er to 
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Here the two groups of bracketed terms rep- 
resent the corresponding two terms of (3), 
with the last three members of the first brack- 
eted group representing the temperature advec- 
tion by the divergent wind. Integration of (ro) 
over the atmosphere’s mass as before yields 


[ford [foxdpds (11) 


The left-hand side here represents the variation 
of the total potential and internal energy of the 
atmosphere in hydrostatic equilibrium, and 
since total energy is conserved in adiabatic 
frictionless flow, we may regard the right-hand 
side of (r1) as the transformation of total poten- 
tial and internal energy into kinetic energy, 
accomplished by the familiar systematic rising 
of warm air (WHITE and SALTZMAN, 1956). 

We note that both temperature advection 
by the divergent wind and a variable static sta- 
bility are present in (10), and hence are implicit 
in the integral relation (11). If the former is 
neglected, say through the geostrophic approxi- 
mation, and the static stability 7 = 9T/dp - ac, 
is assigned a constant value in the manner of the 
conventional two-level model, then we obtain 
instead of (11) the result 


) 
IE Tdpds=o, 


which implies that the mean temperature of an 
isobaric surface cannot change. In this case, 
therefore, no net conversions between kinetic 
and potential energy (as here defined) are 
possible. The quasi-geostrophic models with 
constant static stability are thus artificially con- 
strained in the sense that during the develop- 
ment stage (wa < o) the local conversion of po- 
tential energy into kinetic energy may be syste- 
matically underestimated, while during the 
occlusion and decay stage (wa > 0) the potential 
energy conversion may be overestimated. In 
broad terms, one might interpret this constraint 
as contributing to the characteristic errors of 
types (1) and (2) discussed in § 1. 

It is, however, still possible to preserve the 
consistency condition (11) with a different set 
of simplifications. If one neglects the tempera- 
ture advection with the divergent wind, say by 
the use of (2), and in addition assigns the lapse 
rate JT/dp a constant value rather than the stat- 
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ic stability Z’,! then (11) results directly. Such 


a formulation might therefore serve as a con- 
sistent baroclinic model, although it still neg- 
lects possibly important processes. 


3. The Lorenz two-layer model 


A simple yet consistent model with a mini- 
mum of constraints has been designed by 
LORENZ (1957, 1960). This model is free of 
many of the limitations of the usual quasi- 
geostrophic formulations and may provide a 
significant step short of the primitive equations 
themselves. After a review of the two-layer 
Lorenz formulation, we shall consider its baro- 
clinic stability properties and energy transfor- 
mations as a prelude to actual numerical inte- 
grations. 

From (7) we may note that a consistent 
simplification is to write 

ov? 
a td Ce vy +f) + ¥- (Fv x) =0 (13) 


where V?y = & and where the variation of the 
Coriolis parameter has not been restricted. 
LORENZ (1957, 1960) has shown that consistent 
with this form of the vorticity equation are the 
adiabatic] and “balance” equations in the forms 


Ft (9. 6)+ 9-097) + 5 (8) = 


=0,; (14) 


vig-v-(fvy)=0, (15) 


where © is the potential temperature. Here (14) 
is equivalent to (10) without approximation, 
and (15) represents the “linear” portion of the 
balance equation (5), again without restriction 
of the Coriolis parameter; to this extent the 
system (13)—(15) may be termed a “non-geo- 
strophic” formulation. By applying (13) and 
(14), together with the continuity equation (8), 
to a two-layer atmosphere as illustrated in fig. 1, 
we may obtain the system 
y2 
2 2 ih yew + f) + J (t, v?r) = 0 (16) 


0 v7 


re Pre ver) - 
—v:(fvx)=0 (17) 


1 In the mean winter atmosphere over North America, 
the measure 2T/dp is nearly as constant in the troposphere 
as is the measure J’ (GATES, 1961). 
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0 (w =O) p=0 
YV+T— X —9+e p=Ry, 
2 (6) P=Ry, 
3 Re 020 p= Sry, 
4 =0) P=Pp, 


Figure I. The information levels and primary dependent 
variables of the Lorenz two-layer atmospheric model of 
variable static stability. 


D + J yO) +J (0) - v-(ov#)=0 (18) 


do 

> tJ RC) + J (t, 9) - vO-vy=o (19) 
BEE RON ET Vt) 0. (20) 
Here y = (y + Ws)/2,90 = (0, + O,)/2, 7 = 
(y — Ys) /2, o= (0: - O5)/2, and x= (ys - 41) /2 
where the subscripts refer to the levels shown 
in fig. 1; R is the specific gas constant for dry 
air, and x = R/c,. To obtain (t6)—(19) we 
have applied (13), (14), and (8) to both levels 
ı and 3, and evaluated vertical finite-differences 
over po/2 (Po = 1,000 mb assumed constant), 
with the boundary conditions © = o at both 
p = oO and p = po. We have also used the ap- 
proximation @, = (©, + @;)/2. Equation (20) 
is obtained by differentiating (15) with respect 
to pressure, using the hydrostatic relation (4), 
and evaluating 0w/dp by finite differences over 
Po/2 about the level 2. This is the system given 
by Lorenz (1957, 1960), except for a slight 
difference in the coefficient of V?@ in (20). The 
continuity equation accompanying (16)—(20) 

is 


20 2 
— - vy=0o 
FE CE 


(21) 


where w is the vertical motion at level 2. 
Equations (16) and (17) represent the “mean” 


AE M €. 
Se ee 
JO /dy 00/2} 
do dy 20,/2y 
re er 
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and “thermal” vorticity equations, respective- 
ly, while (18) and (19)-provide for the varia- 
tion of the mean temperature and static stabil- 
ity. Equation (20) represents a form of the 
thermal wind equation. Upon integration over 
the entire atmosphere (or other suitable closed 


0 
we obtain the results > v2y ds=0, 


2 A a N % 
fie [ 04-05, ods = 
r 


“ . od 5 Ls ’ 
Says [ Ov%zdsand 5 | v20 ds=0; 


region), 


this system is therefore consistent in the sense 
of (9) and (x1), and as suggested by Lorenz, 
probably constitutes the simplest such formu- 
lation. It is interesting to note that it is the 
variation of the static stability o in (19) which 
actually provides for the transformation of 
potential into kinetic energy in the mean. 


4. Stability properties of the linearized system 


In order to examine the stability of small 
amplitude perturbations in the baroclinic at- 
mosphere described by this model, we may 
partition the dependent variables into a basic 
state (denoted by the subscript zero) and a 
superposed perturbation (denoted by a prime) 
in the usual manner. For convenience, we shall 
consider the basic or undisturbed state to be a 
uniform non-divergent zonal flow; hence 
% = 0. Moreover, we shall regard Yo, To, Oo 
and oy as linear functions of y (distance north- 
ward) alone, and regard the perturbations as 
functions of x (distance eastward) and f (time) 
only. Although these assumptions somewhat 
oversimplify the problem, they probably do 
not conceal the basic stability properties of the 
present model as compared with those of the 
usual two-level model determined in the same 
manner. 

Assuming perturbations of the form exp 
[ik (x —ct)], where k = 2x/L is the wave num- 

er, L the wavelength, and c the wave speed, 
we find that the existence of non-trivial per- 
turbations in the linearized system requires 


0 0 0 
0 0 — ikf 

- (c=) =) huh Shee ies es) 
U* ze) fo) 

Reste O O 
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Here U = - dy4/dy is the mean zonal wind, 
U* = - 0t/dy is the vertical wind shear of the 
basic state, and 6 = 2f/2y. Upon requiring the 
zonal flow to satisfy the thermal wind equation 
(20), and d?U/dp? =0, we have the relations 


99, À ei 

(aaa (23) 
and 

Cine : > 

rn WP) fF OR (24) 


Expansion of (22) with the aid of these rela- 
tions leads after some manipulation to the 
cubic frequency equation 


(c — U)$+ m, (c- U)?+ m; (c- U) + m3=0, 


where 


(28) 


Here we have introduced the dimensionless 
ratios 
Rk?o, 
ll ae tif? > 
ne eu 
m 1 (30) 
26 


which depend critically upon the stability and 
shear of the basic state. The frequency equation 
(25) reduces to that of the usual two-level mod- 
el when the underlined terms in (22)—those 
dependent upon a variable static stability— 
are omitted. The familiar frequency equations 
for the barotropic and advective models may 
also be obtained by introduction of the ap- 
propriate simplifications. 

The frequency equation (25) is seen to depend 
upon the wave number of the disturbance and 
the parameters of the basic state in an involved 
manner; for the present purposes, however, 
only the critical conditions separating unstable 
and stable (or neutral) perturbations need to be 
determined, and we may proceed directly to 
this result by noting that the unstable waves 
correspond to the existence of complex conju- 
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gate roots of (25). This occurs only when the 
discriminant 


D=27m; - mim + mim; - 18 mı m; m; + 


(31) 


and thus the desired stability criterion is given 
by D =o. This condition, however, is still an 
awkward functional relation, and has therefore 
been determined numerically as follows. For 
each soo km wavelength to 8,000 km and for 
each 2 m sec! shear to 20 m sec-1, the value 
of D was computed and the curve D= 0 
positioned graphically in the stability diagram. 
Figures 2 and 3 show the results of these cal- 
culations for selected o,, and extend the pre- 
liminary results given earlier (GATES, 1961). 
These figures also show the growth rate of the 
unstable perturbations in the form of the time 
required for a doubling of the amplitude, cal- 
culated as above for the unstable waves from 


+ 4m} > 0, 


tp = (2%¢)-1LIn2 


Shear in m see —— 


Vertical 


0 — = L 4 lt 2 L L 


O | 2 3 4 5 6 Te 8 
Wave Length in 10° km —— 


Figure 2. The baroclinic stability properties of the two- 
layer model of variable static stability. The solid curves 
give the doubling time (in-days) of unstable disturbances, 
with the lower curve representing the critical stability 
conditions. Here the vertical shear is measured over po/4, 
the basic state’s static stability o, = 15 deg (measured over 
Pol4) at 45 N, and the basic zonal current U = 20m 
SEC 

The dashed curves are for the conventional two-level 

model of constant static stability. 


> 


Shear in m sec — 
Oo D 


Vertical 


(o) | 2 3 4 5 6 4 8 


Wave Length in 105 km — 


Figure 3. Same as fig. 2, but with og = 30 deg. 


where L is the wavelength and ¢; is the positive 
imaginary part of the phase speed.” 

The case o,=15 deg (fig. 2) is representative 
of the stability in winter for mid-latitudes, 
and shows the presence of a minimum unstable 
shear and a shortest unstable wave, both of 
which, however, are considerably less than 
those given by the usual two-level model 
(PHILLIPS, 1954). Disturbances as short as 2,000 
km wavelength may be unstable in the present 
model with o) = 15 deg (compared with 3,600 
km in the conventional model), and it is just 
for these shorter synoptic waves that the eftects 
of the self-determination of static stability 
could be expected to be the most important. 
For the longer waves, the critical stability con- 
ditions approach those of the usual two-level 
model. The growth rates of the unstable waves 


2 From the algebraic solution of the cubic equation (25) 
we find 


ci = 2-3/3 (a —b), where 


I 
a,b = ——{ (omy my - 2m} - 27m3) + [(om, m; - am? — 
3\2 


— 27M)” + 4 (3m; - my )®] a} 
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are here greater for all values of the shear than 
are those of the conventional model; evidently 
the variation of the static stability has a desta- 
bilizing effect upon the growing baroclinic 
disturbance, and is most pronounced for the 
shorter waves. For the wave lengths greater 
than 3,000 km the critical curve of fig. 2 is very 
close to that of the advective model (Sur- 
CLIFFE, 1951). 

For the larger static stability o,= 30 deg, rep- 
resentative of winter-time conditions at higher 
latitudes, fig. 3 shows an even greater de- 
parture from the stability criteria of the usual 
models. Here waves as short as 2,500 km may 
be unstable (compared with 5,100 km in the 
conventional case), and the minimum unstable 
shear is somewhat less than that of the usual two- 
level model. Comparison with fig. 2 shows that 
the unstable waves’ growth rates are here every- 
where less than those of the case o, = 15 deg, 
as would be expected. Finally, we may note 
that the growth rates of the longer waves are 
relatively insensitive to the static stability of 
the basic state, and the critical stability curves 
in the cases o = 15 and 30 deg are practically 
coincident for wave lengths greater than 
6,000 km. 

In comparing the present model’s stability 
characteristics with those of the usual 2-level 
model, one should also examine the reduction 
of (22) with d6,/dy =o. In this case the left- 
hand portion of the critical curve of figs. 2 
and 3 is shifted toward longer wavelengths 
by about 1000 km, with the right-hand portion 
remaining essentially intact. These new critical 
stability curves are the same as those for the 
usual 2-level model with halved static stability. 

These results suggest that the variation of 
static stability in the present model plays an 
important role in at least the initial growth 
of baroclinic disturbances, and the increased 
growth rates found for the shorter waves may 
help to reduce some of the systematic error of 
the conventional models noted in $ 1. 


5. Energy transformation functions 


Having shown that the stability character- 
istics of relatively young disturbances in the 
present model are significantly different than 
those of the conventional two-level model, it 
now becomes of interest to examine the role 
of variable static stability in the energy trans- 
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formations of the model. To this end let us 
first consider the kinetic energy per unit area, 


(32) 


and seck from (16) and (17) the behavior of 
the spatially-integrated kinetic energ gy for both 
the zonal average motion and for the disturb- 
ances. 


Writing () = () + (), where ( ) denotes 
the zonal average in the usual manner, we may 
average (16) and (17) and form equations in 
dylot, or/ot, Oy’ |dt, and dt" Jat. 

Multiplying these equations by Y, T, y’, and 


t’, respectively, and integrating over the earth’s 
take we find the walk 


— RK aoe Pa ——- \ 5 à a 
at 5 | Ras q fas Ÿy dy (9% Py + Tx T;) ds 


Sede a Vypr ve vt), 


0 I. 
Er Ty By rt Van) ds 
was food (33) 
20 
and 
[2 FEUER ee) 7 
en Kds= q fos (pe py + Tx Ty) ds 
d 
the (nem 
= O’ w’ ds, (34) 


where f Kas =f 2 tw. Vy+Vr.Vr)ds 


is the kinetic ae of the mean zonal (non- 
divergent) oo and 


[Ke 2 [iv y-vy+vrV-vr)ds 
a 


Rp 0 


2X F1 


noting that — =( / 04) =o for adiabatic mo- 


tion with © =o at both p = o and p =pp. 
Evaluating each time derivative in (3 8) from 


(18) and (19), we readily find 
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is the kinetic energy of the disturbances. Here 
use has been made of the zonal and disturbance 
forms of (8), (20), and (21). Equations (33) and 
(34) are the same as those found for the usual 
two-level quasi-geostrophic model by Pumps 
(19 56), and are not directly dependent upon the 
static stability or its possible variations. 

Let us now seek the corresponding equations 
for the available potential energy. As discussed 
by Lorenz (1957), the total available potential 
energy per unit area in the present model is 
given by 

fAd= 


Sides (Leds) + Sords- (J cas 
Jods+ (fads) m 


aeg 


(35) 


where (/ods)„ is the maximum (stable) stabil- 
ity which could be obtained by adiabatic re- 
stratification, and may be considered a constant. 
Here /'Ads is essentially a measure of the total 
variance of the potential temperature in each 
of the two layers of the model. By now writing 
0=9 +0 and 6 = 6 + 0’ 
manner, we have 


| Ads = 
a LEE = (/Ods)? + lois “sl 


in the usual 


2047|) [ods+ (fods), 

' (36) 
RI ese NO: 

er: a | 37) 


for the total available potential energy of the 
zonal flow and of the disturbances, respectively. 

Considering first the zonal case, we may 
differentiate (36) and obtain 


a Fhe |S a [aa - 2 ( [ ads) | | Sods-+ (fads), 4 


[76° ds - (/ ds)” + fords - (fads)? Yo + (f0ds)m | © Jo. 


(38) 


n a O'ds=2f O (Ov + 07, - 


- 0, -0 %) ds, 
Ak ds =2f 0, (ow tO Ca) Ye 
SX) ds-2fa(O' v2 x’) ds. 
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en 
sf ods=S 


vO-vy+v0"-.vy)ds, 
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writing (21) for both the zonal flow and disturb- 
ances, and using the relation 


(fods)a=-(fOds)?+ f (O)’ds+ f (o)ds, (39) 


we find after some manipulation the result 


d 
7 5, | Adm; 


fods - ve] Ve 
ia fou a (fods) 


2*g | sods + 


We may note that if we were to regard o as 
an absolute constant o = o,, then the last two 
terms of (40) would drop out, and we would 


have simply 
Ÿ Be 
oi |: Ads = 


2. Ov JE d'y, dir, | Gods (41) 


2196 


This reduction corresponds to the energy 
transformation obtained with the familiar two- 
level models (Paizirs, 1954 and 1956), when 
A is suitably defined. 

Turning now to the eddy available potential 
energy, we differentiate (37), giving 


Cal 
1/9» ads see oe =f Leas + (fods) À O'«' ds. 


R ; == en om 37 
tea [6 (Ow +07T,-0% 5 + Oy (o Wet OT; - OY ) |a . 


Sods - 


(40) 


De: Er [Sf D [(EYds+ 5 Sf (0) | 
at + Ae 7 ods) | 


_ Rpol /(O)?ds + S(o’)?ds] 0 
2*+1¢ | fods+ (fods) ml? 0 


xf ods. (42) 


From (18) and (19) we then find 


se ar le) pass 


-2 fO,(O'y, + ot, - 0 4.) ds 


= 2f o,(o' vx + OC, - O'y,)ds + pel Ow’ ds, 
0 


with which we may then rewrite (42) after 
some rearrangement as 


— % = Rpo 
ot AUS = 2 “g [fods+ ( Jods) Yl fle 


RK ees: EE = (Jods be 
+ Je w'ds = bes io 88 w ds ur © SRE (eds) 7 | oe O'w’ds. 


It is now convenient to collect the above re- 
sults in the compact form 


RS 2 AUS 
Gf Ris {KR} AR 


(44) 


(45) 


= Rds=— {RK} 4 td KR 


(Oy. +07, — 0'x,) + 0,(0'y + Or. - 


6°7,) |ds+ 


fods - 


(43) 


= i Ads=—{A-K}-{4-a}+{4’. A}, + 
+{A’- A}, (46) 

gf aw-taa-taeR}-{aah- 
- {AA}, (47) 
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where the energy transformation functions are 


given by 


ant = A 
{KK} = 0 uf E jy Wah + tat) + 


d 
+Ty (rip, + vst) | ds, (48) 
ars = Owds, (40) 


PE he 
far K'} si O'wds, (so) 
28 


a {| 6, (Ove 


2*g | fo ds+ (sods) | 


14:47} = 


OT — 043) +0, (pe t+ OT, - 


-0'x,)] ds, (s1) 


De LR leas (f0ds) al > > 
{4 > Ap, 2*g ; + el en (52) 


= 2R Sods _G ER 
MA AN — | Po al Cay ds. 
{ je ai Jods + Baie De 


The transformations (48)—(51) correspond to 
those found in the usual two-level model, al- 
though (51) is here modified by the stability 
variations. The transformations (52) and (53), 
however, do not occur in the conventional 
model, and represent effects directly due to the 
variations of static stability. We may note that 
these transformations modify the exchange of 
available potential energy between the zonal 
flow and the eddies by the same processes which 
are responsible for the exchange between po- 


tential and kinetic energy, i.e., the Ow and 


@'w’ correlations. 

In the middle and high latitudes of significant 
0%, the condition o >/ods will generally 

revail, since the predominantly lower static 
stability of the region below, say, 30 deg N 
will depress the areal mean stability /ods rela- 
tive to o in the remaining half hemisphere. 
Hence the transformation { A’. A}, will serve 
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in favor of that of the zonal mean, and in 
this sense acts as a brake upon the dis- 
turbances’ growth. Further, since fads < (/'ods) m 


generally, the transformation {A' . A}, will 
also act to increase the zonal available potential 
energy at the expense of the eddy available po- 


tential energy when JO ds<o, although this 
effect is probably somewhat less than that of 


{A’ . A}, and may act in the opposite direction. 


From (44)—(47) we may construct an energy 

iagram in the manner of PHicctps (1956), 
shown in fig. 4. Here the indicated direction of 
energy transfer is that characteristic of the 
growing mid-latitude disturbance. As in the 
conventional model, the growth of the disturb- 
ance kinetic energy depends upon the systemat- 
ic rising of the warmer air relative to the cool- 
er air. As noted above, however, this same 
process, through its effective increase of the 
static stability, will simultaneously reduce the 
supply of eddy available potential energy and 
hence inhibit excessive disturbance growth. 
The variation of the static stability thus serves 
both to encourage and then to control the de- 
velopment, and may ultimately suppress the 
development altogether at the end of the dis- 
turbance’s life cycle and thereby set the stage 
for new development. 


[Kk] 


Figure 4. A schematic energy transformation diagram for 
the two-level model of variable static stability. The indi- 
cated directions of transfer are those characteristic of 
growing disturbances in middle latitudes, with the trans- 
formations {A Ê K} and 4 : Ee, possibly proceeding 
in either direction. HereA denotes available potential ener- 
gy, K denotes kinetic energy, and the bar and prime 
denote the zonal average and the perturbations, respective- 
ly. The symbol [ ] denotes spatial integration. 
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6. Summary and conclusions 


The two-layer model designed by Lorenz 
(1957, 1960) appears to be the simplest baro- 
clinic model which possesses consistent relations 
between the total kinetic and available poten- 
tial energy, and between the available potential 
energy and the static stability. The model’s va- 
riation of static stability is perhaps its most 
significant generalization over conventional 
two-level models. The Coriolis parameter is 
allowed full variability, and this feature may 
help to reduce the spurious anticyclogenesis 
characteristic of the quasi-geostrophic formu- 
lations. 

The analysis of the baroclinic stability prop- 
erties of the linearized equations shows that 
the synoptic-scale perturbations are unstable at 
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approximately half the vertical wind shear re- 
quired for instability in the usual two-level 
model; this feature may help the present model 
to remedy the systematic underprediction of 
relatively young disturbances noted with the 
conventional model. The analysis of the energy 
transformations reveals a mechanism whereby 
the supply of available potential energy to the 
disturbances is regulated through the variation 
of the static stability; this feature may help the 
present model to remedy the systematic over- 
prediction of relatively mature disturbances 
noted with the conventional model. In order 
to test adequately these conjectures, a program 
of comparative numerical integration is now 
in progress, using northern hemisphere data 
from the International Geophysical Year. 
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Circulation and Energy Balance in a Tropical Monsoon 


By F. A. BERSON, C.S.LR.O. Division of Meteorological Physics, Aspendale, Victoria, Australia. 


(Manuscript received June 27, 1961) 


Abstract 


Mean monthly meridional and zonal circulations in the region of the Indonesia-Australia 
summer monsoon and associated transfers of heat energy are investigated using upper air sound- 
ings, rainfall and other data from two seasons. Daily sounding data are used to determine also 
the eddy flux of latent heat. A large increase of meridional circulation strength and energy 
transfer from November to December on the northern side of the monsoonal trough suggests 
that establishment of the regime is associated with a comparatively abrupt southward shift and 
intensification of the Northern Hadley cell in the relevant longitude sector. 

When the monsoon has become established latent heat relase in the region exceeds net export 
of sensible heat from it by a comparatively small amount. These quantities dominate the re- 
maining ones in the heat balance, viz. net radiative cooling in the troposphere and sensible heat 
transfer from the ocean. Items in the heat budget are compared with those obtaining during 
the pre-monsoonal regime in the area and on the winter side of the equatorial trough zone at 
large, the latter according to RIEHL and Ma kus (1959). 

Gain of sensible heat is also assessed separately for the lower and upper troposphere in terms 
of gain or loss of internal energy, the work done by horizontal pressure gradient forces and by 
vertical pressure gradient forces (gain of potential energy). The vertical upward flux of energy 
across mid-tropospheric level is deduced and the importance of buoyancy in bringing about 


this flux is indicated. 


I. Introduction 


According to recent compilations of the 
earth-atmosphere radiation budget (cf. J. Lon- 
DON, 1957) poleward transport of heat across 
10° N latitude in late winter amounts to some 
60 per cent of the maximum poleward heat 
flux at 30°N which is 10.5: 1014 cal sec-t. 
These figures illustrate the important role of 
the equatorial trough zone (ETZ) as a major 
source or distributor of heat energy for balanc- 
ing losses in the remainder, more than three- 
quarters, of the winter hemisphere. 

The outflow of sensible heat from the tropo- 
sphere overlying the ETZ occurs chiefly in 
the diverging upper limb of the mean meridio- 
nal (toroidal) mass circulation associated with 
the Hadley cells in the two hemispheres. This 


export of heat occurs as flow of enthalpy and 
of potential energy the sum of which is nearly 
equal to ‘potential heat’ (cf. Batt, 1956) and 
will be referred to as such for ease of style. 
Since radiation losses in the troposphere are 
considerable, the heat budget in the ETZ 
requires an accordingly large rate of latent 
heat release in precipitation. It has long been 
known and demonstrated on detailed obser- 
vations (RIEHL ET AL., 1950, RIEHL and MALKUS, 
1957) that the adjacent tradewind regimes act 
as formidable accumulator of latent heat. In a 
later paper (1959) the latter authors assessed 
latent heat transfer by mean meridional mass 
flow and eddy flux across the winter boundaries 
of the ETZ—i.e. at 10 degrees latitude from 
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the trough—to 6.6 in the above units whereas 
transfer of sensible heat in the lower half of 
the troposphere amounted to less than half 
this value. 

Net import of latent heat into the ETZ 
implies within an excess of precipitation over 
evaporation from the underlying ocean 
surface. The average rate of upward flux of 
total heat energy across the soo mb level 
within the zone must evidently be very high 
as it has to account for the accumulations 
below that level of latent and sensible heat due 
to lateral import, evaporation and sensible 
heat flux from the underlying surface. Allow- 
ing for a reduction due to net radiative flux 
divergence in the lower troposphere the re- 
quired rate of upward flux is about 17 units. 

The authors demonstrated that such rate of 
upward flux cannot be nearly accomplished 
by mean mass ascent and/or vertical diffusion. 
They proposed a mechanism by which power- 
ful updrafts in the protected cores of some 
1,500 to 3,000 large thunderstorm cells would 
carry heat energy to great heights at a rate 
exceeding the required average rate of transfer, 
the compensating somewhat slower and broad- 
er downdrafts enabling dryer air to make 
contact with the surface energy source. The 
hypothesis has as yet to be confirmed from 
special observations and measurements. 

There is the parallel problem of large-scale 
inhomogeneity in longitude, of meridional 
flow, temperature, pressure height and hu- 
midity mixing ratio. This would give rise to 
standing-eddy flux of energy. The cited 
authors have shown that ocean-—continent 
differences in heat energy are in themselves 
rather small, but no attempt has been made to 
evaluate the differences in meridional circula- 
tion strength that occur between longitude 
sectors. Evidently standing-eddy flux of latent 
heat cannot be disregarded in the moisture 
balance on a global scale (Starr and PEIXOTO, 
1958). 

It would seem appropriate to treat mon- 
soonal circulation regimes—the most promi- 
nent among seasonal standing-eddies—as in- 
tegral parts of the tropical circulation. The 
summer monsoons especially deserve attention 
in view of high average rainfall, the year-to- 
year variations of which have such a strong 
impact on the national economies of the coun- 
tries concerned. Thus studies have been made 
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of transfers of heat and zonal momentum 
associated with the Indian monsoons (Rao, 
1959); and of transfer and seasonal balance of 
zonal angular momentum in the region 
dominated by the Indonesia-Australia sum- 
mer monsoon (BERSON and Troup, 1961). 
The objective of the present paper is a quan- 
titative assessment of the circulations and 
associated transfers of heat energy as well as 
heat balance in the latter region. Changes 
from pre-monsoon to established monsoon 
regime and differences with respect to the 
winter ETZ will also be considered. 


Notation 


V horizontal velocity vector 
u, v, w west-east, south-north 
velocity components 


and vertical 


Be sea level pressure; p, constant pressure 
close to tropopause level (~ 100 mb) 
S area of horizontal cross section, do 


element of S 
B boundary of S; dl line element of B 
EB = GT + Agz + La, total heat energy 
E, =¢,T + Agz, ‘potential heat’ approxi- 
mately 
vertically integrated net radiative flux 
divergence in the troposphere 
net radiation on ocean surface 
upward flux of sensible heat from 
ocean surface and in the ocean respec- 
tively 
€ evaporation; P precipitation 
A monthly mean value of A; A’ departure 
of daily from monthly mean 


departure of A from pressure mean, 
Op 


Other symbols above and in the 
following have their usual meaning 


unless explained. 


2. The heat balance equation 


The change of heat energy in unit time 
averaged over a month, in a volume deter- 
mined by S and po - pı assuming hydrostatic 
equilibrium, is given by 
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where the transfer of heat energy from ocean 
to atmosphere H + Le must satisfy the en- 
ergy balance equation at the surface 


R,=H-G+Le (2) 


Equation (1) in which bars have been omitted 
for convenience over all variables except E 
and V, is correct under the assumptions that 

(i) Eddy transfer of energy across the isobaric 
surface p, is negligibly small in comparison 
with transfers at the ocean surface. 

(ii) Horizontal density gradients in the trop- 
osphere are sufficiently small to omit a term 


containing V v, 0 which otherwise would occur 
iy EGP i): 

Owing to the absence of any but weak 
baroclinity in the tropical atmosphere over- 
lying an ocean region assumption (ii) is very 
much justified. Assumption (i) cannot be a 
priori justified but it should not lead to 
discrepancies of a magnitude affecting the re- 
sults and their interpretation. 

The term B represents a change of energy in 
the volume by a change of mass contained in 
it. It is proportional to the area-integrated 
change of sea-level pressure since no mass can 
be transferred across the constant pressure sur- 
face py. 

The term C is the negative horizontal flux 
divergence in the troposphere and represents 
therefore lateral import to, or export of heat 
energy out of the volume according to whether 
it is positive or negative. 

In the expression for flux divergence, the 
contribution by mean meridional motions 
may be separated into two parts, viz. 


Be ASB ERS ON) 


and similarly for-zonal transport. Here À, is 
the transport due to a.net mean mass flux and 
A, is a toroidal transport (PRIESTLEY, 1949), 


which is determined by the covariances of E 


and v or E and u respectively. The asterisk on 
either energy or. velocity component may be 
omitted in the integrand. Introducing eq. (3) 
into (r) it can be shown that the term B cancels 
out with the term containing À, (v) and A, (u). 
Eq. (1) thus simplifies to 


Er FE 
le hae 
Söp 


= I div, (V* E+ WE) dpdo+ 


I 
= 
+ or (H+ Le - ÔR;) do (4) 


where V’ E’ is the residual synoptic-scale eddy 
term. Its evaluation strictly requires to deduce 
daily means from several soundings each day. 
The transfer term would then not include 
mezo-scale eddy flux such as is produced by a 
seabreeze circulation (see section 4). 

By Gauss’ theorem the first integral on the 
right of eq. (4) is equal to the outflow of energy 
through the vertical ‘wall’ bounding the vol- 
ume. The quantities VE, V,'E’ in which 
V, is the outward normal wind component, 
can be determined from zonal and meridio- 
nal transports and the orientation of the sides 
of the selected quadrangle. The implications 
of interpolation between sounding stations will 
be discussed in section 4. 


3. Energy and wind profiles 


The toroidal fluxes of total and potential 
heat are uniquely determined by the vertical 


profiles of u, v and E, or Ej respectively. 
Typical profiles are shown in Figure 1 for 
Singapore and Darwin. Very similar profiles 
of E and E; result for longitude-averaged 
seasonal means at 15° N (PALMÉN, RIEHL and 
VUORELA, 1958) and on the winter boundaries 
of the ETZ (Rieu and Markus, 1959). Heat 
energy decreases upwards above a shallow 
neutral layer and increases above the 700 mb 
level, first moderately, then more sharply in 
the uppermost troposphere. (The sharpest 
increase just below the 100 mb level is a 
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Figure I. Distribution of total heat energy (E), potential 
heat (Eg) and meridional wind component over Singa- 
pore and Darwin, January 1957. 


further complication which may be disregarded 
here.) The increase in the upper troposphere 
of Ex reflects traverses of markedly subadia- 
batically stratified layers by the sounding 
balloons in a great majority of all soundings. 
This apparently persistent stable layer is 
present in the tropics over both oceans and 
land (cf. cited authors; Rao, 1960). It is seen 
that total heat energy at 150 mb exceeds the 
value at the surface where the actual source is 
located. However, energy must be transported 
upwards to great heights because the excess 
amount in the equatorial zone is carried away, 
poleward as potential heat, effectively by 
meridional mass flow in the uppermost 
troposphere as will be seen from the profile of 
meridional velocity. The consequences of the 

resence, in the average profile, of a stable 
ce in the upper troposphere for an effective 
mechanism of energy transfer to great heights 
has been mentioned in the Introduction. The 
compilation of energy balance in the monsoo- 
nal regime (section 8) will give an opportunity 
to assess the position in such a regime. 

A glance at the curve for Ey (at Singapore) 
shows that apart from the mentioned anomaly 
just below the tropopause this can be taken as 
a straight line, to a first approximation. Toroi- 
dal transport of potential heat can therefore be 
determined as a function of vertical wind 
profile. In the analysis of local (station) transfers 
of potential heat it is of advantage to deter- 
mine the regression 
Tellus XIII (1961), 4 


21 = Ja Pr at b S (5) 
where S is the strength of the circulation, 
Sy =Vmax — Vmin being defined as meridional 


circulation strength and S,= Umax — min as 
relative zonal circulation strength. The co- 
efficient b in the regression equation depends 
according to the foregoing consideration on 
the mean slope of the E,-curve in a diagram 
such as Fig. 1. In eq. (5) the subscripts max., 
min. denote the primary maximum and 
minimum in the profiles of the wind com- 
ponents (S,>o in the Northern, S,<o 
in the Southern hemisphere Hadley-type 
circulation). In the subtropics S, is presumably 
closely related to the meridional circulation 
strength in the manner defined by PRIESTLEY and 
Troup (1954). S, is in the nature of a vertical 
shear of zonal wind. Regressions according to 
eq. (5) will be evaluated in section 5. 


4. The region and data considered 


Selection of an oceanic area dominated by 
the circulation associated with the Indonesia- 
Australia northwest monsoon, yet at the same 
time suitable for the calculation of energy 
flux divergence, is dictated by the location of 
the sounding stations. During parts of Decem- 
ber and as a rule during the whole of January 
and February, the regime of the westerly 
monsoon extends from Malaya across the 
Timor, Arafura and Banda Seas to the coastal 
belt of northern Australia, and westward in the 
Indian Ocean to about 90° E longitude (Fig. 2). 

At or near the northern fringe is Singapore, 
near the southern fringe Cocos Island on the 
edge of the Southeast tradewinds, Darwin, 
and Port Hedland on the Australian coast, the 
latter station in the orbit of the circulation 
around the shallow heat low over the western 
parts of the continent. Detailed calculations of 
energy transports and flux divergence were 
made from the data of these four stations. 
Meridional transports and circulation param- 
eters were also evaluated for Koror in the 
Pacific Ocean to the east of Mindanao (Philip- 
pines), in about the longitude of Darwin; and 
for six additional sounding stations in the 
Australian region north of 33° S (Fig. 2 a). 

The periods used for all these stations, 
except Koror, are November to February 
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Figure 2. (a) Normal M.S.L. pressure field for January. 
Position of intertropical convergence in monsoon regime 
indicated by heavy dotted line. Arrows are vector means 
of surface wind in five degree latitude-longitude squares. 
Stations used: @ rawinsonde, © rainfall. 
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(b) Streamlines and isotachs of average 200 mb flow in 
February (After TROUP, 1961 a). 


1956—57 and 1957—58, the latter season during 
IGY. The data used in the cross section showing 
meridional circulation during the December to 
February season (Fig. 3) vary from station to 
station, but none cover less than four seasons. 

The area selected for the assessment of heat 
balance is indicated on Fig. 2. In this region 
with the exception of the extreme northern 
portion, November is a pre-monsoonal month. 
During December marked often abrupt changes 
occur in the circulation and in rainfall. During 
January the monsoonal regime is usually es- 
tablished in the region but is nevertheless ir- 
regular and fluctuating (HANNAY, 1945; TROUP 
1961 a). 

Considerations of moisture balance require 
an evaluation of average precipitation in the 


FA: BERSO!N 


area. This is based on fifteen stations within or 
close to the boundaries of the area, as indicated 
on Fig. 2 a. Radiation data are taken from 
various sources in the literature but mainly 
from Lonpon (1957). However, for an esti- 
mation or radiative flux divergence in the tropo- 
sphere according to a formula by MÖLLER 
(1960) cloud amount and dewpoint depres- 
sion at 300 mb are required. For the former 
recent climatological atlases have been consult- 
ed, the latter has been evaluated from the basic 
sounding data. 

As to the representativeness of monthly 
means evaluated from soundings carried out 
once a day mention was made in section 2 of 
the desirability to exclude seabreeze effects. In 
the event coastal stations figure prominently, 
but the bias for seabreeze components in the 
lowest 3,000 ft is not thought to be serious 
except during the first season (noon soundings) 
at Port Hedland and possibly Darwin. In the 
final assessment of heat balance during February 
(Table 4) both seasons have been averaged so 
that a systematic error would be reduced to 
negligible size. 

Such error could arise also from linear inter- 
polation in connection with the evaluation of 
energy flux divergence from transports at the 
four mentioned stations. However, as seen in 
Fig. 3, the meridional component in the upper 
troposphere changes across the relevant region 
from positive values in the north to negative 


Pressure (mb) 


Port 
Hedland Cocos |. 


Darwin Lae Singapore | Yap Manila 
Koror 

Figure 3. Meridional (toroidal) wind component in 
vertical section from Australian Heat Low across mon- 
soonal trough, December to February average. Speed in 
knots. 
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in the south, presumably in a quasi-linear 
manner. Similarly, in the lowest layer— 
separated from the upper limb of the circula- 
tion by an ‘inactive layer’ (cf. PALMÉN, RIEHL 
and VUORELA, 1958)—the northerly compo- 
nent decreases southwards from the equator to 
about 20° S. 

The zonal flow is strongest in the upper 
troposphere. Here the easterly ‘jet stream’ at 
about 50,000 ft which generally in the region 
is located somewhat north of the equator, 
shifts to a few degrees south of it with the start 
of the monsoon season. Average streamlines 
and isotachs at the 40,000 ft level for February 
according to Troup (1961 a) are shown in 
Fig. 2 b. The maximum easterly component 
at this level remains well over the northern- 
most portion of the area in question. The 
presence of stronger westerlies close to the 
ITC between, say, Darwin and Singapore 
(Fig. 2) may however involve errors in the 
layer surface to 700 mb. 

For completion of monthly series IGY cards 
were used. In both seasons, but especially inthe 
pre-IGY season, soundings were not made on 
single days or during sequences of days. The 
computation of eddy fluxes requires in such 
cases careful interpolation. In the event in- 
terpolation was obviated where possibly by 
treating slightly shorter periods. Processing of 


daily wind components was facilitated by the 
use of CSIRAC Electronic computer. 


5. Transports 
(a) Meridional energy transports 


Toroidal and eddy transports of total heat 
energy at stations between 10° N and 35° S are 
set out in Table 1. To compare these with the 
rate of energy flow across latitude circles as 
required by radiation unbalance on a global 
scale (last column), February has here been se- 
lected; for heat storage in the atmosphere and 
ocean may then be regarded as negligibly 
small. Values north of the equator were de- 
duced from Lonpon’s (1957) radiation bud- 
get, values south of the equator by judicious 
interpolation between these and the distribu- 
tion calculated poleward of 45° S by GABITES 
(1960). 

Table r should be taken as giving merely a 
broad survey of meridional heat transfer in the 
general wider region marked by continental 
cooling in the north and continental heating 
in the south. The northward heat flux at sta- 
tions near the northern edge of the monsoonal 
trough is consistently stronger than required 
by global radiation unbalance. At Cocos 
Island and Darwin, near the southern end of 
the trough, fluxes are southward of consider- 


Table 1. Meridional heat transport during February, northward counted as positive. Unit 10° cal cm! 
sec-!. Toroidal transport refers to total heat energy, eddy flux to latent heat only (except in footnoted 


value). 
Toroidal Eddy Total | Radiation 
Station = a Av. 57—| Unbal- 
3 1957 | 1958 1957 1958 58 ance 

Tropics 
ECOROM EAS tes EN 7.4 N 134.5 E 4.9 1.4 
Singapore. ..... Tea! 100.9 1.0 AT 0.4 0.2 ZEIT 1.0 
Eoeosılsland.....|t 12.09 103.9 — 1.5 — 14 — 0.4 — 0.6 — 1.9 SET, 
Dai ae re 12.4 130.8 — 1.4 — 1.6 — 0.2 — 0.8 — 2.0 
Port Hedland...| 20.3 118.6 — 0.6 EL 0.4 — 0.1 — 0.7 2 
CIONCUEE EEE 20.6 1200 —AT — 4.1 
Subtropics 
Perth 31.9 115.9 110 357 
HOTTES ep 30.9 128.1 — 2.0 0.3 AE 
Lord Howe Isl..| 31.5 159.1 — 0.4 — 11 baat 
Rathmines .....| 33.1 15007 — 0.7 — 1.7 


1 Eddy flux of enthalpy and moisture for Auckland (N.Z.) for period Jan.—Feb., 1950—52, un- 


published data (Priestley & Troup, 1954). 


eee een 0. 
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Figure 4. Toroidal flux of potential heat as function of 
circulation strength (a) Meridional, (b) Zonal. Sloping 
lines indicate regressions according to eqs. $ a and 5 b. 


able magnitude whereas the global radiation 
budget at this time of year does not require 
meridional transfer of appreciable magnitude 
across the relevant latitude circles. Further 
south at Port Hedland, poleward transport 
differs but little from the radiation value. There 
is thus a suggestion of increased meridional 
heat flux divergence in the monsoonal longi- 
tude sector as compared with the belt. 
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Cloncurry in Queensland, Rathmines near 
Newcastle and Lord: Howe Island in the 
Tasman Sea show poleward transports in 
excess of that required by average radiation 
values. If this were interpreted as an additional 
effect of continental heating, it would be 
consistent with the figures for Perth on the 
west coast which indicate a marked equator- 
ward energy transport. 


(b) Circulation strength and potential heat trans- 
port 


In Fig. 4 a and b potential heat transports 
are shown as function of the corresponding 
circulation strengths. As the meridional cir- 
culation reverses from Northern to Southern 
Hadley sense across the axis of the monsoonal 
trough (cf. Fig. 3), the net transports reverse 
from northward to southward. The regression 
equation (5) becomes 


1.040778, Seed 


Ar (v, E,) ie 1.6+ 0.76 Sy S, <8 


(sa) 


where the constants and S, are in the units of 
Fig. 4. Correlation coefficients r (A,, S,) have 
the values 0.95 and 0.91 in the positive and 
negative range respectively. 

The ETZ position in Fig. 4 a refers to the 
values of S, and A, derived from the profile 
of meridional flow and from potential heat 
flux on the winter side at a distance ten degrees 
latitude from the position of the trough during 
the December—February and June—August 
seasons. The flux value is therefore the average 
of northward flux coupled with N-Hadley cir- 
culation at 5° N and southward flux coupled 
with S-Hadley circulation at 2° N. The Singa- 
pore (1.4° N) and Koror (7.4° N) plots for 
November, ice. before the establishment of 
winter monsoon, show agreement in order of 
magnitude with average hemispheric-scale 
toroidal circulation strength and heat flux. 
However, from November to December there 
takes place at Singapore a sharp increase of 
circulation and flux which by February have 
four times the magnitude of the ETZ winter 
values. 

The stations south of 10° S show predomi- 
nantly S-Hadley circulation of moderate 
strength in agreement with the long-term av- 
erage during the December—February period 
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(Fig. 3). However, in November and December 
rather small negative values obtain at Darwin 
and Port Hedland. These reflect low-level 
northerly components associated with sea- 
breezes, an effect of sampling bias which was 
discussed in section 4. This effect too is respon- 
sible for an inconsistency of low-level com- 
ponents as between Figures 2 a and 3. 

The regression equation for zonal transport 
is: 
See 0948,48, > © 


FACE nor SS, <0:6b) 


with coefficients 0.96 and 0.77. Since the high- 
level maximum casterlies are located just south 
of the equator, both northern and southern 
station plots are found in the positive and 
negative range of S,. Absence of noticeable 
trend from November through February at 
Singapore reflects the influence of the semi- 
annual wave in the equatorial upper easterlies 
(cf. CLARKSON, 1955; TROUP, 1961 a). Note- 
worthy is the unexpectedly small scatter about 
the regression line in the range of moderate and 
large positive S,, occurring mainly in N ovember 
and December well south of the equator, where 
deep easterlies decrease with height or shallow 
easterlies are overlain by westerlies. 


(c) Sensible and latent heat transport 


It is informative to plot meridional transport 
of potential heat (hencetorth more suitably 
referred to as sensible heat) against latent heat 
transport in a scatter diagram (Fig. 5). Through- 
out the period November to February at all 
stations except Port Hedland, latent heat is 
transported into the monsoonal trough from 
north and south while sensible heat flux 
diverges, the export of the latter exceeding 
invariably import of the former kind of energy. 
At Singapore meridional energy exchange 
increases abruptly from pre-monsoon to mon- 
soon regime, while the values for Koror 
indicate a similar but more gradual change. The 
data for Cocos Island and Darwin indicate 
reverse changes although they are less abrupt 
and lack consistency. At Port Hedland little 
change is apparent in the meridional energy 
exchange which remains low throughout. 

But for the abruptness and magnitude noted 
at Singapore, these changes would be consistent 
with a gradual seasonal southward shift of the 
Hadley circulations in either hemisphere. 

Tellus XIII (1961), 4 


O 


Northward 


Sensible heat transfer (10° cal cm’ sec”) 


Southward 


Southward Northward 


Latent heat transfer (10° cal. cm.’ sec:') 


Figure 5. Meridional (toroidal) transport of sensible heat 

versus transport of latent heat. Sloping lines indicate zero 

net transport by meridional mass circulation. A. Mon- 

soonal regime. B. Pre-monsoonal or transitional (Koror 
in winter Trades). 


A similar analysis of zonal transports is not 
feasible on account of the limited sampling. 
It has been eastablished in recent studies 
(PALMER, 1958; BERSON and Troup, 1961) 
that frictional control on time averaged large- 
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Figure 6. Components of integrated energy flux diver- 

gence in monsoonal regime. @ meridional, x zonal flux 
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lines see text.) 
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scale zonal fow—i.e. above a shallow plan- 
etary friction layer—is very weak even in the 
tropics to within a few degrees latitude from 
the equator. ‘Uncontrolled’ shifts in longi- 
tude of zonal patterns are indeed not un- 
common. A large and widespread anomaly in 
high-tropospheric zonal flow in the tropics 
occurred e.g. in January 1958 and has been 
analysed in some detail by Troup (1961 b). It 
is reflected in Fig. 6 as an exceptionally large 
magnitude of zonal divergence of both latent 
and sensible heat flux over the portion of the 
monsoonal trough that now will be considered 
in some more detail. 


6. Flux divergence of sensible and latent heat 


The representation in Fig. 6 is similar to 
that in Fig. 5. Export or import of total heat 
by either components of flux divergence can 
be inferred from the appropriate displacement 
of the plots from the sloping line. Consistent 
export by meridional flux divergence is readily 
verified. On the average over the two sampling 
seasons the excess of sensible heat export over 
latent heat import has the same order of 
magnitude as the equivalent net export from 
the winter end of the ETZ!. This agreement 
proves, if nothing else, the adequacy of the 
station network for the method of calculating 
flux divergence. 

Noteworthy is the small overlap, in this 
representation, between meridional and zonal 
flux divergence: latent heat import into the 
monsoonal trough coupled with sensible heat 
export out of it for meridional divergence and 
predominantly vice versa for zonal divergence. 
With the establishment of monsoon and an 
accompanying marked increase in the meridio- 
nal transfer of sensible and latent heat energies, 
the overlap vanishes altogether. Furthermore, 
during the monsoon period, average zonal 
flux divergence is considerably less in magni- 
tude than meridional flux divergence. Since 
toroidal transports do not represent residual 
terms between large quantities (see Fig. 1) and 
therefore are subject to errors no larger than 


compound product errors of v into T, or v 
into q, etc., one must accept as real this lack of 


! Rate of outflow from a sector of the relevant 10° 
latitude belt equal in size to the trial calculation area 
(Fig. 2), so that effectively heating rates per unit area 
are compared. 


FAR ER SION 


compensation between meridional and zonal 
flux divergence with respect to the two forms 
of energy. 

However, net divergence of total heat flux 
emerges as a relatively small residual term be- 
tween meridional and zonal flux divergences, 
similarly as is generally the case with mass 
flux. Although this makes it difficult to assess 
accurately the heat budget, it will be of little 
consequence for determining the broad or- 
ganisation of lateral and vertical energy ex- 


change. 


7. Eddy flux 


According to calculations by PALMEN et al. 
(1958) poleward eddy flux of latent heat 
across 15° N in winter is 2 - 1014 cal sec}, i.e. 
by an order of magnitude smaller than the 
transfer of total energy by mean meridional 
circulations across that latitude. Riehl and 
Malkus computed the moisture eddy flux as a 
residual term in the ETZ heat balance. They 
found that it is directed outward (poleward) 
on the winter side of the ETZ and represents 
slightly more than the net export of total heat 
by mean meridional circulations. 

The monsoonal regime of Indonesia and 
Australia is marked by not infrequent disturb- 
ances such as tropical cyclones, ‘rain depres- 
sions’ and surge-like movements of the ITC. A 
direct calculation of the moisture eddy flux 
seemed therefore desirable. Some results are 
set out in Table 2. The vector fluxes during 
monsoon period indicate moisture export 
from the region. The resultant fluxes at Singa- 
pore, Darwin and Port Hedland are consistent 
with horizontal moisture gradients in accord- 
ance with the large-scale distribution of land 
and sea, but the southward flux at Cocos Island 
appears unduly large. 

As seen from the last row in Table 2 there is 
an increase in the magnitude of eddy transports 
from November to December which is mainly 
a transitional period. Furthermore, during 
monsoon regime eddy fluxes are considerably 
stronger on the southern end of the monsoonal 
trough than on its northern side. This would 
reflect the positions and intensities of the 
synoptic-scale disturbances. 

The effect of including moisture eddy flux 
in the calculation of flux divergence is shown 
in Fig. 6 where horizontal lines indicate ap- 
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Table 2. Eddy-flux of latent heat during 1956/57 and 1957/58 seasons, northward and eastward counted 
positive. Unit 10* cal cm”! sec 1. (‘“Monsoon’ refers here to January and February averages.) 


Saeed | Meridional Zonal Vector 
Noy. Dec. | Monsoon Nov. | Dec; | Monsoon pasos 
Singapore... - — I.I — 1.5 2.0 — I.I | — 3.7 — 2.2 3 310% 
Cocos Island....| 2.5 — 5.8 — 9.0 23 —LI 3.7 10 160° 
Darwin Ser. — 0.2 — 19.2 — 3.7 1.6 8. 7.9 9 120% 
Port Hedland...| —2.5 — 3.1 — 3.8 — 2.1 143 — 4.5 6 230° 
Average of | | | 
Modulus...... 1.6 7.4 AO 1.8 se 4.6 _- | — 


propriate additions or subtractions in excess 
of 0.03 : 1014 cal sec-1. The effect is seemingly 
of some consequence in December and January 
and is to decrease meridional flux convergence 
(import) or to increase zonal flux divergence 
(export) as calculated from mass circulations 
alone. 

Eddy flux of enthalpy has been found to be 
completely negligible at 15° N during winter 
according to calculations by PALMÉN et al. 
(1958). For stations near the equator Rao 
(1960) showed that it is not always negligible 
compared with moisture eddy flux. The 
magnitude of the eddy flux of enthalpy and its 
divergence have not been investigated here; 
however to ignore this eddy flux and its 
divergence is consistent with assumption (ii) 
in section 2. 


8. Sensible heat balance 
(a) General considerations 


Evaluation of the heat storage term in eq. 
(4) showed it to be of smaller order of magni- 
tude than integrated flux divergence. In view 
of this and the minor role of eddies in the 
transfer of heat energy, the case may be treated 
as a quasi-stationary heat-driven mass circula- 
tion. As such it is sustained by surface-sensible 
heat source and by release of heat of conden- 
sation (precipitation heat, LP). The energy 
used in the work done by vertical pressure 
gradient forces against gravity as the air 
ascends in the monsoonal trough zone (“gain 
of potential energy”), is partly off-set by loss 
of internal energy due to expansion and net 
radiative cooling. It is, however, not a priori 
clear in which manner transfers and balance of 
energy are arranged as between the lower 
troposphere, i.e. the monsoon layer proper, 
and the overlying part of the troposphere; 
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and in which layer mainly work is done by 
horizontal pressure gradient forces at the 
boundaries by the surrounding atmosphere on 
the system (the integral of V,po-1=R V,T). 

The work by eddy stresses at the lateral 
boundaries and by viscous stresses at the sur- 
face—the latter tacitly omitted in eq. (4)— 
represent additional losses of energy. The 
former effect is thought to be very small. As 
to the latter a check calculation applying the 
average value of surface zonal stress in the 
region (cf. Berson and Troup, 1961) and 
surface zonal wind components (Fig. 2) showed 
that the frictional stress term is less than 1 per 
cent of horizontal energy flux divergence in 
the troposphere. In the following both these 
works term will be omitted. 

In order to transform eq. (4) into a balance 
equation for sensible heat only, we eliminate 
evaporation heat and latent heat flux conver- 
gence with the aid of the moisture balance 
equation for the troposphere 

s 


L {&@-9 d= = GP ll aaa (7) 


Omitting eddy flux of enthalpy we thus 


obtain 


S Po S 
L [Pdo-, $ | ra [ de (8) 
pı 


where - Q,=H+6R, is the sensible heat 
source for the troposphere (actually a sink), i.e. 
exclusive of conversion of latent into sensible 
heat. 


(b) Formulation for two layers 


The ‘monsoon layer’ (I) may be represented 
by the layer p)- 550 mb, the 550 mb level 
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being but little higher than the equatorward 
ends of the moist tradewind layer. As evident 
from Fig. 1 the moisture content is not quite 
negligible above 550 mb, even on the average 
over a month. In disregarding release of latent 
heat in the upper layer (II) we shall not so 
much introduce an inaccuracy as rather con- 
sider the upper boundary of the monsoonal 
layer as not strictly specified. Release of latent 
heat of fusion in the layer II is negligible and 
the process of fusion is of no consequence for 
the balance except if it is systematically bound 
up with melting in the layer I. This is not 
typical for the tropics (cf. RIEHL and Markus, 
1959). The balance equations for sensible heat 
for the layers I and II may be written in the 
form 


S Po 
far-) 18 | Tr Bjdp dis 
co) ne 
S 


+ ap (ew E),: do 2) 


Ss p 
Jorude=: PART. TE 
R Pi 

- (ev E), do 


where p’ is the pressure at the top of the mon- 
soon layer. 

Consistently with the definition of p’, in the 
expression of the vertical (upward) transfer of 
energy E may be approximated by c, T + 
A gz. 

The average upward transport of potential 
energy Agz at the level p’ represents a 
portion of the total work by vertical pressure 
gradient forces in the monsoon layer, since 


and 


(9 b) 


= = 6 Se TR 
gf evde-glew2o-¢ f 20 = dz, (10) 
o Oz 


by parts. On account of the continuity equation 
the integral on the right, after integrating 
once more over the area S, can be and com- 
monly is, interpreted as outflow of potential 
energy through the vertical surface bounding 


S. Integrating from surface to 2, since (ow), = 


aA aD aS ON 


o by assumption, the first term on the right of 
eq. (10) vanishes and the work by vertical 
pressure gradient forces is equal to the outflow 
of potential energy. 

With respect to each of the layers, I and II, 
the first term on the right of eq. (10) is finite 
and must occur with opposite sign in equations 
(9 a) and (9 b). This is true also for the upward 
transport of internal energy through the level 
p’. In Table 4 the first term on the right in eq. 
(10) figures as “vertical flux divergence of po- 
tential energy”, the second term as “horizontal 
flux divergence of potential energy”. 


(c) Analysis 


Results of the computations on equations 
(7), (8) and (9) are summarized in Tables 3 and 4. 

As will be seen from Table 3 monthly rain- 
fall amounts vary considerably from one 
month to another, even within the monsoon 
period. The variation of precipitation heat is 
however broadly similar to the variation of net 
sensible heat export from the monsoonal 
trough zone. With the exception of December 
1957, the difference between them is positive, 


Table 3. Tropospheric Heat budget. Gains and 
Losses in units 10-* ly min !. For symbols see text. 


ee Sens. Establishment 
Period | Raink. Prec. | Heat of Monsoon 
Heat | Ex- | (cf. Troup 1961 
mm 
port a, b) 


Monsoon onset 
early in month 
Fully established 


regime 
Dec. 57 200 | 27.2 | 27.5 | Onset in N. 
Australia late 
December 
Jan. 58 145 | 19.4 | 12.8 | Anomalous high- 
level zonal 
winds 
Feb. 58 195 | 28.6 | 23.3 | Fully established 
regime 
Average DO OLIS 
Monsoon Os ÖRa A 
2 12.0 7.8 
Winter BLOSS SEITZ 
BAZ 8.4 13.5 5.1 
Pre- 8.3 — 4.2 
Monsoon DS 
(Nov. 57) 
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‘Table 4. Two-layer sensible heat budget for Indo- 
nesia-Australia summer monsoon (February). Unit 
10-2 ly min-!, 


Flux Divergence 


Heat 
Hori- | Ver- Source 
zontal] tical 202 
| LAYER I 
| Internal 
| Energy... 3.4] 25.5] 28.9)-ÖRa| — 6.4 
| Hor. pressure 
HOrces - x 1.4 TA 5.2 
Potential 
| Energy....|— 16.9] 17.4 0:51 LP 31.6 
= ial 42.9 30.8| | 30.4 
I LAYER II 
| Internal 
Energy... .|— 44.0|— 25.5|— 69.5) ÔRa | — 7.1 
Hor. pressure 
Forces... .|— 17.8 — 17.8 
Potential 
Energy.... 99.1|—17.4| 81.7| LP o 
re 5.6| | — 7.I 
TROPO- 
SPHERE 
Net loss and 
SO 25.2 23.3 


and its average value is slightly more than 15 
per cent of precipitation heat. The position 
appears to be significantly different from 
winter ETZ as there sensible heat export 
amounts to only 60 per cent of precipitation 
heat. The data in Table 3 for a pre-monsoonal 
regime indicate that with markedly less rainfall 
in the region horizontal heat transfer can even 
reverse, since in this particular month (Novem- 
ber 1957) import obtains. This is not confirmed 
by November 1956 for which available rain- 
fall data were however indadequate for an 
assessment of P. 

If as suggested in Table 3 during the mon- 
soon regime by and large, Q, is positive and 
small, if not negligible, in magnitude com- 
pared with precipitation heat, then according 
to eq. (8) net radiative cooling would somewhat 
exceed surface heat transfer. According to 
London’s tables the summer value of ÔR, in 
the belt o—20° N is 13.5 units. It is based on 
average cloudiness and vertical moisture distri- 
bution in the belt. Applying Möller’s formula 
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(see section 4) to the monsoonal data, a value 
of 12 obtains. With this estimate H can be 
evaluated from eq. (8) and turns out to be 8 
units. 

Evaporation heat has here been determined 
for individual months from P and resultant 
moisture flux convergence according to eq. 
(7), the average value being 14 units and 
differing little from Riehl & Malkus’ evaluation 
by means of an aerodynamic bulk formula. 
Accepting H and Le as representative, the 
Bowen ratio would be 0.55 and so considerably 
higher than the ETZ value (0.4) obtained from 
net radiation and evaporation according to 
eq. (2) assuming G= o.! If for consistency 
London’s value for R, (19 units) is taken here 
too, then from eq. (2) one finds G = — 3. In 
view of the uncertainty of the precise magni- 
tude of the other terms in eq. (2) no signif- 
icance can be attached to this result except 
that it indicates a correct order of magnitude 
of ocean heat flux G (MALKUS, 1960). Estimates 
of Ro, € and H from various other sources 
such as ALBRECHT’S (1940) measurements in the 


Java Sea and Budyko’s Atlas of Heat Balance 


would not help to clarify the position; for a 
ten per cent error in P would change H bya 
large fractional amount. 


Detailed heat budget for monsoon and overlying 
layer.—The above considerations give some 
assurance that the following analysis will be in 
order, with tolerance of the above discrepan- 
cies. Table 4 summarizes the two-layer budgets 
for February selected as only suitable sampling 
period on the basis of Table 3. The right-hand 
side of Table 4 serves again as a check on the 
values for net flux divergence, within the 
limits of the stated tolerance. H has here been 
obtained from the calculated evaporation 
amounts with the assumption that the Bowen 
ratio is 0.4. The discrepancy between the left 
and right sides of the Table is 2 units, i.e. for 
the whole troposphere. 

As might have been anticipated the lower 
layer gains sensible heat at a high rate while the 
upper layer loses at a much lower rate, the net 
gain being 25 units. To understand how these 
gains and losses are brought about, one must 


1 For oceans in general a range of Bowen ratio 0.1 to 0.2 
is held to be valid but the higher value in the waters of 
the equatorial trough zone is supposed to be consistent 
with the downdraft hypothesis mentioned in Introduction. 
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consider more closely the work by vertical 
pressure gradient forces (gain of potential 
energy). As will be seen in layer I the relevant 
term is negligibly small. This would indicate 
that the effect of buoyant forces is there large 
even when averaged over the whole area under 
consideration; for we may write 


fi ewde=* [wir f edo fetta 
(11) 


where 0** is the departure of density from the 
area mean value. Where 0** is negative 
buoyancy is positive and so the second term 
on the right must be large negative in the 
monsoon layer to reduce the term on the left 
side to a small residual. The discriminating 
effect of buoyancy in effecting vertical heat 
transfer introduced by PRIESTLEY and SWINBANK 
(1947) has been proposed also by Riehl and 
Malkus in the mentioned updraft hypothesis. As 
corollary the downdraft mechanism is quali- 
tatively verified by the consideration of large 
negative buoyancy forces required to explain 
the very high gain of potential energy in the 
upper layer (II), viz. 82 units. 

A relatively large negative value of the 
horizontal pressure gradient work term in the 
layer II implies that work is done by the 
environment on the system. This is consistent 
with a meridional equatorward drift (amount- 
ing to about 1 m sec~) the presence of which 
was found necessary to maintain zonal angular 
momentum balance in the monsoon-dominat- 
ed portion of the ETZ in the eastern hemisphere 
(Berson and Trovp, 1961); for the horizontal 
pressure gradient south of the high-level 
maximum equatorial easterlies is directed al- 
most due northward and is much stronger 
than gradients in the lower troposphere. 


Comparison with non-monsoonal regime.—A 
glance at Table 5 will show that in regard to 
the equatorial regimes organization of transfers 
and sensible heat gain is very much alike: the 
lower branch of the circulation is self-main- 
taining and in addition provides more than 
the energy necessary to sustain the circulation 
in the overlying troposphere where the net 
loss of sensible heat is incurred. The remaining 
excess energy is transferred in the upper limb 
to regimes outside the equatorial trough zone, 
i.e. in the case of the monsoon regime, to 


Er A. BERSON 


Table 5. Sensible heat gain in equatorial and tropi- 
cal regimes. Unit 10-1 ly min !. For last row see 


text. 
Southern Winter a 
| 
Monsoon 18494 en 
Layer (mb) 1,000| 550 —| I,000 500) 1,000 
— 550] 100 500 125 700 
Internal 
Energy 2.9 7 0, 15 930 0.3 
Hor-pressuzel 20.2. |7 1.815. 9.02 20:6 0.1 
Forces 
Potential 
Energy 0.1 8.2| 0.4 3.0| — 0.2 
St | 27 1.9 | u 0.2 
Upward 
flow of 
Energy 4-3 2.3 
w* 10.6 6.5 km/day 


north- and southward with rather indefinite 
net zonal transfer (cf. Fig. 6). Excess energy is 
there made available to be consumed in reduc- 
ing radiative cooling or by release of potential 
energy during descent. 

As to the magnitude of the net heat gain, in 
the monsoonal regime this is almost twice as 
large as in the winter ETZ. The corresponding 
factor in turnover of energy is however in 
excess of two. 

The last two rows serve to indicate the 
magnitude of the buoyancy effect discussed 
earlier. For this can also be gauged by the 
enormous value for the rate average vertical 
mass flow by which the energy transfer across 
mid-tropospheric level could be accomplished 
(w* = f owEdo[/ o Edo). 

The comparison with heat gain in the trade- 
wind regime is illuminating. According to the 
analysis by Rien and Markus (1957) the 
North Pacific tradewind layer during spring 
gains sensible heat at the rate 0.2 in the units of 
Table 5, or only about 5 and 10 per cent of the 
gain in the lower troposphere of the equatorial 
regimes. In the Trades potential energy is lost 
as the air descends through the top of the 
tradewind inversion. As mentioned earlier, 
the role of the tradewind system in the general 
circulation is that it acts as an accumulator of 
latent heat which is made available to the 
equatorial trough zone, but the rate of supply 
to the monsoon regime is greatly stepped up 
on account of the enhanced mass circulation. 
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Conclusion 


It has been shown that the mean mass 
circulations associated with the monsoon are 
the determining factor in the horizontal trans- 
fer of heat energy. The effect of synoptic-scale 
eddies is to counteract the mean-flux conver- 
gence of latent heat there being evidence that 
this is mainly due to ocean-land moisture 
gradients. The effect is much smaller than on the 
winter side of the equatorial trough zone at 
large where it represents a substantial item in 
the heat budget, according to the study by 
Riehl and Malkus. This difference may not be 
a real one since in the latter case eddy transfer 
of moisture has been deduced as a residual 
term in the heat balance. 

The check on quasi-balance between the 
rates of net sensible heat export and latent 
heat conversion by an estimate of radiative 
sources and sinks was of necessity a crude one 
but served nevertheless a useful purpose. Thus 
there is some assurance on the magnitude of 
vertical heat transfer at mid-tropospheric level 


and the necessity of invoking buoyancy effects 
for its mechanism. The problem of vertical 
heat exchange in an equatorial regime has 
here been formulated in a somewhat different 
way than earlier by those authors and shown 
up on mainly independent data. It remains a 
problem to be solved both theoretically and by 
means of observational techniques as yet not 
elaborated. 
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Abstract 


The TIROS II meteorological satellite, containing two television cameras and a family of 
electro-magnetic radiation experiments, was placed into orbit on November 23, 1960. The 
medium resolution radiometer is a cluster of five sensors, their optical axes inclined at 45° to 
the spin axis of the satellite. The spin of TIROS provides the scanning motion. Channels are 
sensitive to the following spectral bands: (1) 6 to 6.5, (2) 8 to 12, (3) 0.2 to 6, (4) 8 to 30, and 


0.55 to 0.75 microns. 


Introduction 


TIROS II, the meteorological satellite launch- 
ed last year on November 23 (1960), carried 
two television cameras and various detectors to 
sense infrared emission and reflected solar 
radiation from our planet. The systems aspect 
of the radiation experiments, the design of the 
radiometers, and the physical significance of 
the measurements have been discussed pre- 
viously (HANEL and Stroup, 1960; BANDEEN 
ET AL, 1961., ASTHEIMER, DE WAARD and JACK- 
SON, 1961; HANEL and Wark, 1961). The 
quoted publications also contain results in 
form of analog presentations of individual scan 
lines, of sample data points which show the 
apparent blackbody temperatures of certain 
geographical locations, and of maps where 
radiation emerging from the 8 to 12 micron 
atmospheric window has been plotted. Most 
of these data have been reduced by analog tech- 
niques. The main bulk of the information 
which is contained in more than 600 reels of 
magnetic tape, cach 5,000 feet long, has to be 
processed by computers. Hand analyses of 
samples must precede automatic data reduc- 
tion mainly to judge the validity of the experi- 
ment and to confirm proper functioning of the 
instruments. We are now at the point where 
we have gained confidence in the measure- 


ments, and where the computer programs are 
working. The five radiation maps shown in 
this paper are pilot runs of the automatic data 
reduction system. 


The Five Channel Radiation Experiment 


To give a better understanding of the mean- 
ing of the maps presented, the instruments 
which collected the information shall be 
discussed briefly. The medium resolution radi- 
ometer scans earth and outer space as TIROS 
spins about its axis. Each detector measures 
the radiation flux within its five degree instan- 
taneous field of view and within the spectral 
range determined by filters and other optical 
elements. The nominal bandwidths of the 
five channels are: | 


1) 6 to 6.5 microns, water vapor absorption | 

2) 8 to 12 microns, atmospheric window 

3) 0.2 to 6 microns, reflected solar radiation 

4) 8 to 30 microns, terrestrial radiation 

5) 0.55 to 0.75 microns, response of TV 
system | 


The optical axis of the instrument is inclined 
at 45° to the spin axis of TIROS. The radi- 
ometer responds to the intensity difference 
between two diametrically opposite directions. 
This configuration gives maximum coverage of 
the earth from a spin stabilized satellite and uses 
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Fig. 1. Radiation map 
from Tiros II scanning 
radiometer: Channel 1 
(6.0—6.5 y) Orbit 30, 
Nov. 25, 1960, 1238 
GMT to 1250 GMT. 
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Fig. 2. Radiation map 
from Tiros IT scanning 
radiometer: Channel 2 
(8—12 u) Orbit 30, 
Nov. 25, 1960, 1238 
GMT to 1250 GMT. 
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Fig. 3. Radiation map 
from Tiros II scanning 
radiometer: Channel 3 
(0.2—6.0 u) Orbit 30, 
Nov. 25, 1960, 1238 
GMT to 1250 GMT. 
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Fig. 4. Radiation map 
from Tiros II scanning 
radiometer: Channel 4 
(7.3—30 u) Orbit 30, 
Nov. 25, 1960, 1238 
GMT to 1250 GMT. 


Fig. 5. Radiation map 
from Tiros II scanning 
radiometer: Channel 5 
(0.55—0.75 u) Orbit 
30, Nov. 25, 1960, 1238 
GMT to 1250 GMT. 


THE TIROS II RADIATION EXPERIMENT 


Fig. 6. 
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Surface weather chart for 1200 GMT Nov. 25, 1960 showing Orbit 30 


of Tiros II. 


outer space as a radiation reference. The out- 
put of the detectors is stored on magnetic tape 
in the satellite and transmitted to earth once 
each orbit. The receiving stations at Pt. Mugu 
and Ft. Monmouth interrogated the TIROS II 
radiation experiment close to 2,000 times 
successfully. 


Data Reduction 


Each ground station sends the magnetic 
tape to the Aeronomy and Meteorology Divi- 
sion, Goddard Space Flight Center, in Green- 
belt, Maryland. The analog information is then 


screened and translated into digital form 
suitable for the IBM 7090 computer of the 
Meteorological Satellite Laboratory, United 
States Weather Bureau. The computer output, 
also on magnetic tape and in digital form, 
gives radiant emittance in watts per square 
meter, position of the satellite, geographical 
locaton of the area observed, viewing angle, 
sun angle, and Greenwich mean time for al 
data points. The maps shown in this paper have 
been hand drawn from computer printouts. 
Eventually mapping will be done autom- 
atically by contour plotters. 
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Fig. 7. Nephanalysis and surface weather chart for 1200 GMT Nov. 25, 1960, showing 
Orbit 30 of Tiros II. 
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Calibration 


The maps show radiant emittance W in 
watts per square meter as seen through the 
instrument’s filter. 


W = | Wipid? 


The spectral sensitivity of a channel is called 
pa and W; is the spectral radiant emittance of 
the area within the field of view. Ideally #: 
should be constant within the band specified 
and zero at all other wavelengths. Actual 
filters used approach this condition. Spectral 
sensitivities have been discussed in detail pre- 
viously (BANDEEN ET AL., 1961; HANEL and 
Wark, 1961). During calibration channels 1, 
2 and 4 were exposed to a blackbody of known 


emittance W,(T). Relations between W and 
the blackbody temperatures are shown in 


Table ı. 


Table 1. 
Blackbody W Watts per square meter 
Temperature for channel 

in deg. K I 2 4 
180 0.028 SR 9-5 
200 0.10 6.5 15 
220 0.28 12 23 
240 0.67 20 35 
260 1.4 32 49 
280 2.6 47 67 
300 45 66 89 


The interpretation of the results of channels 3 
and s has to take spectral sensitivities into 
account also. À perfectly diffuse surface irra- 
diated by one solar constant illuminates the 
detectors of channels 3 and $ with 1,070 and 
250 watts per square meter respectively. All 
scales show higher accuracy at higher radiation 
levels. Channels 2 and 4 exhibit a better signal 
to noise ratio then channel 1, and channel 3 
data are superior to data from channel 3. 
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The radiation maps shown cover the Atlantic 
Ocean from Newfoundland to the Iberian 
Peninsula and North Africa. A surface pressure 
map and a cloud cover analysis based mainly 
on ship reports are shown in the same scale on 
Figs. 6 and 7. . 

The three thermal channels, Figs. 1, 2 and 4, 
show a general temperature gradient from left 
to right which is caused by latitudinal and 
diurnal effects. The local time over Labrador 
is 9 A. M. and over France about 1 P. M. 
Superimposed over this gradient are synoptic 
phenomena. The relatively clear area west of 
Spain can be identified as a warm area on 
channel 2 and 4. In contrast to this, dense clouds 
associated with frontal system (e.g., 35° W 
45°N) appear cold in thermal channels 1, 2 and 
4 and bright in channels 3 and 5. 

As expected, bright areas in channels 3 and 5 
which indicate solid overcast are very well 
correlated to low temperatures in channel 2 
which measures radiation from the atmospheric 
window. This correlation will be used to 
extend cloud analysis to the nighttime side 
of the Earth, unobserved by channels 3, 5 and 
television. 

All maps shown display data uncorrected 
for viewing angle and solar zenith angle. Low 
temperatures at the southern and northern 
boundary of the scan area in channels 1 and 4 
are caused by limb darkening. 

Calculations of the outgoing radiation based 
on temperature and humidity profiles obtained 
from balloon soundings and measurements 
using a TIROS radiometer on a balloon will 
provide additional information on the accuracy 
of the scales obtained in preflight calibrations. 
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On the Geostrophic Flow at the Surface of the Pacific Ocean 


with Respect to the 1,000-decibar Surface 
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Abstract 


The geopotential anomaly at the surface of the Pacific Ocean with respect to the 1,000-decibar 
surface is computed. Horizontal variations in geopotential at greater depth are shown to be 
smaller and are neglected. The geostrophic flow at the sea surface with respect to the 1,000- 
decibar surface is shown to agree in general with the averages of currents estimated from the 
set and drift of vessels and tabulated in various atlases. The major gyres of the ocean are clearly 
shown and one new feature, a South Equatorial Countercurrent, is revealed in the western 
Pacific. 

Since there is meridional flow in the ocean, the flow is not entirely geostrophic and the 
contours of geopotential anomaly may lie at some angle to the streamlines. This is most evident 
where the contours intersect the coast and cross the equator. A qualitative consideration of wind 
stress indicates that where wind and current are in the same direction geopotential anomaly will 
rise along a streamline and where they are opposed geopotential anomaly will decrease along a 
streamline. This effect can at least qualitatively account for the rise of geopotential anomaly to 
the west in the region of the trades and toward the east in the northern hemisphere westerlies. 


Introduction 


For several decades oceanographers have 
been making some use of the geostrophic rela- 
tion to estimate flow in the ocean. Although a 
few systematic attempts have been made to 
compare the geostrophic flow with other 
current measurements, they have generally 
been confined to small areas and short periods. 
Wiist’s (1924) comparison of geostrophic flow 
in 1914 in the Florida Straits with direct 
measurements in the period 1885 through 1889 
is the classic example in the literature. It deals 
with a very fast flow in a narrow, shallow 
passage with a sharply sloping reference surface. 
Although the agreement is excellent, it does 
not follow that a useful reference surface can 
always be chosen or that relative velocities in 
the slower moving deeper ocean will always 
be so accurate. 
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Of the various charts and studies of the 
geostrophic current in the Pacific, the only 
chart covering the greater part of the area 
was that prepared from the CARNEGIE data 
(FLEMING AND OTHERS, 1945), but the data 
were too sparse to permit a detailed study of 
the circulation. Other large areas were studied 
by Upa (1955), Mao AND YosHIDA (1955), 
Dog (1955) and DODIMEAD (1958). 

In this paper it seemed worthwhile to com- 
pare the geostrophic flow of a substantial part 
of the ocean with the long-term averages of 
set and drift of vessels. The results are encour- 
aging, and they imply that in certain areas 
where other data are lacking (Bering Sea, 
central South Pacific) the geostrophic currents 
can be accepted with some confidence. In 
particular, the indication of a geostrophic 
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South Equatorial Countercurrentinthe western 
Pacific is of interest. Since the rest of the 
equatorial current system is found in the 
geostrophic system, this may be accepted as a 
real current. 


I. The geostrophic approximation 


Although the variation of the geopotential 
difference between two pressure surfaces may 
be interpreted in various ways, such as the 
slope of the sea surface under wind stress or a 
simple comparison of the average sea level 
and density of different areas and oceans, its 
most common interpretation is in terms of 
the geostrophic current. 

If certain terms are neglected, the equations 
of motion may be written 


adp/ax= fvt+adty/dz (ra) 
adp|/dy = - fu + adr,[dz (1b) 
adp/dz= -g (1c) 


where x, y, and z are directed toward the east, 
the north, and upward, u and v are the velocity 
components toward east and north, « is specific 
volume, f is the Coriolis parameter 2 © sin & 
(w the earth’s angular velocity, ¢ the latitude), 
p the pressure, and +, and +, are the vertical 
shearing stresses. 

The geostrophic approximation neglects the 
friction terms in equations (1a) and (1b). The 
anomaly of the geopotential difference between 
the 1,000-decibar surface and the sea surface 
is determined from (1c) and shown in Figure 1. 
The gradient of the anomaly is then used to 
estimate the geostrophic flow for comparison 
with various atlases (Section IV). It is well to 
point out that although the neglected terms are 
usually small compared with those retained 
and (as will be seen) the geostrophic current 
may be at any place approximately equal to 
the actual current, the neglect imposes certain 
important limitations which lead to more or 
less obvious contradictions. 


2. The data 


The data used in Figure 1 are listed in the 
appendix. Among these were the three co- 
operative expeditions NORPAC (June—Octo- 
ber 1955), EQUAPAC (July—November 
1956) and EASTROPIC (September—Decem- 
ber 1955), each of which attempted to obtain 
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an approximately isochronal sampling of the 
hydrography of the upper 1,000 meters over 
large areas of the Pacific. 

Most of the northern observations were 
made in northern summer and the southern 
observations in. southern summer, but there 
probably remain some artificial features on 
the chart which are the result of seasonal and 
nonseasonal variations of the hydrographic 
regime. 


3. The reference surface 


The 1,000-decibar surface was chosen as a 
reference because most of the NORPAC, 
EQUAPAC and EASTROPIC measurements 
did not reach much deeper, and these were 
the most suitable materials for an isochronal 
treatment. Since all of the results are to be 
referred to the 1,000-decibar surface, it is 
necessary to describe at least approximately its 
variation with respect to some deeper surface. 
From the CARNEGIE, DANA, DISCOV- 
ERY, TRANSPAC, SNELLIUS, NORPAC, 
VITYAZ 25 and 26, OB 3, DOWNWIND, 
MUKLUK, and CHINOOK materials a chart 
of the geopotential differences between the 
1,000- and 2,000-decibar surface was prepared. 
Its details will not be discussed here other 
than to say that it had a systematic variation 
somewhat similar to the shallower chart, but 
that the range of variation was much less. 
North of 50° S the total range is about 0.22 
dynamic meters, and over the whole Pacific 
the range is about 0.50. The range of the 
geopotential difference in the upper thousand 
decibars north of 50° S is 1.40, or six times as 
great, and over the whole Pacific is 2.00, or 
four times as great. With these limitations in 
mind the anomaly of geopotential difference 
between the sea surface and the 1,000-decibar 
surface will be discussed as if it represented 
north of s0°S the actual shape of the sea 
surface. 


4. The surface currents 


A. Direction 


The conventional descriptions of surface 
currents are derived from measurements of set 
and drift of ships compiled into atlases and 
from versions of these that have been smoothed 
and elaborated by various authors, in particular 
SCHOTT (1935). 
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Fig. 1. The anomaly of geopotential distance betweenthe zero- and the 1,000-decibar surfaces in the Pacific Ocean, 
in dynamic meters. 


The major surface currents of the Pacific 
Ocean as given in the atlases are all indicated 
by the field of geopotential anomaly. The 
northern subtropical anticyclone, and its com- 

onents (the Kuroshio Current, the west wind 
drift, the California Current, and the North 
Equatorial Current) are all indicated by the 
appropriate direction of the contours. 

Tellus XIII (1961), 4 


The southern subtropical anticyclone, which 
consists of the East Australia Current, part of the 
west wind drift, the Peru Current, and the South 
Equatorial Current, is also quite clearly defined. 
The geographical variation in slope of the sea 
surface is about the same as in the northern gyre 
except that across the East Australia Current 
it is not so great as across the Kuroshio. 
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Along the coast of South America there is a 
well-defined equatorward eastern boundary 
current. To the south there is some evidence 
of a poleward flow analogous to the poleward 
flow along Canada and Alaska. The data do not 
give an explicit answer, and the Atlas of Pilot 
Charts (U.S.N.H.O., 1950) gives only mar- 
ginal evidence of such a poleward flow. 

The northern high-latitude cyclone, which 
consists of part of the west wind drift, the 
Alaska Current and the Oyashio Current, is 
also quite clearly defined, as is the great cir- 
cumpolar flow of the Southern Ocean. In 
these data little evidence is found of any west- 
ward flow along the coast of Antarctica. 

In the equatorial region the west-flowing 
currents are seen to be separated by the east- 
flowing (North) Equatorial Countercurrent, 
and an east-flowing current south of the equator, 
equally well-defined in these data in the 
western ocean. If this slope downward to the 
south is geostrophically balanced, then there 
is an eastward flow in this region, approxi- 
mately symmetrical to the (North) Equatorial 
Countercurrent (REID, 1959). It is strong in 
the west and becomes weaker toward the east. 
The data in the eastern Pacific are not spaced 
in a fashion to answer very well whether it 
extends to the longitude of the Galapagos 
Islands, but there is a suggestion of a weak 
eastward flow almost that far. The various 
atlases of sea surface currents are of limited 
value in this region because the data are so 
sparse. The Current Charts of the Southwestern 
Pacific Ocean (U.S.N.H.O., 1944a) show an 
average of about one observation per month 
per one-degree square in the region between 
5° and 10° S and between New Guinea and 
165° W. The monthly pilot charts (U.S.N.H. 
O., 1950) show no good indication at any 
time of year of an eastward flow at 10° S in 
the eastern Pacific. 

The charts prepared by the Arr Ministry 
(1939) include all of the South Pacific Ocean 

y seasons. An eastward flow at 10° S is seen 
from November through April starting from 
the Solomon Islands. It reaches 140° W on 
the November through January chart, when 
some west winds are indicated at 10°S in 
the western Pacific (McDonatp, 1938), and 
160° W on the February through April chart. 
Although this chart includes the Southeastern 
Pacific there are areas covering more than 20 
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degrees of latitude and longitude which con- 
tain no data. The triangle with corners at 1° S 
and 55°S on the American coast and at 
50° S 158° W contains practically no data at 
all except in a narrow band along the coast of 
South America.- 

The distribution of eastward flow near the 
equator (Figure 1) is not entirely inconsistent 
with the conventional concept of the counter- 
current as a return flow in the doldrums. 
The data on Figure 2 reveal that when wind 
stress is averaged all across the ocean through- 
out the year the doldrums are centered very 
nearly at the equator. But if the area is broken 
into eastern and western parts at the 160° W 
meridian, then the average doldrums position 
is found about five degrees north in the 
eastern part and five degrees south in the 
western part. In the west the doldrums are 
even weaker in the southern part than in the 
northern part. 

The seasonal variation of the winds is such 
that in the eastern Pacific the doldrums are, 
on the average, north of the equator in all 
seasons. In the west, however, the average 
position of the doldrums is slightly north of the 
equator in the season June through August. 

The eastern doldrums do lie very near the 
equator in the season December through 
February when minimum wind stress com- 
puted by Hipaxa (1958) at five degree intervals 
is at 2.5°N. It is possible that this south 
equatorial countercurrent extends farthest east- 
ward at this period and in some years might 
carry warm central-ocean water across the 
cold Peru Current to the coast. It might thus 
be the source of the warm water phenomenon 
known there as El Nifio, which occurs occa- 
sionally in December or January. 

Over large areas, however, the limitations of 
the geostrophic approximation are obvious. 
That the contours alone do not describe the 
flow at the equator is to be expected; the 
arrowheads are placed there to emphasize the 
discontinuity. But even in higher latitudes 
only a few of the contours are closed curves; 
the others begin and end at the continents. 


B. The speed of the surface flow. 

The Figure ı is a graphical interpretation of 
the computed values of geopotential anomaly. 
A few of the points have been violated by one 
dynamic centimeter in the interest of smooth- 
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Fig. 2. a. The anomaly of geopotential distance (AD) between the zero- and the 
1,000-decibar surfaces along the eastern and western sides of the Pacific Ocean, 
in dynamic meters. The western side is taken here to be Kamchatka, the Kurile 
Islands, Japan, and Ryukyu and Philippine Islands, New Guinea and Australia. 
Roman numerals indicate the month in which each datum was taken. 

b. Difference in geopotential anomaly at the eastern and western boundaries of 
the Pacific Ocean, and the total zonal wind stress (dynes per cm of latitude) across 


the Pacific Ocean. 


ness. WOOSTER AND TaFT (1958) have estimated 
the errors in measurement of temperature, 
salinity, and depth and concluded that they 
result in a standard error of 0.004 dynamic 
meters in the geopotential anomaly at the sea 
surface with respect to 1,000-decibars. This is 
much less than the periodic and non-periodic 
variations reported by Rem (1956), in which 
values differing by as much as 0.08 dynamic 
meters were found at one position one week 
apart. The data on Figure I suggest that such 
wide fluctuations are rare, since the field is 
relatively free of local disturbances of this 
order over most of the area. Because there is 
some roughness which may not be real, the 
geostrophic speeds in the following discussion 
have been measured from the spacing of 
adjacent contours rather than from the differ- 
ence from station to station within the two 
contours. 

A comparison of the geostrophic speeds with 
the atlases (U.S.N.H.O., 1944b and 1947) 
has been made and the results are summarized 
as follows. 

1. In the Kuroshio the geostrophic speed of 
1.3 to 1.6 knots is slightly higher than most of 
Tellus XIII (1961), 4 


the atlas averages. The averages may have been 
damped by the shifting of the axis of the stream 
(MASUZAWA, 1957, and Fukuoka, 1958). 

2. In the slower west wind drift and Cali- 
fornia Current the agreement is to about 0.1 
knot, averaged over large areas. 

3. In the North Equatorial Current the 
agreement is to about 0.1 knot averaged over 
large areas except in the west where Figure 1 
shows a southward shift of the main stream, 
resulting in a reduction in speed at about 15° N 
of 0.4 knot. 

4. The northward flow past Formosa is 
computed as about 1.3 knots, while the atlas 
average is about 1.15. 

5. The westward flow south of the Aleutian 
Islands appears to be about 0.1 knot, from 
both methods. 

6. In the southern part of the Oyashio 
Current the atlas indicates about 0.2 knot to 
the east, but the geostrophic flow, as contoured, 
is about 1 knot. 

7. In the (North) Equatorial Countercurrent 
at latitudes below 5° and in the westward flow 
nearer the equator the geostrophic speed is 
consistently higher than the atlas speed. 
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5. The effect of wind stress 


One consequence of the neglect of the other 
terms (probably, in this case the stress associated 
with the vertical shear of velocity) is that the 
difterence in geopotential anomaly between 
the eastern and western sides of the ocean 
(Figure 2) must be interpreted as north-south 
flow at the surface. In Figure 2 a mean equator- 
ward flow is indicated between 37° N and 
42°S, and a mean poleward flow north of 
there. 

Along the east side of the ocean (Figure 2a) 
the geopotential anomaly is greatest near the 
equator. This is in accordance with what one 
might expect from the latitudinal distribution 
of solar radiation. Along the west side (here 
defined as Kamchatka, the Kurile Islands, 
Japan, the Ryukyu and Philippine Islands, 
New Guinea and Australia) the variation is 
much greater, with maxima near the tropics 
as well as at the equator. Although the varia- 
tion is also a consequence of the temperature 
distribution (since salinity is highest near 
the tropics and low near the equator and in 
high latitudes) some other factors must operate 
as well as local heating or the maxima would 
not be found so far from the equator. Among 
these one might consider the direct effect of 
the zonal wind stress (MONTGOMERY AND 
PALMEN, 1940). Values of wind stress have 
been computed by HıparA (1958) for each 
5° of latitude and longitude by seasons and 
by annual mean. From these the total zonal 
wind stress across the ocean has been calculated 
as a function of latitude and plotted over the 
latitudinal range of the data (Figure 2b). 
North of the equator the seasonal value for 
June through August is plotted also, since most 
of the oceanographic data north of the equator 
were from the period June through September. 
The oceanographic data south of the equator 
are from various months on either side of the 
ocean, and therefore no seasonal curve has 
been included. There is enough similarity in the 
curves to suggest that the wind stress is of 
great importance in accounting for the differ- 
ence in geopotential anomaly. 

If one integrates equation 2a zonally from 
the western to the eastern edge of the ocean 
and from some depth Z where the vertical 
shearing stress can be neglected, the velocity 
term will be negligible (since one may assume 
a negligible net meridional transport) and 
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fw: = PTS f (rodx (2) 


where (tx), is the surface wind stress. The 
term on the right can be evaluated from Hi- 
daka’s results and is plotted in Figure 2. The 
term on the left is the vertically integrated 
pressure difference across the ocean. In the 
data a value of Z can be found for each zone 
such that the pressure gradient and wind stress 
are balanced. If one assumes that the horizontal 
pressure difference is constant down to Z and 
negligible below, then the data in Figure 2 
suggest that an average value of Z would be 
about 100 meters. Both MONTGOMERY AND 
PALMÉN (1940) and AUSTIN (1958) have shown 
that the pressure differences extend to greater 
depths, however. 

More specifically, the velocity computed 
from the simplified equations does not satisfy 
the equation of continuity with a steady distri- 
bution of density unless the flow is east-west 
only. Since north-south flow also occurs in 
the ocean, the actual transport contains non- 
geostrophic flow, and transport calculated 
from the geostrophic equation will be in- 
complete. The same degree of incompleteness 
will be present in the geostrophic flow at the 
surface and is most immediately apparent in 
the appearance and disappearance of the con- 
tours at the coasts and in the indicated con- 
vergence at the equator. 

These discrepancies, however serious, are not 
unexpected, and their nature might have been 
anticipated, in part, from a qualitative con- 
sideration of some of the neglected terms. 
MONTGOMERY AND PALMÉN (1940, Figure 41), 
considering wind stress, the earth’s rotation, 
and the known surface currents, have prepared 
a schematic representation of the resulting shape 
of the sea surface near the equator that com- 
pares quite well with the corresponding part 
of Figure 1. : 

If one includes the vertical shear stress terms 
in equations ta and ıb, and multiplies the 
first by u and the second by v and substracts, 
one obtains 


uadp/dx + vadp|dy = uadr.|dz + vadty/dz. (3) 


If the local rate of change of geopotential 
difference is zero, the left side of (3) represents 
the time-derivative of the geopotential differ- 
ence or sea level along a streamline (in this 
case, with respect to the 1,000-decibar surface). 
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If 7, and ty are neglected, there is no change 
along a streamline. If they are included (and 
one notes that near the sea surface Or,/dz and 
Jty/dz will have the same signs, respectively, 
as (Tx)g and (ty), the surface wind stress), then 
it is seen that along a streamline sea level will 
rise where the products of wind stress and 
velocity are positive. That is, in downwind 
flow there is an uphill component and in 
upwind flow there is a downhill component. 

If one considers that the North Pacific is a 
series of narrow gyres, each extending zonally 
across the ocean, then the centers of the anti- 
cyclones would be expected to stand higher 
and the centers of the cyclones lower than the 
adjacent areas. Except in the western boundary 
currents where MUNK (1950), CHARNEY (1955) 
and MORGAN (1956) have shown that more 
complex processes dominate, the wind field 
that lies above these gyres will drive the water 
uphill or downhill as it reinforces or opposes 
the flow, and the contours will therefore 
undergo certain excursions with regard to the 
actual flow. 

In the northern subtropical anticyclone (Ca- 
lifornia and Kuroshio Currents and parts of 
the west wind drift and North Equatorial 
Current) the wind is everywhere in approxi- 
mately the same direction as the flow. In 
the greater part of this anticyclone (except in 
the western boundary current) the streamlines 
must be directed uphill and the contours must 
spiral outward from the center of the gyre. 
Those contours (1.5 to 1.8 dynamic meters) 
which spiral outward to the south far enough 
to cross the North Equatorial Current and 
reach the North Equatorial Countercurrent 
will have entered that part of a cyclonic gyre 
which flows upwind. Since the water will 
then be moving downhill, the contours must 
spiral outward again from the center of the 
cyclone, and some (1.6 and 1.9 dynamic meters) 
extend southward across the countercurrent to 
the Equatorial Current. They are now again 
in an anticyclone in which the westward flow 
is downwind and uphill. Therefore the contours 
spiral outward again, and some of them (1.4 
through 1.9 dynamic meters) intersect the 
equator, meeting their counterparts from the 
southern hemisphere. 

By a similar qualitative argument the dis- 
appearance of the contours (1.1 through 1.2) 
in the westward flowing limb of the subarctic 
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cyclone might have been anticipated. This ow 
is upwind, the height along a streamline must 
decrease, and if one considers a streamline 
near the coast in the Southeast Alaska Current 
this can occur only if the contours intersect the 
coast. 


6. Comparison with the Atlantic 


It is interesting to compare this chart with 
the corresponding one for the Atlantic Ocean 
prepared by Drrant (1941, Beilage XIV). He 
used METEOR data in the South Atlantic and 
a combination of data from several expeditions 
in the North Atlantic. 

One first notes that the Pacific values seem 
to be higher than the Atlantic values. It has 
been calculated (Reid, 1961) that the average 
value of the geopotential anomaly at the sea 
surface with respect to the 1,000-decibar sur- 
face in the Pacific is higher than in the Atlantic 
by 40 dynamic cm, and that with respect to the 
1,000-decibar surface the Pacific must stand 
40 cm higher than the Atlantic. 

A second feature of interest is the difference 
in complexity of the two charts. Part of the 
additional complexity of the North Atlantic 
may be the result of having to use data taken 
at quite different times. The South Pacific 
chart is also composed of data from many 
different years, yet does not have so many 
irregularities. The relatively simple pattern of 
large cyclones and anticyclones, roughly 
symmetrical about the equator, which charac- 
terizes the Pacific is apparent in the Atlantic 
also, but is confused by many small hills and 
valleys. 

The appearance and disappearance of con- 
tours at the coast that characterize the Pacific 
chart are found also in the Atlantic chart, to 
about the same extent. 
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AUESBER TNEDITERE 


The materials used in preparing the chart of geopotential 
anomaly at the sea surface with respect to the 1,000-decibar 
surface. 


From the available materials choices were made on 
the basis of the area and time of each expedition. The 
NORPAC, EQUAPAC, and EASTROPIC expeditions 
each covered large areas in relatively short (although 
different) periods. Stations from these three were chosen 
in preference to those stations of other extensive cruises 
(the Carnegie, for example) which overlapped them. 
From expeditions such as the Discovery II, Snellius, 
William Scoresby, and Derwent Hunter that covered 
the same area more than once, selections were made 
on the basis of season and need. In areas where data 
were scanty, overlapping was permitted. 

It is not possible to explain in detail each choice, partly 


because there are so many, and partly because even 
though some more suitable data became available during 
the work, the older choices were in some cases allowed to 
stand. 

The list contains those stations which were used and 
those stations which were considered and rejected without 
apparent fault except that the value seemed not to fit 
the field. For any cruise the stations not listed were 
either too shallow (a few of the listed ones have been 
extrapolated the last hundred meters), redundant, or 
obviously in error (showing instability or extremely 
unlikely values of temperature or salinity). 

T. S. Austin (in press) has reported a consistent depth 
error in the data from the HucH M. SmitH. The error 
has been corrected in the data used in this study. Part of 
the corrections were made through the facilities of the 
department of Oceanography, Texas A. and M. College. 
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The NORPAC Stations used (June—October 1955) 
es ee ee ee ei 


Number 


Stations Number | Station 
Vessel pe to : Reference 
considered SS used rejected 

PAN CX CAO} Eee Mrs I—162 143 142 No. 154 NORPAC COMMITTEE (1960) 
Oceanic Observations of the 

Tenyo Marne nee. ee 6—36 38 38 Pacific: 1955, The Norpac 

AS—5T DATA, pp. I—582. Univer- 

TMCS Stes Therese... I—39 82 82 sity of California Press and 

41—85 University of Tokyo Press, 

(Wanttaka Mar. 1:24 5—59 44 43 No. 47 Berkeley and Tokyo. 

Kagoshima Maru.......| 521—531 II II 

Ketten Maru, euere 565 I I 

USAR VOR Mara ee Dı—Dıı 7 Mi 

Brown BEAT eee I—69 43 43 

EIDME SWithe Sen cetera I—II4 76 76 

@Shorom Markus ER I—I9 29 29 

CE mes VA 
Saesumaree le 2—19 47 47 
21—49 
SOU EMTATU sare. sue alee 0e 448—496 33 31 No. 458, 
493 

Pasir Wiarih ec. 5 ose Ss IOOI—IO3I 22 22 

Shumpuy Mari’ =... 2... 1— 85 32 2 

Meryor Mar seen nee Cı—C5 16 16 

B3—Bi3 
303) 
30372305 
Total NORPAG. <6. 152. | 620 4 | 
The EASTROPIC Stations used (September—December 1955) 

ERENTESMIEH 10 er se I—24 20 20 kane, J. E., Austin, MS", 
and Dory, M. S. (1957) Pre- 
liminary Report on Expe- 
dition Eastrropic. U.S. Fish 
Wildiiie Serves pec moc 
Rep., Fish. No. 201: 1—155. 

TE IG YAS THe ace le Aare CREME 7—2: 76 76 Univ. CALIF., SCRIPPS 

26— 34 INST. OF OCEANOGRAPHY 
36—42 (1956) Data collected by 
A) Scripps Inst. Vessels on 
INN EASTROPIC Expedition, 
September— December 
1955. Unpub. Rep., SIO 
Ref. 56—28. 
Bike hove, OO One I— 24 70 70 No. 5 
3233 
3542 
4455 
58—74 
76—01 
HROtAÏMEASTROPIC LE oo 166 I 
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The EQUAPAC Stations used ” 
SS 
Stations ae Number Station 


considered used rejected 
1,000 m J 


Vessel Reference 


Jets Whe Sion ea ET I—79 79 79 g Austin, THoMAS S. (1957) 
Summary Oceanographic 
and Fishery Data, Marque- 
sas Islands Area, August — 
September, 1956 (EQUAPAC) 
U.S. Fish Wildlife Serv. 
Spec. Sci. Rep., Fish. No. 
217: 1—186. 


S Grate EI ne I—46 45 45 UNIV. CALIF., SCRIPPS INST. 
OF OCEANOGRAPHY (1957) 
Data collected by Scripps 
Inst. vessels on EQUAPAC 

Horizon ne seule I—47 43 43 Expedition, August 1956. 
Unpub. Rep., SIO Ref. 57 
—25. 


OTSOMEITIEEN EEE ENT 19 18 18 INSTITUT FRANGAIS D’OCEA- 
NIE, Orsom III, Croisiere 
Eguapac, September—Oc- 
tober 1956. Unpub. Rep. 


Satsuma tyes lern I—46 44 44 Supa, K., Eguapac Oceano- 
graphic and Meteorological 
Data (Preliminary Report), 
Satsuma, July—August 
1956. Japanese Hydro- 
graphic Office, Unpub. Rep. 


DTA AMATU en ee I—20 8 8 Tokyo UNIV.FISH., EQUA- 
PAC: Umitaka Maru, Egua- 
pac Expedition, October— 
November 1956. Unpub. 
Rep. 


SHOVORMArUE EE citer 5—15 10 10 JAPAN FISHERIES AGENCY, 
EQUAPAC Oceanographical 
and Meteorological Data, 
Shoyo Maru, August 20- 
November 3, 1956, Unpub. 
Rep. 


Kagoshima Maru...... 614 19 19 KAGOSHIMA UNIV. FACULTY 
616—618 OF FISH., Oceanographical 
620—634 Observations made during 
the International Coopera- 
tive Expedition EQUAPAC in 
July—August, 1956, by M. 
Wertens Marta. enter 601—622 19 19 S. Kagoshima Maru and by 
M. S. Keiten Maru, Unpub. 

Rep. 
HOLAlREOUAPACEE Re | 285 o 
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The stations used from expeditions made by single vessels 


Stations nt Number Station 
considered Don used rejected 


Vessel Reference 


Ghallénaen ere sees 220,270 10 10 OFFICERS OF THE CHALLENG- 
(8 Sept.—23 Oct. 1875)} 280—281 ER EXPEDITION (1884) Re- 
283—287 port on the Deep-Sea Tem- 
289 perature Observations of 
Ocean Water, Report on 
the scientific results of the 
voyage of H.M.S. Challeng- 
er during the years 1873— 
1876. Physics and Chemistry 
Vol. 1, Part III, pp. 1—2, 
pl. I—CCLVIII, Tables I— 

VII. 

also 
Wüst, G. (1929) Schichtung 
und Tiefenzirkulation des 
Pazifischen Ozeans: Institut 
für Meereskunde, Berlin, N. 
F., A. Geograph.-Natur- 
wissensch. Reihe, Heft 20, 

pp. I—64. 

IDÉES ee dar get eae ert LRU 2 3 REICHAR, A. C. (1911) Tem- 

(28 May—3 July 1908) 22 peratur- und Salzgehaltbe- 
stimmungen im südwest- 
lichen Stillen Ozean, 1910. 
Ann. der Hydrogr. und 
Mar. Meteor., XX XIX, Heft 
X, PP: 521— 527. 

Barnesie.err. IN 330.0% 35—97 62 61 No. 92 FLEMING, J. A., ENNIs, C. 
(26 Oct. 1928— C., SVERDRUP, H. U., SEA- 
28 April 1929) TON, S. L., and HENDRIX, 

W. C. (1945) Observations 
and results in physical 
oceanography, Scientific 
Results of Cruise VII of the 
Carnegie during 1928—1929 
under Command of Capt. 
Jin aeAulte Oceanogr, IB 
pp. I—315. 

Danae sso enges 3592— 43 43 CARLSBERG FOUNDATION 
(9 Nov. 1928—10 Aug. | 3655 (1937) Hydrographical ob- 
1929) 3666— servations made during the 

3788 Dana expedition 1928— 
1930. Dana-Report Vol. II, 
No. 12, pp. 1—46. 


SREINDSMER Scere sce ee Se 30 40 40 VAN RIEL, P. M., GROEN, P. 
(29 July—27 Nov. 33-35 and WEENINK, M. P. H. 
1929) 90, 95, 100 (1957) Oceanographic re- 

104—113 sults, quantitative data con- 

115—118 cerning the statics of the 

122, 128 East-Indonesian waters; 

Fa, 140 : depths of standard pressures 

147 and stability values. The 
Snellius Expedition in the 
eastern part of the East 
Indian Archipelago 1929— 
793508 Vol 1,2 earte7 pp: 
I—45. 
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tn nt D na re 


5 Number : 
SER Stations pre Number Station Röference 
aaa considered onan used rejected 


157 I 

(7 March— 5 Nov.1930)| 201, 205 
208 2210 
240, 262 
203 0272 
276, 301 
303, 319 
312 
324—327 
334, 349 
359, 353 
354a, 382 

William Scoresby....... 597, 606 27 27, DISCOVERY REPORTS (1949) 
(19 May—27 Aug. 609 Hydrographical observa- 
1931) 612—616 tions made by R.R.S. Wil- 

629, 638 liam Scoresby 1931—1938. 
646, 653 Vol. 25, pp. 143-280, 
668 Cambridge. 

671—688 

694, 701 

793 

705—7II 

719—722 

713472131. 

DISCOVERY elle en ee 873—926 150 150 DISCOVERY REPORTS (1941) 
(8 May—2 July 1932) Discovery : Investigations 
(9 Sept.—29 Oct. 1932) 956—994 Station List, 1931—1933, 

Vol. 21, pp. I—226, Cam- 
bridge. 

(28 Jan.—7 Feb. 1936) | 1662—1679 (1944) Station List, 1935— 

1937, Vol. 24, pp. I—196. 
(12 Jan.—g March 1938)| 2174—2220 (1947) Station List, 1937— 
2220—2280 1939, Vol. 24, pp. 198—422. 
(20 Oct. 1950— 22 June 2728 (1957) Station List, 1950 — 
1951) 2734—2741 1951 Vol. 28, pp. 300—398. 
2768— 2771 
2780—2782 
2789—2798 
2817—2821 
2831—2841 
USSIGANRELEE cose 26, 28 4 4 Barnes, C. A., and THoMp- 
(10—11 Aug. 1933) 2, 34 SON, T. H. (1938) Physical 
and chemical investigations 
in Bering Sea and portions 
of the North Pacific Ocean. 

USEGEICHELANERTERER 102, 105 6 6 Univ. Wash. Pub. in Ocea- 
(18—24 Aug. 1934) 106 nogr., Voll"3;4Noyseepp: 

I17—119 35—79 and Appendix, pp. 
I—164. 

NLDATLOSS ACER IR 43—I01 34 33 No. 93 BRUNEAU, L., JERLOV, N.G. 
(29 Aug. 1947— 103— 112 and Koczy, F. F. (1953) 
12 Jan. 1948) 114—125 Physical and chemical 

127 methods, Appendix Table 2, 
130—132 Rep. Swed. Deep-Sea Ex- 
134— 136 ped., Vol. III, Physics and 
138—142 Chemistry, No. 4, pp. XLII 
144— 150 —LV. 

377 2 
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Vessel 


Stations 


considered 


Number 


to 


1,000 m 


Number 
used 


Reference 


(Cruise 5) 
(30 June—6 Aug. 1950) 


USS EPCE(R) 857 
(Shuttle) 


@SDoro, MAR... 
(31 May—20 June 
1954) 


Derwent Hunter..... 
(23 Nov.—13 Dec. 
1954) 


(29 Aug.—13 Dec. 
1955) 


(13 April—5 May 1956) 


Orsom, WEE 2.2.2.0 aoe 
(Astrolabe) 


(5—17 May 1958) 


Do mule ment 6 25 


(565) 
(13 Oct.—14.Nov. 


1956) 
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I—12 
50—51 


39554 


II—28 


160—195 


96—194 


79—98 


I—I4 


I—23 


16 


29 


12 


22 


14 


2I 


29 


22 


520 


CROMWELL, T., and AUSTIN, 
T. S. (1954) Mid-Pacific 
Oceanography, Parts II and 
III, Transequatorial Wa- 
ters, 1950—1951. U.S. Fish 
Wildlife Serv. Spec. Sci. 
Rep.,, Fish. No. 137. 


U.S. NAVAL ELECTRONICS 
LABORATORY, Shuttle Ex- 
pedition, 24 April to 4 June 
1952. Unpublished mate- 
rial. 


MısHIMA, S., and NISHI- 
ZAWA,S. (1955) Reporton the 
Hydrogr. Invest. in Aleuti- 
an Waters and the South- 
ern Bering Sea in the Early 
Summers of 1953 and 1954. 
Bull. Fac. Fish., Hokkaido 
Univ V.01.26, N0.22,5pP. 
85—124. 


COMMONWEALTH SCIENTIFIC 
INDUSTRIAL RESEARCH 
ORGANIZATION, AUSTRALIA, 
DIV. FISH (1956) Oceanogr. 
Station List, Onshore Hy- 
drological Invest. in Eastern 
and Southwest Australia, 
1954. Vol. 24, F.R.V. Der- 
went Hunter, pp. 54—100. 
(1957) Vol. 27, 1955, PP- 
66—145. 

IGY REPORTS OF THE COM- 
PLEX ANTARCTIC EXPEDI- 
TION OF THE ACADEMY OF 
SCIENCES OF THE USSR 
(1958) Hydrological, Hy- 
drochemical, Geological and 
Biological Studies, Re- 
search Ship Ob 1955— 19506. 
Hydro-Meteorological Pub- 
lishing House, Leningrad, 
pp. I—214 

Rotscu1, H. (1958) Orsom 
III, Océanographie Phy- 
sique, Croisiére Astrolabe. 
O.R.S.T.O.M., I.F.O. Rapp. 
Sie, INS, 1806 AO) 
Nouméa. 

RorscHı, H. (1958) Or- 
som III, Océanographie 
Physique, Croisière 56—5. 
OR SOMME O7 Rapp: 
Sci Non, pp 1-34, 
Noumea. 
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Stations er Number Station 


considered used rejected 
1,000 m J 


520 2 
NAtyazı2 er ae sities 3652— 3683 21 21 : 136 stations in the W. Paci- 
i fic, Vityaz Cr. 25, 28 June— 
11 Oct. 1957. Unpublished 


Vessel Reference 


report made available 
through World IGY Data 
Center A. 


"102 stations in the W. Cen- 
tral Pacific, Vityaz Cr. 26, 
5 Nov. 1957—27 Feb. 1958. 
Unpublished report made 
available through World 
IGY Data Center A. 


Wikies Coles ete 3801—3875 57 55 No. 3855, 
3856 


UNIV. CALIF., SCRIPPS INST. 
OF OCEANOGRAPHY (1958) 
Physical and chemical data. 
Chinook Expedition, 1956, 
Mukluk Expedition, 1957, 
and Downwind Expedition, 
1957—1958, Unpub. Rep., 
SIO Ref. 58—85. 


Cruise 3 of the Ob in the 
Indian, Pacific,and Atlantic 
Oceans, 15 Jan.—18 June 
1958. Unpublished report 
made available through 
World IGY Data Center A. 


COMITE LOCAL D’OCEANO- 
GRAPHIE ET D’ETUDE DES 
COTES DE NOUVELLE-CALE- 
DONIE (1958) Resultats des 
observations scientifiques 
du Tiare, Croisière Bounty, 
20—29 Juin 1958, Rapp. 
Sci.I.F.O., No. 7, pp. I—20, 
Nouméa. 


FONZONE Erler Io—II 31 29 No. 42 
(Downwind) 13—44 43 
(5 Nov. 1957—23 Feb. 
1958) 


OBICTUISSS 3 ern 334—351 43 43 
423—450 


(Bounty) 4—14 


Umitaka#Marukens «ae I—27 13 13 


Coral Sea Cruise, 7—14 
(Coral Sea) 


Jan. 1959. Data kindly 
made available through 
Commonwealth Scientific 
Industrial Research Orga- 
nization Marine Labora- 
tory, Australia, and Tokyo 
Univ. of Fisheries, Japan, 
Unpublished Material. 


ASSOCIATION OF AGRICUL- 
TURAL TECHNOLOGY (1954) 
Northern area oceanogra- 
phic data 1887—1953. Re- 
port of the invest. of long 
range cold weather fore- 
casting. No. I, pp. I—556. 
= Tokyo. 

Total, single vessel | | rg 


EXPECIMONS ca emronene 


Grand total of all stations 1,792 | II | 
_—_—— 0 S 
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The Stability of Thermoclinic Jets 
By MELVIN E. STERN, Woods Hole Oceanographic Institution 


(Manuscript received February 6, revised version August 28, 1961) 


Abstract 


A laminar jet, in a “shallow layer’’ of water of density @ — Ag is in geostrophic equilibrium 
above a much deeper and resting layer of water of density 0. We compare the quasi-hydrostatic 
energy releasing process in this thermoclinic model with other baroclinic and barotropic models. 
It is first shown that a small amplitude disturbance can grow only if there is a transfer of kinetic 
energy from the horizontal variations in the mean jet and that conversions of mean potential 
energy can only occur concomitantly with, or as a finite amplitude result of this process. We 
also examine the perturbation equation for values of the stream Rossby number which are much 
less than unity and show that the finite depth of the thermocline inhibits inflectional instability 
of the jet. The generalization of Lord Rayleigh’s theorem shows that the jet is stable if the 
gradient of potential vorticity does not vanish. In general, there is a critical value of the ratio of 


the radius of deformation (fo VeHAolo) to the half width of the jet above which quasi- 
geostrophic disturbances will start to grow. It is proposed that meanders of the synoptic Gulf 
Stream will develop when the jet becomes ‘‘excessively”” deep and that the effect of such insta- 


bilities is to limit the depth of the main thermocline in the interior of the ocean. 


1. Statement of the Problem 


We shall call the geostrophic laminar cur- 
rent, which is shown schematically in Fig. 1, a 
thermoclinic jet since it is a useful abstraction 
(RossBy, 1936); STOMMEL, 1953) of the ther- 
mal structure of the Gulf Stream. It consists 
of a homogeneous layer of light water of 
mean depth H, lying above a much deeper 
layer of water which is denser by an amount 
40<o. The latter is resting in a coordinate 
system which is rotating at w=f/2 rads. /sec. 
The velocity in the upper layer is V*(x) 
directed along the y axis (or into the page) and 
this is balanced by the hydrostatic pressure 
gradient due to the elevation of the free 
surface. The characteristic width of the jet is 
denoted by 2L and V is the maximum speed. 

The mutual adjustment of the free surface 
and the interface, which is necessary to bring 
about this inviscid equilibrium, allows us to 
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express the hydrostatic pressure gradient in 
the upper layer in terms of its thickness [h* (x)] 
as: gdh*/dx where g’=gA@/o is a “reduced” 
value of gravity. Under the conditions stated 
below! a similar simplification may be made for 
a perturbation [h’(x, y, t)] which is introduced 
upon this mean field, so that if hr=h* +h’ 
denotes the total thickness of the upper layer 
then g’vhr is the horizontal pressure gradient 


1 The hydrostatic (or shallow water) approximation 
may be defended when the mean depth of the top layer 
(H,) and the bottom layer (Hp) are both small compared 
to the horizontal scale of the perturbation (L). On the 
other hand it may be shown that the ratio of the pressure 
gradients in the two layers is small to order H/H4 
provided the stream Rossby number (R) and the rota- 
tional Froude number (F) are of magnitude unity or 
less (see eq. 10 for definitions of R, F). Thus the validity 
of this approximation and the analogy of the two layer 
model in Fig. 1 with the one layer model in Fig. 2, 
depend on the inequalities H,/L<Hp/L<ı when R <1 
and F<ı. 
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Fig. 1. A schematic diagram of a two layer thermoclinic 
jet. The lower layer of density o is infinitely deep and 
resting. The geostrophically balanced velocity field V * (x) 
is directed along the y axis (into the page) and its maxi- 
mum value is V. The section shown is supposed to occur 
at very large distances from the axis of the rotating system 


(f-plane). 


therein. When these conditions are fulfilledthen 
the dynamical model in Fig. 1 is analogous to the 
single fluid model shown in Fig. 2, in which 
gravity is reduced to g’. Although short-wave 
Helmholtz instabilities, associated with the 
interfacial shear in Fig. 1, have been filtered by 
the hydrostatic approximation it will be seen 
that either model has available potential energy 
as well as kinetic energy of horizontal shear. 
Both of these energy sources can be tapped by a 
superimposed perturbation. In Section 2 we 
shall examine certain integrals of the motion 
and show that the kinetic energy releasing 
process is the more fundamental one in the 
thermoclinic model. If an initially small disturb- 
ance releases potential energy from the mean 
flow it can only do so with a concomitant 
release of kinetic energy by the action of the 
Reynolds stress. This model then is funda- 
mentally different from, and a logical extreme 
to, models having a uniformly distributed 
baroclinic zone (CHARNEY, 1947; EADY, 1949) 
and which release only potential energy in 
the limit of small Rossby number. Kortchm’s 
(1932) model of a Margules front separating 
two layers of the same thickness exhibits a 
similar potential energy conversion for long 
wave-length disturbances. SOLBERG (1930) has 
also considered the polar front instability in 
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which both fluids are of infinite vertical extent. 
My understanding is that the instability mech- 
anism here is due to vertical shear across the 
interface. 

The classical Rayleigh-Tollmien theory (Lin, 
1955) of two-dimensional inviscid shearing How 
has established that a symmetric jet is unstable. 
Haurwıtz and PANOFSKY (1950) have suggest- 
ed this as a mechanism for meanders in the 
Gulf Stream ‘(see FUGLISTER and WORTHING- 
TON, 1951, for the observational data which 
supports this picture). For zonal barotropic 
currents on a rotating sphere Kuo (1948) has 
generalized Lord Rayleigh’s criterion to show 
that an extreme value of absolute vorticity is 
necessary for instability of the profile. The 
thermoclinic model allows a modification of 
this energy releasing process by virtue of the 
vertical velocities which work against the 
gravitational field. Nevertheless, a simple gen- 
eralization of the foregoing criteria is possible 
which shows that “large-scale” instabilities 
will occur only when the profile has a maxi- 
mum (or minimum) of potential vorticity. We 
point out the stabilizing influence of the finite 
depth of the thermocline and suggest the 
relevance of this quantity to the problem of 
the stability of the Gulf Stream. 


2. The Energetics of Thermoclinic Jets 


If V(x, y, t) denotes the total horizontal com- 
ponent of velocity in Figure 2 (or in the 
upper layer of Fig. 1) and hr(x, y, t) the 
thickness of the layer then the momentum 
and continuity equations are— 


N 4 V-0V +20xV= = g'vhr (1) 
Ohr 
> sain -Vhr=o (2) 


Le 


FREE SURFACE 


® h*(x) 


ke 


RIGID BOTTOM 
Fig. 2. A single layer system which is equivalent to Fig. ı 
with respect to the dynamics of quasi-hydrostatic 
disturbances. g’=gAo/o is the “reduced” value of 
gravity. 
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Let us resolve the velocity into a mean com- 
ponent V*(x) directed along the y axis and a 
perturbation V’(x, y, ft) which has magnitudes 
(u’, v’) along (x, y). With a similar decomposi- 
tion of hp=h* +h’ the linearized perturbation 
equations are— 


2 


ON 
> +2® xV’+V*.vV’ + V.vVt= -g’vh 
(3) 
oh 


ee (V*h’ + V’h*) = o (4) 


An interesting energy equation may now be 
formed by multiplying (3) by h*V’ and (4) 
gh’. Denoting an average value over the 
entire horizontal plane by a bar we then get 

AI 


4 I [4 
REN ce ENE UNG 4 
dt 2 2 


+h*V’.V’-vV*= -gh*V’.vh (5) 


: = gh?+ghv-V*h’+g’hv-V’h*=o (6) 
In proceeding further we make use of the 
boundary condition that the normal compo- 
nent of V’ vanishes at the x-boundaries and 
the downstream variation is well behaved 
(e.g. periodic in y). Note that the second term 
in (5) vanishes because the undisturbed velocity 
field has no horizontal divergence and also 
v -h*V* =o. Likewise, the second term in (6) 
vanishes and these equations become— 


u h*¥V2= - h*V -V-vV*+ghv-V'h* 
(7) 

23 yin - go VF 
33 £h re ees v-Vh (8) 


Equation (7) may be interpreted as a statement 
of the rate of increase of the perturbation 
kinetic energy in the undisturbed depth h*. 
The first term on the right represents a con- 
version of kinetic energy as the result of the 
Reynolds stress acting on the horizontal shear 
of the mean field. The second term represents 
work done against the gravitational field and 
will have a counterpart in the increase or 
decrease of the potential energy of the mean 
field. Note that in (8) there is an equal and 
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Opposite transfer to perturbation potential 
energy so that the sum of (7) and (8) is— 


a(t ——— 1 
Sean. roles TJ". 
xa V ré] h*V’.V’.vV* 


= — h*yv Boe (9) 


Since both integrands on the left side of (9) are 
positive definite it follows that spontaneous 
instabilities can only occur if - h*V’- V’- vV* 
is positive, i.e., the mean field must do work on 
the perturbation through the Reynolds stress. 
If this term were zero or negative neither V’ 
nor h’ could increase with time. This energy 
equation implies that the simple meander 
theory of STOMMEL (1953) is untenable, for 
in his model h*(x) is linear and JV*/0x =o. 
In that case 1/,h*V’? +1/,g’h’? is conserved 
and a small disturbance could not grow unless 
the boundaries were squeezed. We also con- 
clude that in the thermoclinic model, which 
we are about to consider in more detail, the 
shearing mechanism is fundamental to insta- 
bility despite the availability of potential energy 
as represented by the departures of the surfaces 
from equi-potentials. This statement does 
not imply that there is no conversion of mean 
potential energy to the perturbation, but 
rather that if such occurs it can only be as a 
concomitant to the shearing mechanism. It 
should be borne in mind that this conclusion is 
based on the initial hypothesis that the bottom 
layer is very deep. Accordingly, this approxi- 
mation may have the effect of eliminating cer- 
tain modes of instability which are of interest. 


3. The Stability of Quasi-Geostrophic 
Disturbances 


The stability of the dynamical system (eq. 
(3) and (4)) depends on two “external” non- 
dimensional parameters: 


R=V/fL (10) 


The stream Rossby number (R) measures the 
ratio of inertial to coriolis force, for pertur- 
bations which have characteristic horizontal 
dimensions of order L (the half width of the 
jet) or greater. The rotational Froude number 
(F) is a measure of the vertical velocity of the 
perturbation, since the ratio of total horizontal 
convergence to either component (du'/Ax) is 


F=fPL2/¢ Hy 
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of order RF. One would like to discuss the 
stability characteristics for all values of R and F, 
using suitable profiles which describe the form 
of the mean flow. However, the analytical 
difficulties are such that we must satisfy our- 
selves with partial solutions, each of which is 
valid over certain ranges of the parameters. 
For example, when F < r the barotropic non- 
divergent theory is applicable and the necessary 
condition for instability requires that the jet 
have an inflection point. The case which is 
considered below is equivalent to an asymp- 
totic? expansion of (3) and (4) in the Rossby 
number and is valid for R<ı and F of order 
unity or less (as compared with the conclusion 
in Sec. 2 where the Rossby number may be as 
large as unity). 

Equations (1) and (2) imply the conservation 
of the vertical component of potential vortic- 
ity. If &*(x) denotes the undisturbed value of 
relative vorticity and £’(x, y, t) the perturbation 
vorticity then we have: 


0 Demo TC 


++ 
(rer) ° = 


h*+h 
(11) 


When linearized eq. (11) becomes 
AR TE A eas 
FG) x (x) AL Be 


RN 
(+7 IL tu LE Lo (32) 


At this stage we introduce the geostrophic 
approximations : 


oh’ , 4 oh’ 
ete 
olf dy 


and therefore “correct to order R” eq. (12) 
becomes 


? The y-wavelength is restricted to be of order L or 
greater and the time scale (period) of order L/V or 
greater. These assertions lead to a self consistent expansion 
of (3) and (4) whose zero order term is identical to that 
(15a) which is obtained by the simple derivation given 
below. 
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Let 


h' = o(x)e#@@-%, c= a +ib (14) 


and denote P(x) =(f+ &*)/h* as the undisturbed 
potential vorticity. We then obtain the eigen- 
value equation 


8°) = [at +fPO)IEIB- page =o 


p(+ ©) =0 (15) 


where a prime denotes differentiation with 
respect to x. We now parallel Lord Rayleigh’s 
technique (LIN, 1955) by multiplying (15) by 
the conjugate eigen-function (6) and average 
the result over all values of x. When the im- 
aginary part is set equal to zero one obtains 
+ © = 

ip MTP (x) pe at 

b| (Wr 4 


(16) 


It then follows that if there are unstable eigen 
functions (i.e. b >o) then P’(x) must vanish at 
some point which is at a finite distance from 
the centre of the jet. If P’(x) were always of 
the same sign then (16) could not be satisfied. 
It is therefore necessary that the potential 


vorticity 
P(x) =A RS 


(17) 


have a maximum (or minimum) interior to 
the endpoints of the jet, in order that the latter 
be unstable. The derivation of this criterion is 
such that it can only be relied upon in the 
limit of small Rossby number (£* <f). Utilizing 
this approximation we may alternatively con- 
clude that the flow is stable if 


~ rer] V* = fi/g'h*) 


is never zero. Let us consider the class of veloc- 
ity profiles for which V*"/V’* is always finite 
and let max. V*"/V* denote its greatest value 
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(for example: V*(x) #0 except at x=+0 
where V* goes to zero as some positive power 
of 1/x). From the above equation it then 
follows that the mean field is stable if 


fi ver) 
En > max. Ty (x) 
> ee h* (co) 


(18) 


Therefore a given V *(x) profile (corresponding 
to sufficiently small Rossby number) can be 
stabilized with respect to quasi-geostrophic 
modes by a uniform reduction of the thickness 
of the upper layer (other things being equal) 
to the point where (18) is satisfied. 

On the other hand, it is easy to establish 
instability at sufficiently large h* by noting 
that in this limit (F <1) P(x)>o while h* P’(x) 
— d£*/dx. In this case (15) reduces to 


#6) a8 - TO) =o 


b(s =) =0 (19) 


which is mathematically identical to the 
classical equation for two-dimensional inviscid 
instability. Tollmien’s theorem (LIN, 1955) 
tells us that (for symmetrical V*) there are 
marginally unstable solutions in the vicinity of 
certain discrete wave-numbers («). Therefore, 
with a given profile, there must exist a maxi- 
mum critical depth [(H:%)]|c. above which 
the mean field is unstable and this number is 


bounded by: 


Fe 
g(Hıo)e. 


The phase speed of the neutral disturbance 
is equal to the mean speed of the jet at the 
potential vorticity extremum. 

More precise specifications of these quanti- 
ties require the solution of the eigenvalue 
equation (15), which is equivalent (correct to 
order R) to the following non-dimensional 
equation: 


OFS < max. 


W'(&)-F 


ve Fin 
d(=+%)=o (15a) 
where I 
f= xT, V*(x) = VW(E) 
(> clV % = aL 
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Although (15a) is difficult to solve in general, it 
is possible to obtain particular solutions for the 
case of the “Bickley Jet”: W=sech?x & (Savic, 
1941). By direct substitution one finds that 
&=sech? € satisfies (15a) provided 1/, «2 = 
= Co = 1/6 (4 — F). Therefore long quasi-station- 
ary waves (%~0, € — 0) are neutral at F~ 4. 
Using Tollmien’s perturbation technique (Lin, 
loc. cit.) one can readily show that unstable 
solutions exist in the neighborhood of the 
neutral ones and that F<4 is a sufficient 
condition for the instability of a Bickley Jet. 


4. Discussion 


Let us use the Bickley Jet profile to make 
some order of magnitude estimates of the 
depth of the upper layer which is necessary for 
instability. Consider water with a vertical 
temperature difference of 10°C (or air with 
a 1° C increase in potential temperature), f= 
=10-* sec}, and a jet half width of L= so 
km. Then F=4 when H,.=250 meters. 
When this critical depth is exceeded the jet 
will be unstable to long waves provided the 
Rossby number is sufficiently small. To ful- 
fill this, and validate the asymptotic theory, 
the amplitude of the jet must be less than 
V<fL=s meters/sec. 

As regards the Gulf Stream, it should be 
pointed out that observations of the mean 
shear indicate that the characteristic Rossby 
number is not small but of order unity. 
Furthermore the restriction on the depth of 
the bottom layer in this model is such as to 
preclude other types of modes which are 
associated with the baroclinic structure of the 
jet. Nevertheless it is believed that the criterion 
for instability which is, expressed by an ex- 
tremum in the potential vorticity has a range 
of applicability which is greater than that which 
is indicated by the formal limits of validity of 
this theory. The relevance of this criterion to 
the stability of the Gulf Stream is supported 
by the observation of STOMMEL (1958 Fig. 65, 
based on a section at Lat. 38° N) that the layer 
between 17°C and 19°C is approximately 
one of constant potential vorticity. The 
disruption of the Gulf Stream north of Cape 
Hatteras may act as a control on the maximum 
depth of the thermocline in the interior of the 
ocean as well as along the western boundary. 
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Size Distribution of Large Droplets in Precipitating Clouds 


By T. OKITA, Geophysical Institute, Asahikawa Branch, Hokkaido Gakugei University, Japan. 


(Manuscript received October 28, 1960, revised version July 10, 1961) 


Abstract 


Measurements have been made of the size distribution of large droplets in various clouds at a 
station situated on Mount Kuro. The results show that most nimbostratus, cumulus congestus 
and cumulonimbus clouds contain large droplets of radii above so u. The spectral distribu- 
tions fit fairly well the relation N= Nor~* where « has values between 4 and 11. The concen- 
tration of large droplets varies widely and values above 10° per m? are frequently observed in 


precipitating clouds. 


1. Introduction 


Measurements of the size and number of 
cloud droplets have been made in a number of 
investigations, however, most of these have 
been concerned with non-precipitating clouds 
and there is still a lack of data with regard to 
precipitating clouds. The investigation describ- 
ed here reports on a number of measurements 
of the size distribution of large cloud droplets 
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observed on Mount Kuro (altitude 1,900 m, 
sce Figs. 1 and 2) in Hokkaido, Japan, during 
the summers of 1955—58. 


2. Gravity settling method 


In order to avoid the difficulties encountered 
when using impactor or oil or magnesium 
oxide coated slides (exposed vertically) for 
measuring sparsely distributed large droplets, 
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the following technique was adopted. Magne- 
sium oxide coated glass slides or water-blue 
dye coated celluloid film strips (OKITA, 1958 
a), having an area of 2 x 4 cm?, were placed 
in a rectangular hole in the center of a thin, 
streamlined wooden plate of the size 30 x 30 
cm’, Fig. 3. During the measurements the 


Fog Slide or film 
—h ee 
| ! 
R——30 Cm — + 


Fig. 3. Horizontal droplet sampler. 


plate was kept in a horizontal position in the 
cloud air. When large raindrops fell, the plate 
was covered with a woolen cloth to prevent 
splashes from influencing the measurements. 
In these cases the slides or the strips were 
placed on top of the cloth. 

The counting and measuring of the stains 
was made with the aid of a microscope, having 
a magnification of 1: 20 or I: 50. 

The size distribution of raindrops with 
radii larger than 0.2 mm was also measured, 
using sooted paper. The art-papers, sooted 
with smoke from kerosene, were horizontally 
exposed to the rain. Upon contact with the 
paper, the raindrops produced circular stains, 
whose size was a function of the diameter of 


~ 
S 


Stain Diameter mm 
G 


/ 2 
Raindrop Diameter mm 


Fig. 4. Relationship of raindrop diameter versus stain 


diameter on sooted paper. 


TO Rei: 


the original drop: A calibration curve is 
shown in Fig. 4. 

In evaluating the results it has been assumed 
that the raindrops and the cloud droplets 
will fall on the slides with terminal velocity 
v. If the number of drops or droplets hitting 
1 cm? of the slide or the film is n,, the corre- 
sponding number of droplets per cm? of air is 


Ue ns|v (1) 


Naturally, turbulent motion in the air will 
to some extent influence the results of these 
calculations. The influence may not be too 
critical, however, as a laminar boundary layer 
of a thickness of several millimeters will be 
formed along the plate. In this layer the 
vertical velocities will be considerably less 
than the terminal fall velocity of the droplets. 
According to wind tunnel measurements the 
boundary layer is always laminar at wind 
speeds below 4 m/sec. At velocities above 8 
m/sec. transition to a turbulent boundary 
layer will occur, but as the sampling was 
usually made at wind speeds below 5 m/sec. 
this case may be excluded (GOLDSTEIN, 1950). 


3. Liquid water content check 


In order to check the reliability of the sam- 
pling techniquea comparison was made between 


Rain shelter 


INTIME] ITTTTTTT 
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Fog 


Q/2 mm 
Cupper Wire 
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Fig. 5. Cylindrical wire screen. 
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the cloud liquid water contents determined by 
this method and those obtained using a wire 
screen. A cylindrical wire screen shown in 
Fig. 5 was used to collect the cloud or fog 
water (Tabata, Fujioka and MATSUMURA, 
1953). The amount of cloud water AM 
collected in time At is represented by the 
equation: 


AM = EW, VAAt (2) 


where W,, is the cloud liquid water content, 
A the cross-section area of the wire screen 
(174 cm?) and V the wind velocity. According 
to MATSUMURA (1953) the collection efficiency 
E is given by E= 0.11 V°-® where Vis measured 
in m/sec. 

If n; is the concentration of cloud droplets of 
radius r; per cm? of the air derived from the 
equation (1) the total liquid water content 
W, is calculated from the equation 


4 
34. 
“700 Zr; n; 
1 


W.= 


where o is the density of water. 

A simultaneous measurement of W,, and 
W. was made from 1605 to 1615 JST, August 
19, 1958, at a station on Mt. Kuro. The 
water-blue coated celluloid films were horizon- 
tally exposed in the cloud for ten seconds 
every minute. The amount of cloud water 
collected by the wire screen and the wind 
velocity were also measured for each minute. 
The calculated mean value of W. was 0.39 
g/m3, while the wire screen method gave 
the mean value of W, of 0.45 g/m?. Consid- 
ering the low accuracy of cloud liquid water 
content measurements, the results of these 
two methods show quite a good agreement. 
The mean wind velocity was 2.3 m/sec. The 
typical cloud droplet size distributions during 
the measurements are shown in Fig. 15. 
Droplets with radii 5 to 20 microns yielded 
most of the cloud water. 

Measurements were also made in sea fog 
from 1300 to 1400 JST on July 26, 1952, at 
Ochiishi on the eastern coast of Hokkaido 
(Fig. 1, Oxrra, 1953). The results shown in 
Table ı also indicate close agreement between 
the fog water contents as determined by 
these two methods. The wind velocity during 
the measurements was about 5 m/sec. The fog 
contained droplets of radii ranging from 2 to 
Tellus XIII (1961), 4 


Table 1. Fog liquid water content g/m? 


No. I | No. 2 


Wire screen method Wy, 
Slide method W, 


0.096 


0.097 
0.109 


0.093 


60 microns, but the droplets of radii 5 to 25 
microns yielded most of the water. 

The above measurements confirm the 
validity of the new sampling technique for 
droplets of radii between $ and 25 microns 
when the wind velocity is 5 m/sec. or below. 
In the case of large droplets of radius about 
so microns, no direct confirmation was ob- 
tained. It may, however, be supposed that 
the influence of vertical air motion, even if 
existing, would be completely negligible for 
these larger droplets. 

In order to determine the limitations of the 
method for small droplets, a comparison was 
also made with results obtained when using 
the impactor method. An impactor with a 
hole diameter of ı mm was used. Fig. 6 shows 


2 Gravity Settling Method 
©] Impactor Method 
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Fig. 6. Droplet size distribution of radiation fog at 
2250 JST, Dec. 10, 1959 in Asahikawa. 


O51 


$12 


one example of measured fog droplet size 
distributions. Both techniques give the same 
concentration for droplets of radii between 3 
and $ microns. For droplets with radii larger 
than 5 microns the gravity settling gave 
larger concentration than the impactor indi- 
cating that the collection efficiency of the 
impactor was reduced for large droplets. The 
gravity settling method gave inaccurate results 
for droplets of radii less than 3 microns, ob- 
viously due to the fact that the fall velocities of 
these droplets are too small to give accurate 
determinations. 


4. Measurements in July and August of 1956 


During the periods July 2s—29 and August 
II—I4, 1956, an extensive study of cloud 
droplet size distribution was made on Mt. 
Kuro. Photographs ot the cloud covering 
Mt. Kuro were also taken at Kamikawa, 
Kiyokawa or Rikuman (Fig. 2) in order to 
determine the form and the thickness of the 
clouds. 


(a) July 25 

In the early morning (0700 JST) an alto- 
stratus cloud with base as high as 2,600 m 
covered the whole sky but there was no cloud 
around Mt. Kuro. At about 0810 JST a cumulus 
cloud formed around Mt. Kuro and soon 
covered the mountain. At about 0820 JST 
the cloud tops intruded into the base of the 
altostratus. On Mt. Kuro drizzle began to fall 
at 0835 JST. At 0905 JST the altostratus cloud 
began to break and the tops of the cumulus 
cloud became visible at 0935 JST. The thick- 
ness of the cumulus cloud was about 950 m. 
From 1000 to 1005 JST the drizzle intensified 
again. After 1030 JST the drizzle ceased but 
Mt. Kuro was still covered with cumulus 
cloud until evening. The vertical development 
of the cloud was suppressed by a subsidence 
inversion and the thickness of the cloud varied 
between 680 and 1,350 m. Drop measurements 
during this time showed that the most frequent 
droplet radius was below 25 microns. Three 
examples of the cumulus cloud droplet size 
distributions (0845, 1005 and 1320 JST) are 
shown in Fig. 7. 

The observations suggest that the drizzle 
around 0845 JST was induced by the intro- 
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Fig. 7. Droplet size distributions on July 25 and 26, 1956. 
Weather at Mt. Kuro: 


Jul. 25; Air temperature at 0800 JST was 9.5° C. 
Wind direction was ENE to ESE in the morning 
and N to NW in the afternoon. 
Wind velocity was below 2 m/sec. 
Jul. 26; Air temperature at 0900 JST was 12.5° C. 
Wind direction was SW to NW. 
Wind velocity was 0.3 to 2.8 m/sec. 


duction of large cloud droplets at the top of 
the cumulus cloud. It is uncertain, however, if 
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such a mechanism was also effective for the 
drizzle occurring around 100$ IST. 


(b) July 26 

From early morning a light rain fell from 
nimbostratus clouds, the amount of rain fall 
being o.1—0.3 mm/hr. The cloud base was 
generally several hundred meters above Mt. 
Kuro and the mountain was hardly covered 
with cloud. One example of the cloud- and 
rain-drop size distributions (0920 JST) is 
shown in Fig. 7. 


(c) July 27 

Mt. Kuro was covered with nimbostratus 
cloud throughout the day. Drizzle with 
maximum drop radii of 0.3 to 0.4 mm was 
falling continuously now and then accom- 
panied by light rain. Three examples of the 
droplet size distributions (1025, 1220 and 1435 
JST) are shown in Fig. 8. 

Simultaneous samplings of cloud- and rain- 
drops were made from 1025 to 1530 JST at 
five minute intervals. Fig. 9 shows the varia- 
tion of the concentrations in four groups of 
droplets with radii of 4—8.5 microns, 55—90 
microns, 0.I—0.2 mm and 0.2—0.3 mm. 
The rain intensity (above 0.1 mm/hr) and the 
radii of largest raindrops are shown at the 
bottom of the figure. 

Observations of the vertical variation of 
raindrop size distribution (Oxita, 1958 b) 
have shown that, when the rate of rainfall 
exceeds 1 mm/hr, the concentration of small 
drops with radius 0.05—0.25 mm _ decreases 
rapidly from the cloud base downwards due 
to mutual coalescence. As this type of process 
would certainly also take place in the cloud 
itself, one may from Fig. 9 infer that a con- 
siderable number of the raindrops of this size 
were formed within a layer of perhaps 1,000 
m thickness above the station. As the air 
temperature at Mount Kuro on this occasion 
varied between 13.0 and 14.4° C, the production 
of the drops occurred at temperatures well 
above freezing. . 


(d) July 28 

In the early morning Mt. Kuro was covered 
by nimbostratus cloud, which gradually trans- 
formed into cumulus around 0900 JST, the 
thickness being between 1,100 and 2,000 m. 
Tellus XIII (1961), 4 
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Fig. 8. Droplet size distributions on July 27, 1956. 
Weather at Mt. Kuro: 


Air temperature was between 13.0 and 14.4° C. 
Wind direction was W but it turned to NW after 1500 
Sk 


Wind velocity was 2.8 to 5.8 m/sec. 


The size distribution of the cloud droplets at 
0900 JST is shown in Fig. 10. 
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Fig. 9. Variation of droplet concentrations on July 27, 1956. 
R. R.: Rate of rainfall 
M. D.: Maximum drop diameter 
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From early morning Mt. Kuro was covered 
by cumulus or cumulonimbus clouds. With 
rising cloud base the station came into cloud- 
free air at 1100 JST. Rain started at 1150 JST 
and continued until 1800 JST. Two examples 
of the cloud droplet size distributions (0945 
and 1505 JST) are shown in Fig. 10. 


(f) August 12 


From early morning cumulus or cumulus 
congestus clouds, formed over the west side 
of Mt. Asahi and moved in over Mt. Kuro. 
From 1235 to 1300 JST drizzle fell on the 
mountain. The thickness of the cloud esti- 
mated from photographs was about 3,000 m. 
The air temperature at the cloud top, estimated 
from Wakkanai and Sapporo radiosonde data 
was about 0°C, and the drizzle was thus 
produced in the cloud at temperatures above 
freezing. The droplet size distribution at 1235 
JST is shown in Fig. 11. Drizzle with maximum 
drop radius of 0.3 mm also fell from 1740 JST, 
when the thickness of the cumulus congestus 
cloud was about 3,000 m. 


Mt. Kuro was covered with nimbostratus 
clouds throughout the day. Light rain or 
drizzle fell continuously with a rate of 0.2 to 
0.3 mm/hr. Two examples (1155 and 1415 
JST) of the drop size distributions are shown 
in Fig. 11. 


(h) August 14 


Mt. Kuro was covered with nimbostratus 
clouds as on the preceding day, the precipita- 
tion continuing in the form of intermittent 
rain. Two examples (0455 and 1505 JST) of 
drop size distributions are shown in Fig. 11. 

A continuous sampling of both cloud- and 
rain-drops was made from 0500 to 0710 JST 
as on July 27. The result is shown in Fig. 12. 
It is of interest to note that the increase in 
concentration of large cloud droplets as well 
as of small raindrops occurred concurrently 
with increasing rainfall intensity (e.g. 0530 to 
0545 JST and after 0640 JST). These large 
cloud droplets or small raindrops are supposed 
to be formed at the lower part of the cloud 
under a rain generating cell. 
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Droplet Radius (A) Fig. 11. Droplet size distributions on August 12, 13 and 14. 
Weather at Mt. Kuro was as follows: 

Fig. 10. Droplet size distributions on July 28 and 29, 
1956. 


Date Aug. 14 


Weather at Mt. Kuro: 


Date | July 23 | July 29 


Air temperature °C 
Wind direction 
Wind velocity m/sec 


13.0—15.5 | 12.0— 14.5 
WSW—W | W—WNW 
3-3—5-9 1.9—4.3 
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| Aug. 12 | Aug. 13 


Air tempera- 


ture °C 15.0—20.4]14.2—17.0 | 11.5—17.5 
Wind direc- SW— SW— 

tion W—SW WSW WSW 
Wind veloc- 

ity m/sec 1.7—3.5 | I.I—3.0 | 0.7—5.2 


Concentration of Droplets cm” 


4-4 420 rain» 


054 084 025 10 07, 012 U8 


RR. nz 029 1.8 1.2 z 2 Ben 
Gan OM Za Ge aL OOS. ON], À. 
MD mm 6 LE ft fe 42 “Ub 2226 


Fig. 12. Variation of droplet concentrations on August 
14, 1956. Symbols are the same as in Fig. 9. 


5. Measurements in August of 1955, 1957 and 
of 1958 


(a) August 17 and 19, 1955 


Simultaneous samplings of cloud- and rain- 
drops were made on several occasions on Mt. 
Kuro on August 17 and 19, 1955. Some of the 
cloud- and rain-drop size distributions are 
presented in Fig. 13 (cf. Oxrra, 1958 b). 


(b) August 23, 1957 


Light showers from cumulonimbus clouds 
fell from midnight of August 22. Samplings 
were made from 0600 to 0900 JST. Four size 
distribution curves are presented in Fig. 14. 


(c) August 19, 1958 

Fig. 15 shows two examples of cloud 
droplet size distributions on August 19, 1958, 
when measurements of cloud liquid water 
content also were made (see Paragraph 3). 
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Fig. 13. Droplet size distributions on August 17 and 19, 


195$. 
Weather at Mt. Kuro: 


Date | Aug. 17 | Aug. 19 
Air temperature °C — 12.6 
WNW— 
Wind direction NW WSW—W 
Wind velocity m/sec 1.3—4.7 5.5—9.9 
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1957 Aug. 23 
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x07roor 
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Fig. 14. Droplet size distributions on August 23, 1957- 


6. Some characteristics of drop size distribu- 
tion curves 


The cloud- and rain-drop size distributions 
represented by Figs. 7, 8, 10, II, 13, 14 and 15 
show several characteristic features. In general 
the distribution curves for each day have 
similar forms. The size distributions on August 
17 and 19, 1955 and on August 23, 1957 
are represented by almost straight lines in the 
log N-log r graph, where N is the concentra- 
tion calculated for 5 micron radius intervals 
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Fig. 15. Droplet size distributions on August 19, 1958. 


Weather at Mt. Kuro: 
Air temperature was 12.0° C. 
Wind direction was S. 
Wind velocity was 2.3 m/sec. 


and r the droplet radius. These distributions 
may therefore be represented by the relation 


NaN, he (3) 


The values of N, and « of some representative 
size distributions are listed in Table 2. It is of 
interest that on August 23, 1957 « had constant 
values of about 5.1 up till about 0800 JST 
but then suddenly changed to a value around 
6.8. 

‘The distribution curves for the other days 
deviate somewhat from a straight line, especi- 
ally the curves for July 25 and 26, 1955 and for 
August 19, 1958, show a sharp bend. It has 
been assumed, however, that these distributions 
should also be expressible by the equation (3) 


and the corresponding approximate values of 
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Date 


Aug. 17, 


Aug. I9, 


July 25, 


July 26, 
July 27, 


July 28, 
July 29, 


Aug. 13, 


Aug. 14, 


Aug. 23, 


Aug. 19, 


1955 


1955 


1956 


1956 
1956 


1956 
1956 


1956 
1956 


1956 


1957 


1958 


Time 


(JST) 


I145 
1230 
1431 
1522 
1218 
1225 


1255 


0835 
0845 
0930 
1005 
IOI5 
1320 
1540 
1655 
0920 


1025 
1050 
1145 
1220 
1250 
1405 


‚1435 


1510 
0900 


0910 
0945 
1025 
1155 
1505 
1605 


1235 


1155 
1415 
0455 
0605 
0840 
0930 
1025 
1335 
1505 
0637 
0700 
0820 
0840 


1610 
1612 
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w 
- 


wr 
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No 


+ 10-13 
+ 10-1? 
+ 10-13 
. 10-14 


* 10-15 
+ 10-12 
+ 10-14 


- 10716 
+ 10-14 
+ 10-17 
+ 10-14 
+ 10-1 
+ 10-21 
+ 10-27 
+ 10-26 


+ 10-13 


- 10-15 
+ 10-13 
- 10-17 
« 10-18 
+ 1016 
- 10-17 
+ 10-15 
» 10-16 


+ 1033 


. 10-20 
+ 10-29 
+ 10-24 
. 10-24 
+: To-14 
+ 10-18 


- 10-17 


+ 10-17 
+ 10-14 


. 10-18 
+ 10-17 
+ 10-15 
+ 10-15 
. 10-20 
+ 10-15 
+ 10-16 


+ 10-14 
+ To-14 
+ 10-19 
+ 1018 


+ 10-18 
* 10-14 


II.O 


6.6 
Io.I 
8.0 
7.6 
4.9 
6.2 
5.5 
5.8 
4.9 
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Table 2. 
Number of droplets m= 
Range M elo ud 
of » R mm/hr Radius (u) en 
(¥) 50—125 125—250|250—400 
8.5—600| 1.20 | 2.6 48,000] 15,400 730 Ns 
3—750| 0.83 | 3.0 57,000| 15,300 630 » 
3—750| 1.03 | 1.0 168,000| 4,500 180 » 
3—750| 2.56 | 0.46 * 49,000} 1,300 150 » 
3—800 1975301027 27,000| 21,000 5,100 Cb 
3—850| 0.85 | 1.6 14,200| 3,200 280 » 
8.5—700| 4.5 0.68 280,000| 11,700 430 » 
3—175| 0.14 | — 350 ign — Cu 
3—175| 0.12 | — 1,600 60 — » 
3—125| 0.18 | — 520 — — » 
3—175| 0.025] — 320 18 — » 
3— 50] 0.18 | — = == = » 
3— 17| 0.021] — — = == » 
8.5— 30] 0.15 | — — — — » 
Sr] OD Dee a = = , 
3—750| 0.016} — 150 60 37 Ns 
3500| 0.43 | 0.014 77° 35 4-9 » 
3—700| 0.89 | 2.0 3,400 270 13 » 
3—300| 1.7 [8.5 - 1073 1,560 28 — » 
3—500| 0.87 | 2.6 11,700 430 8.1 » 
3—400| 0.94 | 0.46 9,400 440 II » 
8.5—600| 0.79 |2.8 - 103 2,600 2 0.25 » 
8.5—600| 1.5 0.062 2,600 124 18 » 
3—350| I.o 0.45 8,400 900 7.5 » 
3— 25| 0.056, — — — — Cu 
3— 50| 0.022] — == = = Cb 
3—100| 0.95 = 150 — = » 
3— 30] 0.2 - u — - » 
She SONT ZI = = ve 4 
3—175| 1.2 = 4,600 280 —- » 
Sam) ON | 650) — = » 
3—150| 0.016] — 60 3.8 — » 
3—250| 0.043/5.0 » 1074 720 5-4 —— Ns 
3—250| 0.047| 0.28 5,000 136 — » 
3—800| 0.067] 0.31 870 49 2.5 » 
3—300| 0.048]4.4 + 1074 107 25 — » 
10—500| 0.058| 0.23 1,130 450 20 » 
3—400| 0.081] 0.016 1,460 7 5.2 » 
3—175| 0.043] 0.12 4,600 80 — » 
3—125| 0.048] — 480 — — » 
3—300| 0.13 | 0.015 2,400 113 — » 
3—250| 0.86 — 47,000 410 — Cb 
3—300| 1.34 — 38,000] 1,320 44 » 
3—125| 1.57 — 1,700 aa — » 
3—100| 1.17 — 2,200 2 == » 
302170 | S358] Se 17,000 105 5 Ns 
3—225| 0.46 — 40,000 1,220 — » 


N, and « are also listed in Table 2. The cloud fall intensities by using the sooted paper 
liquid water contents in Table 2 are calculated 
from the size distribution curves and the rain- 


method. 
The effects of cloud form and rain on the 
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Table 3. 
À 


Maximum drop radius (1) Cloud 
< 50 50—250 250—500 | 500 > io 
MCA mE allen. eee, ern ae - 5.4 5.7 4.8 
RETRO as, Aa == 4.6—6.8 4.9—7.3 4.I—6.2 Ns 
Number ot Samples. 0... o 10 14 10 
VIG UMNe CSOD Te ne mo 8.4 6.4 se 5.0 
REISEN AR SE RR 6.6— 10.0 4.9—I0.I 5.2 4.5 —5.3 Cb 
Number Of Samples... ..... 13 8 I 3 
IMI Samara Or Ge ee. 8.1 5.8 — — 
Ranee ofen. seele ne series 4.6—13.1 4.7—6.7 — — Cu 
Numberonr samplesijo- 6. 4... <0. - 16 7 o o 


value of a are demonstrated in Table 3. In 
this table 82 size distribution samples including 
those shown in Table 2 are classified according 
to cloud form and maximum drop radius in 
the distribution. The table indicates that the 
mean value of xis smallest with rain of moder- 
ate intensity, in which raindrops with radius 
larger than 0.5 mm exist. This is because the 
concentrations of large cloud droplets or small 
raindrops are large (August 17 and 19, 1955; 
July 27, 1956) or because the clouds have 
relatively few small droplets (July 26, 1956). 
When drizzle or light rain falls inside the cloud 
the value of « is between 5.2 and 6.4. For 
non-precipitating clouds its value is largest 
(8.1—8.4). There is of course a large varia- 
bility of « between individual distributions and 
there seems to be no distinct difference in the 
a-values for various cloud forms. 

Large cloud droplets with radii above 50 
microns are always found in nimbostratus 
clouds (August 17, 1955; July 26 and 27, 
August 13 and 14, 1956; August 19, 1958). In 
cumulonimbus or cumulus congestus about 
half of the size distributions show such large 
droplets (August 19, 1955; July 29 and August 
12, 1956; August 23, 1957). These were also 
found in cumulus cloud when the top had 
been seeded from higher clouds (July 25, 
1956). On the other hand, in cumulus clouds 
of thickness less than 2,000 m and with no 
surrounding clouds no droplets of radius 
above 25 microns were found (July 28, 1956). 

The observations on July 27 and on August 
12 and 14 in 1956 further suggest that the large 
droplets are formed in the part of the cloud 
where the air temperature is above freezing. 

Brown and BRAHAM (1959) made aero- 
plane observations of large droplets in tropical 
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cumuli. The curves A and C in Fig. 16 show 
their narrow and broad distributions respec- 
tively. The size distribution B—R in the 
figure is the one obtained by BATTAN and 
REITAN (1957) in tropical cumuli. No deter- 
mination was made of droplets with radii 
between so and 100 microns. If it is assumed 
that the size dristibution of these droplets be 
represented by the broken line in the figure. 
The corresponding values of « in tropical 
cumuli are 6.7 and §.0 for A and C distributions 
respectively. These values are in good agree- 
ment with those given in Table 3 for cases 
when the largest droplet radius is between so 
and 250 microns. In Fig. 16 the size distri- 
butions in fair weather cumulus and cumulo- 
nimbus clouds obtained by WEICKMANN and 
AurMm KAMPE (1953) are also shown. The 
corresponding values of « are about 5 and 3 for 
fair weather cumuli and cumulonimbus clouds 
respectively. These values are much smaller 
than the values of « given in Table 3 in the 
corresponding cases. 


7. Concentration of large cloud droplets 


The concentrations of large cloud droplets 
with radii above so microns are shown in 
columns 8—10 of Table 2. There is great varia- 
bility in the concentration and no significant 
difference seems to exist between cumulonim- 
bus and nimbostratus clouds. In the cumulus 
clouds of July 25, 1956 the concentration of 
droplets of radius so to 125 micron varies 
between 320 and 1,200 m®, 

The concentrations of large droplets in 
tropical cumuli and in cumulus and layer 
clouds over England are given in Table 4 
(BLANCHARD, 1957; MURGATROYD and Gar- 
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Table 4. Concentrations of large droplets in tropical cumuli and in cloud over England. 


Investigator Concentration m? Size intervalmm . . Type of Cloud 
— 

Bowen 17,150 0.2—1.6 Tropical cumuli 
Brown and Braham 1,492 (mean) 0.15—0.65 ~ Tropical cumuli 

7,581 (max.) € 
Murgatroyd and Garrod 91,600 (max.) = ON Cold cumuli, England 

6,690 (max.) = OT . „| Warm cumuli 

16,800 (max.) >o.1 Mixed Cu and Sc 

Singleton 40—18,000 O.1—0.25 Layer clouds, England 
Blanchard 59,000—83,000 0.1—0.34 | Jawaiian cumuli 


Concentration of Droplets cm? (sp)! 


10,9 70° 10 
Droplet Radius (pr) 


ROD, 1960; SINGLETON, 1960). The concentra- 
tions found in Table 2 are generally of the 
same order of magnitude. On „four days 
(August 17 and 19, 1955; August 23, 1957; 
August 19, 1958) the concefitrations exceeded 
40,000 m~® and on August 1g,and 19 in 1955 
concentrations above 100,000:m~*. were found 
in precipitating clouds. Blanchard also found 
numerous large droplets in Hawaiian, warm 
cumuli, and according to the experience of the 
author such high concentrations ef large 
droplets frequently occur in heavily or moder- 
ately precipitating clouds. * 
€ 

8. Conclusion 

The observations of cloud- and rain-drop 
size distributions reveal that even in a tem- 
perate latitude a considerable number of large 
droplets of radii above so microns are con- 
tained in nimbostratus and cumulonimbus and 
even in cumulus clouds. Some observational 
results further suggest that these large dropicts 
are formed in the warm part of the cloud 
where the air temperature is above freezing. 

These droplets would be formed by con- 
densation-coalescence process. WELANDER (1959) 
calculated the effect of coalescence on the drop- 
let size distribution and deduced that both in 
the “time-dependent” and in the “steady 
state” cases the spectra could be approximated 


Fig. 16. Droplet size distributions. 
W—K Cb: WEICKMANN and AUFM KAMPE, cumulonim- 
bus cloud 
W—K Cu: WEICKMANN and AurM KAMPE, cumulus 
cloud 
B—R: BATTAN and REITAN 
B—B: Brown and BRAHAM 
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by a power law, r® and r? respectively. He 
further stated that the transient model leading 
to the r-5-law should apply best to “young” 
clouds such as fair weather cumuli while the 
steady state model should hold for “older” 
clouds such as cumulonimbus or cumulus 
congestus. The theory was confirmed by 
observations made by Weickman and Aufm 
Kampe. However, our observations indicate 
that the value of the parameter « lies between 
5.2 and 6.4 for “. Ider” clouds (nimbostratus 
and cumülonimbus clouds) and between 8.1 


and 8.4 for “young” clouds (cumuli). The 
discrepancy may be due to the presence of gi- 
antnuclei or to other factors. 
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SHORTER CONTRIBUTION 


Tropospheric Geopotential and Velocity. Variances 


By C. EUGENE BUELL, Kaman Nuclear, Colorado Springs, Colorado 
(Manuscript received May 29, 1961) 
In connection with a recent publication of relationship between (sin y) o, (timewise vec- 


Gopson and MACFARLANE (1958), we would tor standard deviation of wind at a point) and 
like to present Figs. ı and 2 which show the S, (timewise standard deviation of height at a 
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point) over North America. The data were 
obtained from the computations made at the 
University of Wisconsin (LAHEY ET AL., 1958; 
1959). The point symbols indicate the latitude 
and season to which the data pertain. It seems 
quite evident that the relation proposed, valid 
for the northern latitudes (larger values of S,), 
does not describe the situation in general. In 
particular, the vertical tangent at S, = o and 
the sudden increase in values of (sin @) o4 for 
small S, are absent. A general description was 
attempted in earlier analyses (BUELL, 1957) 
where real winds were used in place of geo- 
strophic winds and where winds and height 
(pressure) data were by no means synchronous. 
(This heterogeneity of the data naturally in- 
creases the depression.) 

That such a general description is possible is 


523 
connected with the fact that, though height 


systems vary in size from place to place, this 
variation amounts to only about 20 % of a 
nominal value and may be absorbed in the 
dispersion about the regression line. 

The relation suggested by the Arctic data 
(Gopson and MACFARLANE, 1960) tends to 
imply that the size of height systems decreases 
with latitude, which is not generally the case 
(BUELL, 1954; 1961). This would suggest that 
the relation concerned (Gopson and Mac- 
FARLANE, 1960) is restricted to (near) the area 
of the initial data. 

Data presented herein was developed in con- 
nection with Contract No. AF19(604)-7282, 
Geophysical Research Directorate, Meteoro- 
logical Development Laboratory. 
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